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PREFACE 


Our initial motivation was to provide an up to date translation of the monograph 
m written in french by the hrst author, taking account of more recent developments 
of inhnite dimensional dynamics based on the Lojasiewicz gradient inequality. 

While preparing the project it appeared that it would not be easy to cover the entire 
scope of the french version in a reasonable amount of time, due to the fact that the non- 
autonomous systems require sophisticated tools which underwent major improvement 
during the last decade. 

In order to keep the present work within modest size bounds and to make it available 
to the readers without too much delay, we decided to make a hrst volume entirely 
dedicated to the so-called convergence problem for autonomous systems of dissipative 
type. We hope that this volume will help the interested reader to make the connection 
between the rather simple background developed in the french monograph and the 
rather technical specialized literature on the convergence problem which grew up rather 
fast in the recent years. 
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Chapter 1 


Introduction and basic tools 


1.1 Introduction 

The present text is devoted to a rather specific subject; convergence to equilibrium, 
as t tends to infinity, of the solutions to differential equations on the positive halfline 
{t > 0} of the general form 

U'{t)+AU(t) = 0 

where ^ is a nonlinear, time independent, possibly unbounded operator on some Ba¬ 
nach space X. By equilibrium we mean a solution of the so-called stationary problem 

AU = 0. 

By the equation, taken at a formal level for the moment, it is clear that if a solution tends 
to an equilibrium and if A is continuous : X ^ Y for some Banach space Y having 
X as a topologically imbedded subspace, the ’’velocity” U'{t) tends to 0 in Y. If the 
trajectory U is precompact in X, it will follow that this means some strong asymptotic 
flatness of U (t) for t large. Conversely, systems having this property do not necessarily 
enjoy the convergence property since trajectories might oscillate (slower and slower at 
infinity) between several stationary solutions. 

A well known convenient way to study the asymptotic behavior of solutions is to 
associate to the differential equation a semi-group S{t) of (nonlinear) operators on 
some closed subset Z of the Banach space X , defined as follows: for each t > 0 and 
each z G X for which the initial value problem is well-posed, S{t)z is the value at t 
of the solution with initial value 2 . Since the initial value problem does not need to be 
well-posed for every z G X An general Z will just be some closed set containing the 
trajectory 

r(z) = IJ 5(f)/ 

t>0 

For some results the consideration of r( 2 ;) will be enough, for some others (for instance 
stability properties) it will be preferable to take Z as large as possible. The standard 
terminology used in the Literature for such semi-groups is “Dynamical systems” and 
we shall adopt it. Since the operator A does not depend on time, both equation and 
dynamical system are called autonomous. According to the context, the word ’’trajec¬ 
tory” will mean either a solution of the equation u{t -f s) = S{t)u{s) on the halfline. 
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or the closure of its range. 

The present work concerns dissipative autonomous systems. In the Literature the 
term “dissipative” has been used in many different contexts. Here, dissipative refers to 
the existence of a scalar function <!> of the solution U which is dissipated by the system, 
in the sense that it is nonincreasing; 

Vs >0, Vf > s, <i>([/(f)) < <l>(C/(s)). 

If in addition $ is coercive , this implies that U{t) is bounded in X. The problem of 
asymptotic behavior becomes therefore natural. Such non-increasing functions of the 
solution play an important role in the theory of stability initiated by Liapunov. For this 
reason, in this text, they will be called Liapunov functions (resp. Liapunov functionals 
if X is a function space). 

Let us now define more precisely the main theme of the present text. The structure 
of trajectories to dynamical systems tends to become more and more complicated as 
the dimension of the ambient space X increases. When X = M, ^ is just a scalar func¬ 
tion of the scalar variable U and if A is locally Lipschitz, as a consequence of local 
uniqueness, no trajectory other than a stationary solution can cross the set of equilibria. 
As a consequence all bounded solutions are monotonic, hence convergent. In higher 
dimensions, what remains true is that convergent trajectory have to converge to a sta¬ 
tionary solution. But the equation u” -f u = 0 , which can be represented as a first 
order differential equation in AT = exhibits oscillatory solutions, and even when a 
strictly decreasing Liapunov functions exists, two-dimensional systems can have some 
non-convergent trajectories. Our main purpose is to find sufficient conditions for con¬ 
vergence and exhibit some counterexamples showing the optimality of the convergence 
theorems. Finding sufficient conditions for convergence is a program which was initi¬ 
ated by S. Lojasiewicz when X = M.^ and A = XF with F a real valued function. By 
relying on the so-called Lojasiewicz gradient inequality, he showed that convergence 
of bounded solutions is insured whenever F is analytic. From the point of view of a 
sufficient condition expressed in terms of regularity, this result is optimal; there are 
C°° functions on AT = for which the equation U'{t) + VFU{t) = 0 has bounded 
non-convergent solutions. An explicit example was given by Palis & De Melo in ll73l . 
and in this text we extend their example in such a way that any Gevrey regularity con¬ 
dition weaker than analytic appears unsufficient for convergence. 

This text is divided in 11 chapters; the first 3 chapters contain some basic material 
useful either to set properly the convergence question, or as a technical background 
for the proofs of the main results. In Chapter 4 we fix the main general concepts or 
notation concerning dynamical systems. In chapter 5 a general asymptotic stability 
criterion is given, generalizing the well known Liapunov stability theorem (Liapunov’s 
first method) in a framework applicable to infinite dimensional dynamical systems and 
in the same vein, a finite-dimensional method used by R. Bellman to derive instabil¬ 
ity from linearized instability is applied to some infinite dimensional dynamical sys¬ 
tems. Chapter 6 is devoted to the definition and main properties of a class of “gradient¬ 
like systems” in which the question of convergence appears fairly natural. Chapter 7 
concerns the general invariance principle and its connection with Liapunov’s second 
method. After Chapter 8, in which simple particular cases are treated by specific meth¬ 
ods, Chapter 9 and 10 are devoted to convergence theorems based on the Lojasiewicz 
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gradient inequality, respectively in finite dimensions and infinite dimensional setting 
with applications to semilinear parabolic and hyperbolic problems in bounded domains. 
Chapter 11 is devoted to a somewhat informal description of more recent or technically 
more elaborate results which are too difficult to fall within the scope of a brief mono¬ 
graph. 

We hope that this text may help the reader to build a bridge between the now clas¬ 
sical Lojasiewicz convergence theorem and the more recent results on second order 
equations and infinite dimensional systems. 


1.2 Some important lemmas 


The first lemma is classical and is recalled only for easy reference in the main text. 

Lemma 1.2.1. (Gronwall Lemma) Let T > 0, X G L^{0,T), A > 0 a.e. on {0,T) 
and C > 0. Let ip G L°°(0, T), if > 0 a.e. on (0, T), such that 

ip{t)<C-\- f X{s)ip{s)ds, a.e. on(0,T) 

Jo 


Then we have 


Proof. We set 


Lp(t) < Cexp 



tj;{t) = C + 


X{s)Lp{s)ds, 


a.e. on (0, T) 

Vf G [0, T] 


Jo 

Then ijj is absolutely continuous, hence differentiable a.e. on (0,T), and we have 
= X{t)(p{t) < X{t)tp{t) a.e. on (0,T). 


Consequently, a.e. on (0, T) we find ; 

J X{s)ds^] < 0. 

Hence by integrating 

'ip{t) < Cexp(^ J X{s)ds^ , Vf G [0, T]. 

The result follows, since ip < ij: a.e. on (0, T) □ 

The next lemmas will be useful in the study of convergence and decay rates 

Lemma 1.2.2. (cf. e.g. US 31 / . ) Let X be a Banach space, to € K and z G 
C{{to,oo); X). Assume that the following conditions are satisfied 

z G L^((to,oo);X) 

Z is uniformly continuous on [to, oo)with values in X. 


lim \\z{t)\\x = 0 

t —¥00 


Then 


( 1 . 1 ) 

( 1 . 2 ) 
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Proof. Let e > 0 be arbirary and let J > 0 be such that 

sup \\z{t + h- z{t))\\x <e 

t^[to ,oo) ,h^[0,6] 


Then we find easily 


1 

Vte[io,oo), \\z(t)\\x<e+-J ||z(s)||xds. 


implying 

limsup ||2:(i)||x < e 

t—¥oo 


The conclusion follows immediately □ 

Lemma 1.2.3. Let X be a Banach space, fo G ® tmd u € oo); X) . Assume 

that there exists H G ( 7 ^ ((fo, c»), R), p G (0,1) and c > 0 such that 


H{t) >0 for all t>to. 

— H'(t) > cH{ty~^ ||u'(<)||x for all t > to- 
Then there exists p £ X such that lim uii) = ld in X. 

t^oo 


Proof. By using (ll.41 i. we get for all f > fo 


> cr]\\u{t)\\x- 


(1.3) 

(1.4) 


(1.5) 


By integrating this last inequality over {to, T), we obtain 


ho crj 


( 1 . 6 ) 


This implies u' G oo); X). By Cauchy’s criterion, lim u(f) exists in X. □ 

t—foo 

Lemma 1.2.4. Let T ^ 0, let p be Q. nonnGgQtivG sc^uqkc iHlGgro-blG function on [0,2^)- 
Assume that there exists two constants 7 > 0 and a > 0 such that 

VfG[0,T], p'^{s)ds < ae~'^*. 

Then setting b := — l),/or all 0 < t < T < T we have: 

J{t,T) := p{s)ds < s/abe~^. 


Proof. Assume first that t — t < 1. Then we have 
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If T — t > 1 we reason as follows. Let N be the integer part of t — t, we get 

N-l ft+i+l 

p(s)ds + 

It+N 


N—1 /‘T 

J{t,T) < E p{s) ds+ p(s) ds 

-i~f) Jt-\-N 


2=0 

N-1 


< 


E _ Nj-t+i) j— -yit-^N) 

y'ae 2 + ^/ae 2 


i=0 


e 2 

< \/a— -e 2 . 

6 2—1 

□ 

Lemma 1.2.5. Let p be a nonnegative square integrable function on [1, c»). Assume 
that for some a > 0 and a constant K > ff, we have 


/ ‘2t 

p'^{s)ds < Kt 

Then for all t > t > 1 we have: 


-2a-l 


Pis)ds < ^ “• 


Proof. By Cauchy-Schwarz inequality, for all f > 1 we may write: 

('2,t 


Pis)ds < Vi (iLf-2a-l)l/2 ^ 


hence 


/ r poo 00 p2 ^ t ^ 

p{s) ds< p{s) ds = E . p{s) ds < Vk “ = -—t 


Vk 


□ 

Finally, in the application of the Lojasiewicz gradient inequality to convergence 
results, the following topological reduction principle will play an important role. 

Lemma 1.2.6. Let W and X be two Banach spaces. Let U G W be open and E : 
U —> R and Q : U —> X be two continuous functions. We assume that for all a £ U 
such that Q(a) = 0, there exist Oa > 0, 0(a) € (0,1) and c(a) > 0 

\\Giu)\\x > c(a)\E{u) - \/u : ||w - a||M^ < CTo. (1.7) 

Let T be a compact and connected subset ofG~^{0}. Then we have 

(1) E assumes a constant value on F. We denote by E the common value of E(a), 
a G F. 

(2) There exist ct > 0, 0 G (0,1) and c > 0 such that 

dist(u, F) < CT ==^> ||f/(M)||x > c|ii^('u) — E\^~^ 
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Proof. By continuity of E we can always assume that aa is replaced by a possibly 
smaller number so that \E{u) — E{a)\ < 1 for all u such that ||u — a||vy < (Ja. Let 
a £ r and 

K = {bGr/ E{b) = E{a)}. 

It follows from O that K is an open subset of L which is obviously closed by conti¬ 
nuity and since L is connected by hypothesis we have K — T. 

On the other hand, since L is compact, there exist oi, • • • , Op £ L such that 

rc 

i^l 

The result follows with cr = ^ inf tJo^, c = inf c{ai) and 9 = inf 6{ai). 

□ 


1.3 Semi-Fredholm operators 

Let E, F be two Banach spaces and A : E —^ L" be a linear operator. We denote by 
N{A) and R{A) the null space and the range of A, repectively. 

Definition 1.3.1. A bounded linear operator A £ L{E, F) is said to be semi-Fredholm 

if 

(1) N {A) is finite dimensional, 

(2) R{A) is closed. 

We denote by SF{E, F) the set of all semi-Fredholm operators from E to F. 

Remark 1.3.2. The fact that N{A) is finite dimensional implies that there exists a 
closed subspace X of E such that E = N{A) 0 X (cf Il20l p. 38). Moreover R{A) = 
A{X) is a Banach space when equipped with the norm || • ||i?. 

Theorem 1.3.3. Let A £ L{E, F) and assume that N{ A) is finite dimensional. Then 
we have A £ SF{E, F) if and only if 

> 0, Vu £ X \\Au\\f>p\\u\\e- (L8) 

Proof. (IL81 l implies that R{A) is closed. In fact, let (/„) = (Aun) be such that 
/„ —/ in F. Let (xn) and (j/„) be such that = a;„ + with (a;„) C X and 
iVn) C N{A). So /„ = Axn- Then the inequality \\xn - XmWs < ^\\fn - fmWp 
implies that (xn) is a Cauchy sequence, hence converges. Let x be the limit. We have 
AXn —Ax so / = Ax. 

Conversely, R{A) is a Banach space and C ■.= A/x ■ X —^ R{A) is bijective and 
continuous, by Banach’s theorem we get that C~^ is continuous and (IL81 l follows. □ 

Remark 1.3.4. \f A : E —F is a topological isomorphism, then A £ SF{E, F) 
with N{A) = {0}. Conversely, as a consequence of Banach’s theorem, if A £ 
SF{E, F) with N{A) = {0} , then A : E —> RiA) is a topological isomorphism. 

Theorem 1.3.5. Let A £ SF{E,F) and G £ L{E,F). Assume that G is compact, 
then A + G€ SF{E, F). 
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Proof. We divide the proof into 3 steps : 

Step 1 : If (un) C E with ||u„|| < 1 and {A + G)('u„) —0, then (m„) has a 
strongly convergent subsequence in E. Indeed we can assume Gun —> g G F. Let 
Un = Xn + Un, Xn G X, ijn G N{A) where X is as in the remark fl. 3. 21 Since 
Aun = Axn —> —g, (xn) is Convergent in E. Then (?/„) is bounded in N(A), since 
dimiV(A) < c» we can assume that yn — y'm E with y G N{A). In particular 
Un= XnE yn IS Convergent in E. 

Step 2 : Let (m„) C N{A + G) with ||wri|| < 1. By step 1, (u„) is precompact in E, 
hence the unit ball of iV(A+ G) is precompact and consequently dim A^(A+ G) < oo. 
Step 3 : Let F be a Banach space such that E = N{A + G)^Y. Assuming 
R{A + G) not closed, then by Theoreii lL3.3l we can hnd yn G Y with ||y„|| = 1 
and {A + G)y„ —0. By step 1, up to a subsequence we can deduce yn —>■ y in 
E. We immediately find ||y||_E and y G Y. Hence since {A + G)y„ —0 we have 
y G N{A + G). Since 7V(A + G) fl F = {0}, we end up with a contradiction since 
y G N{A + G) n F and ||y||_E = 1- □ 

For the next corollary, we consider two real Hilbert spaces V, H where V C H 
with continuous and dense imbedding and H', the topological dual of H is identified 
with H, therefore 

V C H = H' cV 
with continuous and dense imbeddings. 

Corollary 1.3.6. Let A G SF(V, V') and assume that A is symmetric. Then A -\- P : 
V —> V' is an isomorphism where P : V —> X(A) is the projection in the sense of 

H. 

Proof. First we have iV(A+P) = {0}. Indeed if Au+Pu = 0, we have Am = —Pu G 
N{A), then Am G N{A) fl P(A) = {0}, so Am = 0, hence u = Pu = —Au = 0. 

On the other hand, since A G SF(V, V), dim N{A) < oo and then P is compact. By 
Theorem lL3.5l A + P G SF{V, V), then P(A + P) is closed. Now since A + P is 
symmetric and A^(A + P) = {0} then P(A + P) is dense in V, hence P(A + P) = V. 
By Banach’s theorem we get that (A + P)~^ G L{V', V). □ 

Example 1.3.7. Let H be a bounded and regular domain of V = Hq (H) 

A = -A+p{x)I, pGL°°{n) 

G := p{x)I : V —>■ V is compact. —A G Isom(F, F') then by Theorem [03] 
A G SF{V, V). Corollary 11. 3. hl imnlies that A + P G Isom {V, V). 

I. 4 Analytic maps 

In this section, we introduce a general notion of real analyticity valid in the Banach 
space framework which will be essential for the proper formulation of many conver¬ 
gence results applicable to P.D.E. One of the difficulties we encounter here is that the 
good properties of complex analyticity cannot be used and all the proofs have to be 
done in the real analytic framework. For example, in this framework the result on com¬ 
position of analytic maps is not so trivial as in the complex framework and its proof 
is generally skipped even in the best reference books. Here we shall give a complete 
argument relying on the majorant series technique of Weirstrass. 
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1.4.1 Definitions and general properties 

Definition 1.4.1. Let X, Y be two real Banach space and a G X. Let U be an open 
neighborhood of a in X. A map f : U —> Y is called analytic at a if there exists 
r > 0 and a sequence of n—linear, continuous, symmetric maps (M„)nen fulfilling the 
following conditions 

(1) ^ \\Mn\\cxx,Y)r‘^ < oo where 

nSN 

\\Mn\\c„(x,Y) = sup{||M„(a:i,a; 2 , • • • ,Xn)\\Y, sup \\xi\\x < !}• 


(2) B{a,r) C U. 

(3) yhGB{Q,r), /(a + /i) = /(a) + ^M„(/i(”))w/zere/i(”) = (/i,--- ,h). 

n>l 

n times 

Remark 1.4.2. Under the previous definition, it is not difficult to check that 

- V& G B{a, r), f is analytic at b. 

- f G C^{B{a, r), Y) with /^"/(a) = n!M„. 

- A finite linear combination of analytic maps at a is again analytic at a. 

Definition 1.4.3. / is analytic on the open set U if f is analytic at every point ofU. 

Exampie 1.4.4. It is clear from the definitions that any bounded linear operator, any 
continuous quadratic form and more generally any finite linear combination of restric¬ 
tions to the diagonal of continuous Ic-multilinear maps: X^ -G Y (usually called a 
polynomial map) is analytic on the whole space X. 

Proposition 1.4.5. Let f G C^{U,Y). The following properties are equivalent 

(1) f : U —> Y is analytic; 

(2) Df : U —> £(X,Y) is analytic . 

Moreover if 

n>l 

is the expansion of f{a + h) for all h in the closed ball B{0, r) G U — a, then 
Df{a + h) = Mi + J2 •) 

n>2 

is the expansion of Df(a + h) for all h in the open ball B(0, r). 

Proof. First let us explain the meaning of the formula for the derivative. It involves an 
infinite sum of expressions of the form 
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Indeed, since Df{a + h) is for all vectors h an element of C{X, Y), the formula really 
means 

Ve G X, Df{a + hm = Mi(e) + ^ 

n>2 

and for any n > 2 fixed we must identify nMn(h^^~^\ •) as the trace on the diagonal 
of of an n — 1-linear symmetric continuous map with values in C{X^ Y). The 

corresponding map is just 

Xn—l (^1 j ■ 5 ^n— 1) (0 — {x\ , • ■ •, Xji— \ , . 

Assuming 1), Let us consider a and r > 0 with B{0,r) C U — a. The expression of 
the norms of in the space of n — 1- linear symmetric continuous map with values 
in £{X, Y) shows that the formal series given by 

G X, Df{a + hm = Mm + ^ 

n>2 


satisfies \\Kn\\cn{x,cix,Y)Y'^ < oo for any r' G (0,r). The summation formula 

n^N 

for the derivative is now obvious when the expansion is finite. The general case is more 
delicate and is in fact related to the formula permitting to recover / from the knowledge 
of Df. This formula: 

f{a + h) = f{a)+ f Df{a + sh){h)ds 

Jo 

is classical and valid for any function /. When we substitute the expansion of Df in 
this formula, the summability of its terms transfers easily to yield the desired expansion 
for /. We skip the details which are classical for this part of the argument. □ 

1.4.2 Composition of analytic maps 

Let Z be a Banach space, V be an open neighborhood of /(a) and g : V —Z be 
analytic at f(a). This means that for some p > 0, we have 

g(f(a) + k) = gif (a)) + ^ Pmik*''^'^) 

m>l 


whenever ||fc||F < P and E \\Pm\\Cm{X,Z)P^ < OO- 

m^N 

Theorem 1.4.6. The map go f is analytic at a with values in Z. More precisely, setting 

m<d EjLi Pf'j—d 

(the sum is finite for any d) we have 

\\Rd\\cd{X,Z)0''^ < oo (1.9) 

d>l 


as soon as 
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go f(a + h) = go f(a) + ^ Vh, \\h\\x < a. 

d>l 

Proof. We have the obvious estimate : 

\\Rd\\c4X,Z) < II Pm\\c,r,{Y,Z) E ll^«J|---||M„^|| 

m<d \^^\—d 

where (ni,*** ,n^), |/x| = ni +- Vum and = \\Mn,\\c^.{^x,Y)- Indeed 

Rd ^ ^ ^ ^ Pm {h\ ’ ‘ ‘ ) ; * * ‘ 7 ^^Tlm -—l + l 5***7 ) 

m<d |/J.|=d 

Therefore 

Eii^^ii^.i^.^)'^' ^ EEii^-iiEii^^".ii---ii^"™ii'^' 

d>l l<m <d 

= E E E iiJ^-iiiiM„jk"^ • • • iiM„„iia"- 
= Ell^-ll E l|M„J|a"^...||M„^||a"- 

m d>m, \ii.\=d 

< Eii^-ii(Eii^^"ii^i"- 

m 

Then (11.9b follows. Concerning the convergence of the series to g o f, we notice that 

ig ° /)(a + h)-{go f){a) = ^ Pm((/(a + M “ /(a))^"*^) 

m>l 

Hence 

M 

II (5 o f){a + h)-{go f){a) - ^ Pm{{f{a + h) - /(a))('"))||z 

m=l 

< E ii^mii(Eii^^-ii'^")™ 

m>M+l 

< e for M > M{e). 

Then for M > 1 fixed 

M M 

^ P^iifia + h)- /(a))M = E E 

m=l ci>l m=l 

with Q^iihr) = • • • , 

M MM 

II P^((/(a + h)- f{a))^^^ -EE 

m=l d—1 m—1 

M 

< ^ ^ ||( 5 ^((ft.)^‘^^)|| —>■ 0 as M —>• 00. 

m=l |;i|=d>M+l 
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Finally 


M M 


Ugof)ia + h) - {gof){a) - XI X! X! < 2^ 


d—1 m—1 n\—, 


d—1 m—1 ^\—d 


for M large. But 


MM M 


E E E = 



M 


since ^ ^ = Erf for all d<M. 


since 


□ 


m=l /x|=d 


1.4.3 Nemytskii type operators on a Banach algebra 

Let yt be a real Banach algebra and / be a real analytic function in a neighborhood of 
0, which means that for some open subset U of R containing 0 we have / G C°°{U, R) 
and for some positive constants M, K 


VneN, |/(")(0)| < ME"n! 


It is clear that for any n G N the map u —u" is the restriction to the diagonal of 
of the continuous n-linear map 


n 



It follow that the map 



is analytic in the open ball Bq = B{0, A) in the sense of Subsection 1 1.4. II This map 
will be called the Nemytskii type operator associated to / on the Banach algebra A. 

Example 1.4.7. Let us consider the special case A = L°° (S) where S is any positively 
measured space. Then for any / as above the operator defined by 



for all G B{0,^) C L°°{S) and almost everywhere in S is usually called the Ne¬ 
mytskii operator on L°°{S) associated to / and is an analytic map in a ball centered at 
0. The same holds true if we replace L°°{S) by the set of continous bounded functions 
on a topological space Z or more generally any Banach sub-algebra of it. 

Remark 1.4.8. (i) The Nemytskii operator N/(u)(s) = f{u{s)) makes sense in other 
contexts, for instance from a Lebesgue space into another assuming some growth re¬ 
strictions of the generating function /. 

(ii) We shall use this operator exclusively in the case where / is in fact an entire func¬ 
tion, i.e. K can be taken arbitrarily small. 

(iii) Moreover, in the applications we shall usually need some growth restrictions on / 
or even its first derivative. 

(iv) In our applications to convergence, N/(u)(s) = f{u{s)) will usually appear as the 
derivative of a potential function G{u) = fg F(u(s)ds where F is a primitive of /. 
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I. 4.4 Inverting analytic maps 

Let X, Y be two real Banach space and a G X. Let U be an open neighborhood of a 
in X and / G C^{U, Y). The well known inverse map theorem says that if Df{a) G 
Isom {X, Y), there exists a possibly smaller neigborhood W of ainX such that f{W) 
is open in Y and f : W ^ fO^) is a C^-diffeomorphism. Moreover we have the 
formula 

yy G fiW), Dif-^){y) = [Dfir\y))]-^ 

We note that in order for / to be a diffeomorphism, we need the existence of a linear 
topological isomorphism between X and Y, namely L = Df{a), so that diffeomor- 
phisms can be reduced to the case X = Y hy replacing the general function / by the 
’’operator" g = L~^ o f. By combining (1 1.4. 5b with the fact that the map T T~^ is 
analytic on the open set Isom (X, X) C C{X, X), it is easy to prove the following 

Theorem 1.4.9. Giving a function f G C^{U,Y) which is analytic at a G U, if 
Df[a) G Isom (Jf, Y), the inverse map f~^ is analytic at f{a). 

Proof. By construction, g -.V ^ X k analytic with V an open ball of X contained in 
U and centered at a, so that we may assume V = U. As a consequence of Proposition 

II. 4.51 Dg is analytic : V —>■ C{X) and we have Dg{a) = Idc{X). Then Dg~^{x) = 

{Dg)~^og~^[x) throughout g{V), so that Dg~^ appears as a composition of 3 analytic 
maps by reducing if necessary L to a small ball around a in which Dg is sufficiently 
close to Idc{X) in the norm of C{X) to use the formula (/ — t)“^ = X)where 
T{y) = Idc{X) — Dg{y) . Finally by using once more Proposition ! 1.4. 51 the gradient 
Dg~^ is lifted to g~^ which is therefore also analytic . The details are essentially 
classical and left to the reader. □ 








Chapter 2 


Background results on 
Evolution Equations 

2.1 Elements of functional analysis. Examples of un¬ 
bounded operators 

Throughout this paragraph, X denotes a real Banach space. The norm of X is denoted 
by II II. The results will generally be stated without proof. For the proofs we refer to 
the classical literature on functional analysis, cf. e.g. Il20ll82l 

2.1.1 Unbounded Operators on X 

Definition 2.1.1. A linear operator on X is a pair (D, A), where D is a linear sub¬ 
space of X, and A : D ^ X is a linear mapping. We say that A is bounded if ||>lu|| 
remains bounded for u £ {x £ D, ||a;|| < 1}. Otherwise, A is called unbounded. 

Remark 2.1.2. \f A is bounded, then A is the restriction to D of some operator A G 
L(y, X), where F is a closed linear subspace of X containing D. On the other hand 
if A is unbounded, then there exists no operator A G L{Y,X) with Y a closed linear 
subspace of X and D CY such that A\D — A. 

Definition 2.1.3. If {D, A) is a linear operator on X, the graph of A and the range of 
A are the linear subspaces G{A) and R{A) of X defined by 

GiA) = {{uj) G X X X,uG D,f = Au} and R{A) = A{D). 

As it is usual, we shall frequently call the pair {D, A) as ”A with R>(A) = D ”. 
However one must always keep in mind that when we define a linear operator, it is 
absolutely crucial to specify the domain. 

Definition 2.1.4. A linear operator Aon X is called dissipative if we have 
VuGi:>(A),VA>0,||u-AAM|| > ||u||. 

A is called m-dissipative if A is dissipative and for all A > 0, the operator I — XA is 
onto, i.e 

V/ G X,3u G D{A),u — XAu = f. 
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Proposition 2.1.5. Let A be a linear dissipative operator on X. Then the following 
properties are equivalent. 

(i) A is m-dissipative on X. 

(ii) There exists Aq > 0 such that for each f G X, there exists u G D{A) with : 
u — XqAu = f. 


2.1.2 Case where X is a Hilbert space 

Let us denote by (•, •) the inner product of X. If A is a linear densely defined operator 
on X, the formula 

G{A*) = {{v,g)GXxX, V(u, /) G G{A), {g, u) = {v, /)} 

defines a linear operator A* (the adjoint of A), with domain 

D{A*) = {vGX,3C < 00 , \{Au,v)\ < C'||u||,Vm G D{A)} 

and such that; {A*v,u) = {v^Au),'iu G D{A)^'iv G D{A*). Indeed the linear 
form u {v,Au) defined on D{A) for each v G D{A*), has a unique extension 
p G X' = X, and we set: p = A*v. 

Obviously, G(A*) is always closed. Moreover, it is immediate to check that if 
B G L{X), then {A + B)* = A* + B*. 

In the Hilbert space setting , m-dissipative operators can be characterised rather easily 
. First the following proposition follows from elementary duality properties 

Proposition 2.1.6. A linear operator Aon X is dissipative in X if and only if 

Vu G D{A), {Au, u) < 0. 

In addition if A is m-dissipative on X, then D{A) is everywhere dense in X. 

The following result is often useful, especially the two corollaries: 

Proposition 2.1.7. Let Abe a linear dissipative operator on X, with dense domain. 
Then A is m-dissipative if and only if A* is dissipative and G{A) is closed. 

Corollary 2.1.8. If A is self-adjoint in X, in the sense that D[A) = D{A*) and 
A*u = Au, for all u G D{A), and if A < 0 (which means {Au,u) < 0 for all 
u G D{A), Then A is m-dissipative. 

Corollary 2.1.9. If A is skew-adjoint in X, in the sense that D(A) = D(A*) and 
A*u = —Au, for all u G D{A), then A and —A are both m-dissipative. 

2.1.3 Examples in the theory of PDE 

In this paragraph, we recall some basic facts from the linear theory of partial differential 
equations which shall be used throughout the text. The definitions of Sobolev spaces 
and the associated norms are the standard ones as can be found in |[3l. In particular, H 
being an open set in R^, we shall use the spaces 


= {uG L^{D), D,u G L^{D), Mj : |jj < m} 
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endowed with the obvious inner product 

H^(fl) = completion of C°° functions with compact support in fl for the iJ™ 
norm. 

We recall the Poincare inequality in i?o (12) when 12 is bounded : 

Vw G HqID), f \Vw\'^dx > Ai / \w\'^dx, 

Jq Jq 

where Ai = Ai (12) is the first eigenvalue of (—A) in Hq (12) . We are now in a position 
to describe our basic examples. 

Example 2.1.10. : The Laplacian in an open set of theory. 

Let 12 be any open set in R.^, and H = L^(12). We define the linear operator B on H 
by 

D{B) = {uG iLo(12), Au G L^(12)}, 

Bu = Au, Vm G D{B). 

Then B is m-dissipative and densely defined. More precisely B is self-adjoint and 
B < 0. In addition if the boundary of 12 is bounded and C^, then 

D{B) = 2/2(12) niJo^ (12), 

algebraically and topologically. 

Example 2.1.11. : The Laplacian in an open set of : (7° theory. 

Let now 12 be any open set in R^. We consider the Banach space 

X = C'°(12) = {uG C{Ti),u = 0on912} 

endowed with the supremum norm and we define the linear operator A by 

D{A) = {u G XnH^{n),Au G X};Au = Au,yu G D{A). 

Then if the boundary of 12 is Lipschitz continuous, A is m-dissipative and densely 
defined on X. 

Example 2.1.12. : The wave operator on H^{n) X L2(12). 

Let 12 be any open set in R^ and X — //q (12) x L‘^{Tl). The space A is a real Hilbert 
space when equipped with the inner product 

{{u,v),{w,z)) = / {XuS/w + vz) dx, 

Jn 

inducing on A a norm equivalent to the standard product norm on Hq{0) x 
We define the linear operator A on A by 

D{A) = {(u,r;) G A, Am G L^{Ul), v G i2o(12)} 

A{u,v) = (m, Am),V(m,m) G D{A). 

Then A is skew-adjoint in A, and in particular A and —A are both m-dissipative with 
dense domains. 
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2.2 The semi-group generated by m-dissipative opera¬ 
tors. The Hille-Yosida-Phillips theorem 

2.2.1 The general case 

Let X be a real Banach space and let A be a linear, densely defined, m-dissipative 
operator on X. The following fundamental Theorem is proved for instance in ll74l[82l . 

Theorem 2.2.1. There exists a unique one-parameter family T{t) G L{X) defined for 
t > 0 and such that 

(1) Tit) G LiX) and ||r(f)|U(x) < l,Vf > 0. 

(2) r(0) = I, 

(3) Tit + s)= T(f)r(s), Vs, t > 0. 

(4) For each x G DiA), it(f) = Tif)x is the unique solution of the problem 

M G (^([O,+oo); I?(A)) n +oo); X) 

< u'it) = Auit), Vf > 0 

m ( 0 ) = X 

Finally, for each x G DiA) and t>Q, we have: Tit)Ax = ATif)x. 

2.2.2 Two important special cases 

In this paragraph, we assume that X is a (real) Hilbert space. The following two results 
can be considered as refinements of Theorem l2.2.1l 

Theorem 2.2.2. Let A be self-adjoint and < 0. Let x G X, and u(f) = Tit)x. Then u 
is the unique solution of 

uG C'([0,+oo);X)nC'((0,+oo);L>(A))nC'i((0,+oo);X) 

< u'it) = Auit), Vf > 0 
u(0) = x 

Remark 2.2.3. Theorem 12.2.21 means that T(f) has a ’’smoothing effect” on initial 
data. Indeed, even if a: G DiA), we have T(f)x G DiA), for all f > 0. As a 
basic example, let us consider the case X = A defined by D(A) = {m G 

Tfp (H), Alt G L^(H)}, Au = Alt, Vit G DiA) where H is a bounded open set in 
and the boundary of H is smooth. Theorem l2.2.2l here says that for each itg G L^(H)}, 
there exists a unique solution 

It G ^([o, +oo), L^(H)) n ^(0, + 00 , H'^in) n H^in)) n ^^(o, +oo, 
of: 

Ut = Xu ; u(0) = itQ. 

Actually a much stronger smoothing property holds true since by iterating the proce¬ 
dure we prove easily that ii(f) G Z)(A") for all n G N and f > 0. In particular it(f,.) 
is smooth up to the boundary. 
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A somewhat opposite situation is that of isometry groups generated by skew-adjoint 
operators. 

Theorem 2.2.4. Let A be skew-adjoint. Then T (t) extends to one-parameter group of 
operators T[t) : K —>■ L{X) such that 

(1) Vx e X, T{t)x e CiR,X). 

(2) Vx S X,Vt € K, ||T(f)x|| = ||x||. 

(3) Vs G R,\ft G R,T(t + s) = T{t)T{s). 

(4) For each x G D{A), u(t) = T(t)x is a solution ofu'{f) = Au(t), Vt G K. 

Example 2.2.5. Let X = Hq{TI) x and let A be as in Example 12. 1.1 21 We ob¬ 
tain that for any (uq, vq) € X , there is a solution u G C(]R, i/p (R, L^(n)) fl 

C^iR,H-\Q)) of; 


Utt = Au; u(0) = uo,ut(0) = Vo¬ 


lt can be shown that u is unique. 

2.3 Semilinear problems 

Let be a real Banach space, let A be a linear, densely defined, m-dissipative operator 
on X, and let T(t) be given by Theorem l2.2.1l The following Theorem is quite similar 
to the construction of the flow associated to an ordinary differential system and is the 
starting point of the theory of semilinear evolution equations. 

Theorem 2.3.1. Let F : X ^ X be Lipschitz continuous on each bounded subset of 
X. Then for each x € X, There is t{x) G (0, -|-c»] and a unique maximal solution 
u G t(x)), X) of the equation 

u{t) = T[t)xf T{t — s)F{u{s)) ds 

Jo 

The number t(x) is the existence time of the solution , and satisfies the following al¬ 
ternative: either t(x) = oo and the solution u with initial datum x € X is global (in 
X); or t(x) < oo and the solution u with initial datum x G X blows up infinite time 
(in X). In the latter case we have 

||u(t)|| —> -foo as t — > t(x). 

In the theory of semilinear evolution equations, a basic tool to establish global 
existence, uniqueness, boundedness or stability properties of the solution will be the 
Gronwall Lemma (cf. Lemma [1.2.11 i. 

2.4 A semilinear heat equation 

Let Tl be any open set in R^ with Lipschitz continuous boundary dTt , and let us 
consider the equation 

Ut — Xu -\- f{u) = 0 in R'*" x Tt, u = 0 on R"*" x dTt (2.1) 
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where / is a locally Lipschitz continuous function: R —K with /(O) = 0. It is natural 
to set 

X = C'°(f2) = {u e C{n),u = Oon^fl} 

and to introduce the semi-group T{t) on X associated to the homogeneous linear prob¬ 
lem 


ut — Alt = 0 in R'*" X n, it = 0 on R"*" x dU 

In fact here T [t) is the semi-group generated by the operator A of Examr)Ie l2.1.11l Let 
G A: by Theorem 12. 3.1 1 we can define T{tf) < oo and a unique maximal solution 
solution It G <^([0, T{ip)), X) of the equation 

u{t) = T{t)x-\- f T(t — s)F{u{s))ds 
Jo 

with F : X ^ X given by {F{u)){x) := —f(u(x)) for all x in the closure of W. Then 
u can be considered as the local solution of (12. Il l with initial condition it(0) = ip'm X. 
The following simple result will be useful later on. 

Proposition 2.4.1. Let/satisfy the condition 

Vs G R with |s| > C, f{s)s > 0 (2.2) 

Then we have for any ip G X 

r((/?) = 00 and sup ||it(f)||Loo < Maa;{(7, ||(/3||lcxi} < oo (2.3) 

t>o 

where u is the solution of (EB with initial condition it(0) = (p . 

Proof. Let M = Max{C, Hv’IIloo} and let us show for instance that u(t, x) < M on 
(0, X 11. Introducing z = u — M, we have 

- Az = /(M) - f{u) - f{M) < f{M) - f{u) 

since f{M) > 0. In addition it can be shown that 

u G 0(0,T{ip); H^{n) n H^{n)) n (7^(0, T{:py,L^u)) 

and then 

{d/dt) f \z'^\'^dx = 2 ( z'^ztdx = 2 ( z^{Az + f{M) — f{u) — f{M))dx 

J Q, J Q. J Q, 

<—2 j Vz“''.Vz dx + 2 f z^\f{M) — f{u)\dx 

Jq, Jn 

Because / is locally Lipschitz and it is bounded on (0, t) x ft for each t < t{‘p), we 
have 

\f{M) - < K{t)\z{t,x)\ on (0,1) xfl 

Then by using the identities z = z+ — z“ and z+.z“ = 0, Vz+.Vz“ = 0 almost 
everywhere, we obtain: 

{d/dt) f \z~^ydx <—2 [ \\Xz'^\y‘ dx F2K{t) f \z^y dx 
J VL h Cl J Cl 

The inequality u(l, x) < M on (0, t{‘p)) x H now follows easily by an application of 
Lemma [1.2. II since z“''(0,a;) = 0. Similarly we show u{t,x) > —M on (0 ,t((^)) x 

n. □ 
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2.5 A semilinear wave equation with a linear dissipa¬ 
tive term 

Let SI be any open set in with Lipschitz continuous boundary dfi , and let us 
consider the equation 

Utt — Art + 'yut + f(u) =0 in K"*" x fl, u = 0 on R.’*" x dUl (2.4) 

where / is a locally Lipschitz continuous function; R. —R. with /(O) =0 satisfying 
the growth condition 


|/'(u)| ^^(l + litr), a.e. on R (2.5) 

2 

with r > 0 arbitrary if = 1 or 2 and 0 < r < ——^ if iV > 3. It is natural to set 

a: = hI{vl) X l^{vl) 

Let us denote by /* the mapping dehned by 

f*{{u,v)) = {0,-f{u)),V{u,v) G a:. 

The growth condition (12.51) together with Sobolev embedding theorems imply that 
f*{X) C X-, f* '■ X —X is Lipschitz continuous on bounded subsets. 

We also dehne the operator L G L{X) given by 

r((u,r;)) = ( 0 , 71 ;), V(m,u)gA:. 

Finally let T{t) (cf. Theorem 12. 2 .41 with A as in example 12.1.121 in X = x 

L^(n)) be the isometry group on X generated by the linear wave equation 

Utt — Art = 0 in R"^ x 12, u = 0 on R"*" x dTt 

For each {(fij'ip) £ by Theorem 12. 3. II we can dehne a unique maximal solution 
solution?/ = {u,ut) G C([0, t( 93,'0)); X) of the equation 

uit) = i,) + f Tit - snrauis) - viuismds 

Jo 

The following simple result will be useful later on. 

Proposition 2.5.1. Assume 7 > 0, and let f satisfy the condition 

Vs G R, F(s) > (—^ + e)s^ — C with e > 0, C > 0 (2.6) 

where F is the primitive of f such that X(0) = 0 and Ai is the first eigenvalue of 
—A in H^iU). Then we have for any ip,ip) G X : rip,ip) = 00 and the solution 
U = (u, Ut) of (12.41 ) such that ?7(0) = [p, ip) satisfies : 

sup \\iuit),utit))\\x < 00. 
t>o 
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Proof. The solutions of (I2.41 i satisfy the energy equality 


7 / f u^{t,x)dxdt + E{u{t),ut{t)) = E(ip,'ip) 

Jo Jn 


with 


•= T / \ [ \Jjix)\^dx+ [ F{ip{x))dx 

In particular since 7 > 0, we find E{u{t),ut{t)) < E{ip,ip) and the result follows 
quite easily from (12.61 1. Indeed, from Poincare inequality we deduce 

Vtu £ Tfg (fl), (1 — 7 ) f iVrupda: > (Ai — 2£) f lu'^dx, 

Jq Jq 

whenever p < 2e/Ai. Then 

E{<f,i:)>{n/2) [ \V^fdx + l [ \i;{x)fdx-C\ni V((p,V’)G 2 f, 

Jn ^ Jn 

and a bound on i? implies a bound in AT. □ 


Chapter 3 


Uniformly damped linear 
semi-groups 


3.1 A general property of linear contraction semi-groups 

Let X be a real Banach space and L any m-dissipative operator on X with dense 
domain. We consider the evolution equation 

u' = Lu{t), t>0 (3.1) 

For any uq G X, the formula u{t) = S(t)uo where S{t) is the contraction semi-group 
generated by L defines the unique generalized solution of (13.11 1 such that u(0) = uq. 
We recall the following simple property ; 

Proposition 3.1.1. For all f > 0, let us denote by ||S'(f)|| the norm of the contractive 
operator S(t) in L{X). Then ||5'(f)|| satisfies either of the two following properties 

(1) For all t> 0, \\S{t)\\ = 1. 

(2) 3e > 0,3M > 0,forallt> 0, ||S'(f)|| < 

Proof. The function ||S'(f)|| is nonincreasing. If for some T > 0 we have |lS'(t)|| = 1 
for f G [OjT) and ||S'(T)|| = 0, then Ve > 0,Vf > 0,||S'(f)|| < M{e)e~^* with 
M(e) = Assuming, on the contrary, that for some r > 0 we have 0 < ||S'(r)|| < 

1, for each f > 0 we can write t = nr + s, with n G N, 0 < s < r. Then ||S'(f)|| < 

||S'(t)||” and we obtain (2) with e = — and M = = 1/||5'(t)||. □ 

T 

3.2 The case of the heat equation 

The linear heat equation can be studied in many interesting spaces. Its treatment is 
especially simple in the Hilbert space setting of examr)le f2. 1.101 However, in view of 
the applications to semilinear perturbations the Cq- theory is more flexible. Let us start 
with the Hilbert space setting : following the notation of example 12.1.101 we denote 
by S{t) the semi-group generated by i? in iT = We have the following simple 

result. 
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Proposition 3.2.1. Let Ai = Ai(ri) be the first eigenvalue of {—A) in iJg (12). Then 

\\Smc{H)<e-^^\ Vt>0. (3.2) 

Proof. Let ip € D{B), and consider 

fit) = \\Sit)p>\\H)2, Vt > 0. 


We have 


e ^^^*f'{t) = 2Ai / u{t,x)^dx + 2 / u{t,x)u'{t,x)dx 

JQ JQ 

= 2Ai / u{t,x)‘^dx + 2 / u(t,x)Au{t,x)dx 
Jn Jn 

= 2^Ai J u{t,x)^dx — J \Vu{t,x)\^dx'j <0. 

Hence 

Vt > 0,V(^ e DiB). 

The result follows by density. □ 

We now assume that H is bounded with a Lipschitz continous boundary and we use 
the notation of Example 12. Till Let T{t) denote the semi-group generated by A in X. 
Since X C H with continous imbedding and G{A) C G{B), it is classical, using the 
Hille-Yosida theory, to prove 

Vv? GX,yt> 0, T{t)ip = S(t)ip (3.3) 

In particular we have; ||S'(f)i^||// < for each f > 0 and Lp £ X. The 

following property of uniform damping in X will be more interesting for semilinear 
perturbations 


Theorem 3.2.2. Let Ai = Ai(r2) be the first eigenvalue of {—A) in iJg (12). Then 


wsmax) < Vf>0 


(3.4) 


with 


,Ai|12|2/w, 


(3.5) 


In the proof of Theorem l3.2.2l we shall use a rather well-known smoothing property 
of S{f) in LP spaces. Denoting by || • ||p the norm in Lp{Lt) , we recall 

Proposition 3.2.3. Let 1 < p < q < oo. Then 

||-5'(f)¥3||5 < (^)^^p"’^||‘/3||p,Vf > 0,V(/3 G X. 


A possible proof, omitted here, relies on the explicit form of the heat kernel in 
together with a comparison principle. 
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Proof of Theorem l3.2.2l Let ip € X and T > 0. First for 0 < f < T, we have trivially 

llS'(f)(/3||oo < llv^lloo < 

Then if t > T, we find successively, applying first Proposition 13.2.31 with p = 2 and 
q = oo 


< (by Proposition 13.1.il l 


Then the estimate follows by letting T 


47r 


□ 


Remark3.2.4. Actually (I3.41 i is not valid with M = 1. More precisely, if ||S'(f)||£(x) < 
with m > 0, we must have M' > 1. Indeed, let ip £ 'D{D) be such that 
ip = \ near Xq G Tt and \\p\\x = 1, and let u{t) = S{t)p. It is then easily verified that 
u £ C°“([0, oo) X n). Consequently Mt(0, x) = 0 near xq. Hence, for any e > 0 and 
any x close enough to xq, we find 


u(t, x) > 1 — et, 

for all t sufficiently small; in particular 

||■w(^)||x > 1 - £f 

for t small. This estimate with e > 0 arbitrary small is not compatible with 
||5'(f)||£(x) < e~^^, for whatever value p > 0. 


3.3 The case of linearly damped wave equations 

We have the following result 

Proposition 3.3.1. Let D be a bounded domain in Consider the equation 

Utt — Am + Xut = 0 in x H, u = 0 on R^ x dO (3.6) 

Then, denoting by || • || the norm in iFp (H) and | • | the norm in L^(Ll), for any 
solution u of (O we have 

IKi)ll + kWI<C^(lk(0)|| + k(0)|)e-^‘ (3.7) 


for some C,5 > 0 . 

This result is a special case of the following more general statement. Let A be a 
positive self-adjoint operator with dense domain on a real Hilbert space H with norm 
denoted by |.| and inner product denoted by (.,.). A is assumed coercive on H in the 
sense that 


3a > 0, Vm £ D{A), (Au, u) > a\u\‘^. 












3.3. THE CASE OE LINEARLY DAMPED WAVE EQUATIONS 


29 


We introduce V := the closure in H of D{A) under the norm 

p{u) := {Au,u)^. 

The norm p extends on V and we equip V with the extension of p, denoted by || || so 
that 

Vu G V, Hull = |A^/^u| 

where G L(l/; iJ)nL(I?(A); TG) is the unique nonnegative square root of Gl. The 
duality product between V and its topological dual V extends the inner product on H 
in the following way: 

V(/, v)€ H xV, (/, v)v',v = if, v). 

In particular we have 

V{u,v) G D{A) X V, {Au,v)v',v = {Au,v) = (A^^'^u, A^^'^v). 

In particular by the definition of the standard norm on V we have 

\/uGD{A), \\Au\\v'<\A^/^u\ = \\u\\. 

By selecting v = uwe even obtain 

VugD{A), ||xlu||y'= ||u||. 

By Lax-Milgram’s theorem the extension A of ^ by continuity on V is bijective from 
V to V' and in addition, A satisfies 

V(u,v) G V X V, {Au,v)v',v = [A^/'^u, A^/'^v) 

so that A becomes by definition the duality map from V to 1 ^'. Finally, denoting by 
II • II* the standard norm on of V' we remark that 

v/gf', ll/ll* = ||a-VI|. 

Let now B G L{V-,V')h& such that 

Vu G V, (Bv, v) > 0. 

We consider the second order equation 

u" + Au + Bu = 0. 

and the energy space E = V x H is equipped with the Hilbert product space norm. 
Proposition 3.3.2. The unbounded operator on E defined by 

D{L) = {iu,v)GV xV] Au + BvGH) (3.8) 

L{u,v) = {v,—Au — Bv) V{u,v) G D{L) (3.9) 


is m - dissipative on E. 
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Proof. We denote by (,) the inner product in E. First L is dissipative on E . Indeed 
for any U = {u, v) G D{L) we have 

{LU,U) = {v,u)v + {—-^u — Bv,v)h 

= {—Au — Bv,v)v',v 

= {—Bv, v)v',v < 0- 

In order to prove that L is m-dissipative on E we consider, for any (/, g) G E the 
equation 

{u,v) G D{L); -L{u,v)-\-iu,v) = {f,g) 

which is equivalent to 

{u,v)gVxV; —v-\-u = f; AuBvv = g 

or in other terms 

{u,v) G V xV; u = f -\-v; Aw + Bv v = g — Af 

Assuming we know that the operator C = A + B + Jis such that C{y) = V we 
conclude immediately that 

[I -L)D{L) = E 

and therefore L is m-dissipative as claimed. The property C'(V) = V' is an immediate 
consequence of the following elementary lemma □ 

Lemma 3.3.3. Let V be a real Hilbert space and C G L{V, V'). Assume that for some 
r] > 0 we have 

Mv gV, {Cv,v)v',v > gWvW'^. 

Then C(V) = V! 

Proof. First C(V) is a closed linear subspace of V . Indeed if /„ = Cw„ G C(V) and 
/„ converges to f G V' we have for each (m, n) the inequality 

llt^n 't’m II ^ ~{fn /mi V\V r ||Wyj || ^ Wfn /m |U 

Hence w„ is a Cauchy sequence in V and its limit v satisfies Cv = f. Now if C(V) 

V' there exists a non-zero vector w G V such that 

\/vGV, {Cv,w)v>,v = Q 

By letting v = w we conclude that w = 0, a contradiction. □ 

Proposition 3.3.4. Let A, V and H be as above. Let B G £(V, V') satisfy the follow¬ 
ing conditions 

3 a > 0 , VvGV, {Bv,v)v',v > a|wp 

30 0, \/vGV, \\B{v)\\^y, <C{{Bv,v)v',v + \v\^)■ 

Let u G C^(0, -foo, V) n <3^(0, -boo, V') be a solution of 

u" -b Au -b Bu' = 0. 


There exists some constants C > 1 and 7 > 0 independent of u such that 


V>0, ||M(f),'u'(f)|l£; < Ce '>'*||m(0),u'(0)||£;- 
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Proof. We consider for alH > 0 and e > 0 small enough 

HS) = \\u\t)\\^ +\\Aiu{t)f+e{uit),u\t)) 

and we compute 

HUt) = +s\\u\t)f+e{u"it),u{t)) 

= -{B{u'{t)),u'{t)) +e\\u'{t)r-e\\A^u{t)r-e{Bu'{t),u{t)) 

< -{B{u\t)),u'(t)) + e\\u\t)r - e\\A^u{t)f + r^eWuml + 

1 ] 

Cr C 1 

< [-l + —){B{u'{t)),u'{t)) +e{l + -)\\u'{t)r-s{l-rj)\\A^u{t)\\^. 

rj r] 

Choosing for instance rj = i/e and letting e small enough we obtain first 

On the other hand it is not difficult to check for e small enough the inequality: 

(1 - Me)\\u\t)r + \\Aiu{t)r < HS) < (1 + Me)\\u'{t)r + . 

where M is independent of the solution u as well as t and e. This concludes the proof. 

□ 


Remark3.3.5. If (u(0), w'(0)) G £>(L), then clearly u G C^(0,+oo, C)nC'^(0,+oo, C'). 
By density. Proposition 13.3.41 means that the semi-group generated by L is exponen¬ 
tially damped in E. In particular Proposition l3.3.1l folIows as a special case. 






Chapter 4 

Generalities on dynamical 
systems 

4.1 General framework 

Throughout this paragraph, (Z, d) denotes a complete metric space. 

Definition 4.1.1. A dynamical system on (Z, d) is a one parameter family {5'(t)}t>o 
of maps Z ^ Z such that 

(i) yt > 0, S(t) € C'(Z, Z); 

(ii) S{0) = Identity; 

(Hi) Vs, t > 0, S{t + s) = S{t) o 5'(s); 

(iv) Vz G Z,S{t)z G (^([0, +oo), Z). 

Remark 4.1.2. In the sequel we shall often denote S(t)S{s) instead of S(t) o S{s). 

Remark 4.1.3. If F is a closed subset of Z such that S(t)F C F for all t > 0, then 
{5(f)/F }t>o is a dynamical system on {F, d). 

Definition 4.1.4. For each z G Z, the continuous curve t —>■ S(t)z is called the 
trajectory of z (under S{t)). 

Definition 4.1.5. Let z G Z. The set 

uj{z) = {y G Z, 3tn —>■ + 00 , S(tn)z y as n ^ +oo} 

is called the uj-limit set of z (under S(t)). 

Proposition 4.1.6. 'We also have 

w(2)=n ij{5'(f)0}. 

s>0 t>s 

Proof. Immediate according to Definition |4T3] □ 
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Proposition 4.1.7. For each z G Z and any t > 0, we have 


ui{S{t)z) = uj{z)] 

(4.1) 

S{t){uj{z)) C oj{z). 

(4.2) 

In addition, if is relatively compact in Z, then 

t>0 


S{t){uj{z)) = uj{z) 0. 

(4.3) 


Proof, a) (14. Il l is an immediate consequence of Proposition |4]T]6] 

b) Let y G u!{z). There is an infinite sequence —>■ +oo such that as n —>■ +oo, 

S{tn)z —>■ y. For each f > 0, setting t„ = + f, we find S{Tn)z —>■ S{t)y, therefore 

S(t)y G w(^); hence (I4.21 i. 

c) Finally, assume to be precompact in Z. There is an infinite sequence 

t>o 

tn —>■ +00 and y G Z such that as n —>■ +oo, S{tn)z —>■ y. Thus y G oj{z) and 
oj{z) ^ 0. To establish the inclusion uj{z) C S{t){uj{z)), let us consider y G uj(z) 
and tn —?■ +00 such that S{tn)z —?> y. let Tn = tn — t. By possibly replacing by a 
subsequence, we may assume S{Tn)z w G oj{z). Hence by continuity of S{t) 

S(t)w = S{t) lim S{t„)z = lim S{tn)z = y, 

n—>-+oo n^+oo 

and (14.31 1 is completely proved. □ 

In the sequel, a subset B of Z being given , we shall denote by 

d{z,B) := inf d{z,y) 
yes 

the usual distance in the sense of {Z, d) from a point z G Z to the set B. Using this 
notation we can state 

Theorem 4.1.8. Assume that is relatively compact in Z. Then 

t>o 

(i) S{t){uj{z)) = oj{z) ^ Ih, for each t > 0; 

(ii) uj(z) is a compact connected subset of Z; 

(Hi) d{S{t)z,uj{z)) ^ 0 as t ^+ 0 O. 


Proof, (i) is just (14.31) . Moreover, for all s > 0, |J{S'(f)z} is a nonempty com¬ 
ics 

pact connected subset of Z . Proposition 14.1.61 therefore implies that a;(z) is a com¬ 
pact connected subset of Z as a nonincreasing intersection of such sets: this is (ii). 
To check (iii), let us asssume that there exist -(-oo and e > 0 such that for 

all n, d{S{tn)z,uj{z)) > e. By compactness and by the definition of oj{z), there is 
a point y G uj{z) and a subsequence —>■ +c» for which S{tn')z —>■ y. Hence 

d{S{tn')z, uj{z)) —>■ 0, a contradiction which proves the claim. □ 
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We now introduce the basic example of dynamical systems to be studied in this 
book. Let X he a real Banach space, let A be a linear, densely defined, m-dissipative 
operator on X, and let F : X —> X be Lipschitz continuous on each bounded subset 
of X. As recalled in Theorem 12. 3. II for each x G X, there is t{x) G (0, +c»] and a 
unique maximal solution u G (^([0, t{x)), X) of the equation 

u{t) = T{t)x+ f T{t — s)F{u{s)) ds VfG[0,r(a:)) (4.4) 

Jo 

where T{t) is the semigroup generated by A (cf. Theorem l2.2.1l) and the number r(a;) 
is the existence time of the solution. For x G X and t G [0, T(a;)),we set 

S{t)x = u{t). 

Let Y C X he such that for some M < +oo we have 

r(y) =+oo,Vj/e F; (4.5) 

||F(f)y|| <M,Vj/GF,Vf>0. (4.6) 

We set Z = (^{ 5 (^) 2 ;} and we denote by d the distance induced on Z by the 

y£Y t>s 

norm of X. 

Lemma 4.1.9. We have the following properties 

(i) t{z) = + 00 , Vz G Z; 

(ii) ||S'(f)z|| < M, Wz G Z, Vf > 0; 

(Hi) S{t)z G Z.'iz G Z, 'it > 0. 

Proof. Let y G Y. Then if u{t) = S{t)y is the solution of (14.41 1 with x = y a straight¬ 
forward calculation shows that for any s > 0, v{t) = u{t -f s) is the solution of (14.4b 
with X = u{s). Therefore, 

S{t)S{s)y = S{t){u{s)) = u{t + s), is, t>0. 

Consequently T{S{s)y) = -foo for aWy gY and each s,t > 0 and ||S'(f)5'(s)y|| < M 
for all y G F and each s,t > 0. Now let z G Z. There exists a sequence (f„) in [0, -(- 00 ) 
and a sequence (y„) in F such that S'(t„)y„ —>■ z as n —>■ -foo. Pick T < r(z). Of 
course we have by GronwalTs Lemma nemma fl.2.1b : 

S{f)S{tn)yn —>■ S{f)z as n —>• -f 00 , uniformly on [0,T]. (4.7) 

In particular ||5'(f)z|| < M,it G [0,r]. Since T < r(z) is arbitrary, we deduce 
first (i), then (ii). Finally (iii) follows as a consequence of (14.7b . □ 

Theorem 4.1.10. {S'(f)}t>o is a dynamical system on {Z, d). 

Proof. First S'(O) = Identity. Moreover for each z G F, if z^ G F and Zn ^ z 
as n ^ -foo, as a consequence of the Gronwall Lemma (lemma rL2. lb we obtain 
classically : 

S(f)zn —>■ S{t)z as n —>■ -foo, uniformly on [0,r] 
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for each finite T . In particular S{t) £ C{Z, Z) for all t > 0. Moreover for each y G Z, 
the calculation performed in the proof of Lemma l4.1.9l shows that 

S{t)S{s)y = S{t + s)y 

for all s,t > 0. Finally by construction we have S{t)z £ C{\0, +oo), Z) for each 
z G Z. Hence the result. □ 

As a particular case of Theorem l4.1.101 we can choose X = > 1. For each 

vector field F G we consider the (autonomous) differential system 

u'it) = F{u{t)) (4.8) 

and its integral curves u(t) =: S{t)x defined for t G [0,T(a;)). Theorem 14. 1.1 01 savs 
that if T{y) = +(X) and the corresponding local solution u{t) remains bounded for 
t > 0, then t{z) = +c» for each z G Z := u.(]R+) and the restriction of S{t) to Z 
(endowed with the distance associated to the norm) is a dynamical system. To see this 
we apply Theorem l4. 1. lOl with A = 0 and Y = {y}. 

Other important examples of dynamical systems will be associated to the partial 
differential equations studied in Chapter |2l Their properties will be studied precisely 
in the next chapter. 


4.2 Some easy examples 

In the first section (Theorem l4.1.81 l. we showed that the w-limit set of a precompact tra¬ 
jectory u{t) = S{t)z is a continuum invariant under S{t) and which (by construction!) 
attracts the trajectory as t ^ +c». In some cases this gives directly a convergence 
result. As a first easy case we have 

Proposition 4.2.1. Ifui{z) is discrete, there exists a G Z such that d{S{t)z, a) ^ 0 as 
t —>■ +c» 

Proof. This is an immediate consequence of Theorem 14.1.81 Indeed, uj{z) , being 
compact and discrete is finite. But a connected finite set is reduced to a point. □ 

As an example let us consider the second order ODE 

'' I ' I 3 n 

U +U+U — u = u. 

All solutions are global and an immediate calculation gives: 

(d/dt)l(l/2)u'^ + (l/4)u^ - (l/2)u2] = _u '2 < 0. 

Hence we can define the dynamical system generated on the whole of by setting 
[/(t) = (u(t), u'{t)) and writing the equation as a first order system. The function 1 1 —>■ 
+ (l/4)tt^— (l/2)it^](f) is nonincreasing along trajectories. Consequently it 
has a limit as t tends to infinity and, as a consequence, each trajectory (v, v’) contained 
in the w-limit set of a given trajectory satisfies automatically 


0 = {d/dt)[{l/2)v'^ + {l/A)v^ - (1/2)u 2] = -u'2. 
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It follows, since this implies v' = 0, that the w-limit set of any trajectory consists of 
stationary points and is therefore contained in {0,1, —1} x {0}. By connectedness, 
the w-limit set reduces to a singleton {{z, 0)} with z = 0,1 or (-1). Therefore every 
solution has a limit at +c». 

Actually the argument which we gave above in this special case is general for sys¬ 
tems having what will be called a ’’strict Liapunov function”. On the other hand already 
in there are many examples of systems with non-convergent bounded trajectories. 
For instance the basic second order equation 

u" + uj'^u = 0 

has no convergent trajectory except m = 0. Here instead of a Liapunov function we 
have an invariant energy, and the w-limit set of any solution other than the single equi¬ 
librium point (0,0) does not intersect the set of equilibria. 


4.3 Convergence and equilibrium points 

In this section we introduce some general concepts which will be used throughout the 
text. 

Definition 4.3.1. Let z £ Z. The trajectory t —>■ S(t)z is called convergent if there is 
a G Z such that 

lim d{S{t)z,a) = 0. 

t—¥-\-OQ 

Definition 4.3.2. A point z G Z is called an equilibrium point (or equivalently a 
stationary point) of the dynamical system S(t) if {z} is invariant under S{f), i.e. 

'it > 0, S{f)z = z. 

The following property is now obvious 

Proposition 4.3.3. If a trajectory of the dynamical system S{t) is convergent, the limit 
is always a stationary point. 

Proof. This is an immediate consequence of Proposition 14.1.71 Indeed if a trajectory 
converges, it is precompact and the omega-limit set is an invariant singleton. □ 

Remark 4.3.4. As a trivial consequence of Proposition l4.3.3l a necessary condition for 
a precompact trajectory to be convergent is that its w-limit set be made of equilibria. In 
chapter 6 we shall study an important class of systems for which the w-limit set of all 
precompact trajectories is reduced to equilibria. Then if the set of equilibria is finite, 
convergence follows from Proposition 14.2.11 On the other hand an important part of 
the book will be devoted to the harder case of a continuously infinite set of equilibria. 


4.4 Stability of equilibrium points 

Another important concept concerning equilibria (and more generally trajectories) of a 
dynamical system is the concept of stability as defined by Liapunov. 








4.4. STABILITY OF EQUILIBRIUM POINTS 


37 


Definition 4.4.1. An equilibrium point a of the dynamical system S{t) is called stable 
(under S(t)) if 

Ve >0, 3(5 > 0, Vz G Z, d{z, a) < 6 Vt > 0, d{S{t)z, a) < e. 

Otherwise we say that a is unstable. 

The following result, relted to the concept of Liapunov function, provides a general 
stability criterion applicable even to infinite dimensional systems. 

Theorem 4.4.2. Let a G Z be an equilibrium point of the dynamical system S{t) and 
U be an open subset of Z with a G U such that for some V G C{Z) we have 

VrG(0,ro), min V{u)>V{a) (4.9) 

d{u,a)—r 

Vu gU, V< > 0, V{S{t)u < V{u) 

Then a is a stable equilibrium point of the dynamical system S(t). 

Proof. Let r > 0 be such that B{a, r) C U and let 

c := min V{u)>V(^a) 

d{u,a)—r 


Let 

W = {u G B{a,r),V{u) < c} 

It is clear that W is open with a G W. In addition if uq G W, u{f) = S{t)uo satisfies 

Vf > 0, u{t) G W 

Indeed if this property fails for some uq G W, we can consider 

to = inf{f > 0, u(t) ^ W}. 

We have V (u(to)) < V (uq) < c and since W is open the only possibility is d{u{to), a) = 
r, a contradiction with the definition of c. The result is now immediate since r > 0 can 
be chosen arbitrarily small. □ 

Under the hypothesis that balls with finite radius are compact subsets, we obtain 
the following result applicable in finite dimensions. 

Corollary 4.4.3. Assuming that closed balls with finite radius are compact subsets of 
Z, let a G Z be an equilibrium point of the dynamical system S{t)and U be an open 
subset of Z with a GU such that for some V G C{Z) we have 

'iu G U,u f a ^ V{u) > V{a) 

yuGU, Vf > 0, V{S{t)u < V(^u) 

Then a is a stable equilibrium point of the dynamical system S(t) . 

Proof. Let r > 0 be such that B{a,r) C U : as a consequence of the compactness of 
closed balls we have (14.91 ). The result is now an immediate consequence of Theorem 
l4A2l □ 
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Definition 4.4.4. An equilibrium point a of the dynamical system S(t) is called asymp¬ 
totically stable (under S(t)) if it is stable and in addition 

3^0 > 0, yz £ Z,d{z,a) < So lim d{S{t)z,a)=0 

i^+oo 

Remark 4.4.5. The first order ODE 

u' -lu^ — u = 0 

generates a dynamical system on Z = R which has a set of 3 equilibria {—1,0,+!}. 
It is easy to verify that all trajectories of this system are convergent, positive initial 
data lead to a trajectory converging exponentially fast to +1, negative initial data to a 
trajectory converging exponentially fast to -1. Therefore +1 and -1 are asymptotically 
stable, whereas 0 is unstable. It is not too difficult to check that the equilibria (1,0) 
and (—1,0) are also asymptotically stable for the system generated in Z = by the 
second order ODE 

u" 4- u' 4- — u = 0 

considered in the previous section, whereas in this case the set of initial data leading 
to a trajectory tending to (0,0) is a ID curve separating the attraction basins of the 2 
stable equilibria. Hence (0,0) is also unstable in this case. 

In the case of the basic oscillator governed by 

u” + = 0 

the only equilibrium is 0 which is stable (with 5 = e since we have an isometry group 
on Z = R^) but not asymptotically stable. This result can also be viewed as a special 
case of theorem l4.4.2l with V{u, u') = l(u'^+a;tt^). The same argument holds true for 
the wave equation with V the usual energy functional. We remark that except for the 
initial data (0, 0), the omega-limit set does not cross the set of equilibria. In fact if the 
omega-limit set of a trajectory contains a stable equilibrium point, the trajectory must 
converge to this point. This makes the study of convergence somewhat easier when the 
dynamics is unconditionally stable, a typical case being contraction (or more generally 
uniformly equicontinous) semi-groups which will be studied in Chapter^ 




Chapter 5 


The linearization method in 
stability analysis 


When looking for stability of an equilibrium point a for an evolution equation U' + 
AU = 0, a natural idea is to examine the nature (convergent or divergent) of the linear 
semi-group generated by the linearized operator DA{a). It is intuitively clear that this 
will work only when the spectrum of DA{a) does not intersction the imaginary axis. 
In this chapter, we hrst describe an extension of the Liapunov linearization method to 
establish the asymptotic stability of equilibria. The perturbation argument developed 
here is applicable, in conjonction with the linear results of Chapter 2, to various semi 
- linear evolution problems on inhnite dimensional Banach spaces. At the opposite, an 
argument essentially coming back to R. Bellman ifTSll allows to deduce instability from 
the existence of an eigenvalue with the ’’wrong” sign. We shall also provide an inhnite 
dimensional version of the linearized instability principle. 

5.1 A simple general result 

Let X be a real Banach space, T{t) a strongly continuous linear semi-group on X, and 
F : X —> X locally Lipschitz continuous on bounded subsets. For any x G X, we 
consider the unique maximal solution u G (^([0, t{x)),X) of the equation 



By a stationary solution of (15. Il l we mean a constant vector a G X such that 



The following result is an easy consequence of the general theory of strongly continu¬ 
ous linear semi-groups. Let L denote the generator of T{t). Then we have 

Lemma 5.1.1. A vector a G X is a stationary solution of (15.11 1 if and only if we have 

a G D{L) and La + F{a) = 0. 

We are now in a position to state the main result of this section 
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Theorem 5.1.2. Assume that for some constants i5 > 0, M > 1 we have 

Vi > 0, ||r(t)|| < (5.2) 

Let a € X be a stationary solution of (iriT) such that 

3Ro >0,3r]> 0 : ||F(m) - i^(a)|| < rj\\u - a|| for ||u - a|| < Rq (5.3) 


with 


Then for all x G X such that 


V < 


M' 


\\x 


a|| < i?i = 


i?o 

M 


the solution u of dSB is global and satisfies 


Vi > 0, ||M(i) — a|| < M\\x — a\\e 


(5.4) 


with : j = S — tjM > 0. 

Proof. On replacing u by it — a and F by F — F{a), we may assume a — 0 and 
F(a) — 0 with ||F(it)|| < ? 7 ||it|| whenever ||it|| < Rg. In particular, setting 

T = Sup{i > 0, ||it(i)|| < i?o} < + 00 , 


we find 


Vi G [0,T), ||M(i)|| < M||x||e-^* + t]M [ 

Jo 


||m(s)|| ds. 


Letting (p(t) = e^‘||it(i)||, we obtain 


(p{t)<Ci + C2 [ ^(s)ds ViG[0,T) 

Jo 

with: Cl = M||x||, C 2 = r]M. By applying Lemma fl. 2. II with A(i) = C 2 we deduce 

Vi G [0,r), e^*|lit(i)|| < M\\x\\e^^K (5.5) 

Since 6 > rjM , we conclude that if M||x|| < Rg, then T = +00 and (I5.51 l holds true 
on [0, + 00 ). This completes the proof of (I5.41 i. □ 

5.2 The classical Liapunov stability theorem 

5.2.1 A simple proof of the classical Liapunov stability theorem 

The object of this paragraph is to give a simple proof of the following well known 
result: 
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Theorem 5.2.1. (Liapunov) Let X be a finite dimensional normed space, and f £ 
{X, X) a vector field on X. Let a € X be such thatfia) = 0 and assume 

All eigenvalues {sj, 1 < J < fc} of Df(a) have negative real parts. 

Then a is an asymptotically Liapunov stable equilibrium solution of the equation 

u' = t > 0. (5.6) 

More precisely : for each S < rj = min {—i?e(sj)}, there exists p = p(5) > 0 and 

l<j<k 

M(d) > 1 such that if\\x — o|| < p{S)), the solution u of (I5.61 l such that u{0) = x is 
global with 

Vt > 0, ||u(i) — a|| < M{5)\\x — a\\e~^^. 

Proof. We consider first the case where a = 0 and / coincides with a linear operator 
A. In this case, the question reduces to the following; 

Lemma 5.2.2. Let X be a finite dimensional complex vector space, A £ L{X) and 
u £ a solution ofu'{t) = Au(t). Then we have 

k 

u{t)=Y,PjitV^* (5-7) 

f=i 

where {sj}i<j<k is the sequence of eigenvalues of A and Pj a polynomial with coef¬ 
ficients in X for all j. 

Proof. By induction on dime (26) = p. 

- If dimc(26) = 1, then j = 1 and A = sil, hence u{t) = 

- If dime (26) = p > 1, assuming that the result is true for all complex vector spaces 
with complex dimensions < p — 1, we set 

v(t) = M(f)e“®^*, 

therefore u is a solution of 

v' = (A — sil)v. 

Then setting Y = R{A — Sil), B = (A — si/)|y and w = v', it is clear that rc is a 
solution of 

r(;£C^(]R,F); w'(t) = Bw{t). 

Since by construction ker(A— si/) 7 ^ {0}, we have R{A—siI) 7 ^ 26 and in particular 
dimc(l^) < dimc(26) — 1 = p — 1 . 

By the induction hypothesis we have 

k 

i=i 

because the eigenvalues of B are of the form sj — si. By integrating we obtain 

k 

w{t) = oi + 

f=i 

then on multiplying by we obtain (15.71) . completing the proof by induction. □ 
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Completion of the proof of Theorem l5.2.1[ Since all eigenvalues of Df{a) =: A 
have negative real parts, it follows obviously from (15.7b that ||e‘"^|| < C{S)e for 
a\l S < r] = mini<j<fc{—i?e(sj)}. Then we apply Theorem 15.1.21 with T{t) = 
and F defined by the formula 

F{u) = f{u) - Df{a){u - a). 


The result follows at once. □ 

5.2.2 Implementing Liapunov’s first method 

Theorem 15.2.11 gives an apparently simple and almost optimal way of checking the 
asymptotic stability of a given equilibrium point of a differential system; check whether 
all (complex) eigenvalues of the linearization at this point have negative real parts. 
However in practice we have to check this property on the characteristic polynomial, 
but as soon as TV > 3 in general the roots cannot be computed . 

Definition 5.2.3. We say that a polynomial P with real coefficients 

N 

P{X) = 

3=0 

is a Hurwitz polynomial if all its zeroes have negative real parts. 

Proposition 5.2.4. If P is a Hurwitz polynomial, then po 7 ^ 0 and for each j G 
{0, ...TV}, we have pjpo > 0. 

Proof. We have 

P{X) = PN IK-^ + Afe) IK^ + pj + ii'j){X + pj — ivj) 

k j 

where all numbers Xk , Pj are positive . But 

(2f + Pj + ivj^i^X + Pj — ik'j') = X^ + 2pjX + p^ + Uj 

The result follows immediately by expanding P. □ 

Remark 5.2.5. The converse of Pror)osition l5.2.4l is false if TV > 2 . If all coefficients 
of P have the same sign, of course P cannot have a positive real root but on the other 
hand the polynomial 

Pe{X) = {X + l)iX^ - sX + 1) = X^ + {1 - e)X^ + (1 - e)X + 1 

has all its coefficients positive for 0 < e < 1 , although the two conjugate imaginary 
roots have imaginary parts equal to |. 

It is sometimes useful to remember the following criterion which we give without 
proof: 

Proposition 5.2.6. For N < 4 a polynomial P of degree TV with po > 0 is a Hurwith 
polynomial if and only if the following inequalities hold true 

- IfN = 2: Pi > 0,p2 > 0. 

- IfN = 3.- Pi > 0,P3 > 0,P2Pl > P3Po 

- IfN = 4: Pi > 0,P3 > 0,P4 > 0,P3(P2P1 - P3Po) > PiPl 
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Remark 5.2.7. The general conditions for TV > 5 become complicated and are known 
as the Routh-Hurwitz criterion. The criterion consists in N inequalities which can 
be computed either using the diagonal (N-1) dimensional minors of some N x N 
matrix (cf. 113 ) or through a step by step inductive procedure involving only some 
determinants of order 2. 


5.2.3 Remarks on Liapunov’s original proof of the stability theo¬ 
rem 


The original method of Liapunov consisted in introducing the quadratic form 



\\Tit)ufdt 


where T{t) = exp(fA). For a solution of the equation 


It' = Au + F{u) 


we have 


d 

dt 


r+oo 

$(u(f)) = 2 / {T{s)u{t),T(s)u'(t)) ds 

Jo 

P+OO 

= 2 {T{s)u{t),T{s)Au(t) + T{s)F{u(t))) ds. 

Jo 


But 


P+oo P+oo , . 

J {T{s)u{t),T{s)Au{t))ds = J {T{s)u{t),—T{s)u{t))ds = --\\ 


2 


and 


^+oo 


{T{s)u{t),T{s)F{umds < 2C||u(f)|||li^(u(f))|| 


The result then follows for ||F||Lip small enough. On this proof we want to make two 
observations that will justify our choice of a perturbation argument in integral form : 

1) Even when F = 0, the decay rate obtained by Liapunov’s method is not optimal. 
For instance if X = and we apply the above estimates to the equation 


u" + u + 2u' = 0, 


we obtain 

||r(f)|| < 

which is not optimal since in fact 

||r(f)|| < C{1 +t)exp{-t). 

2) When F = 0, the quadratic form $ does not provide the decay in the correct 
space if X is an infinite-dimensional Hilbert space. If, for instance, we consider the 
heat equation 


Ut — Au = 0 in K.’*" x u = 0 on K"*" x dfl 
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in a bounded open domain of which generates a contraction semigroup T{t) on 
X = L^{Tl), the quadratic form $ does not control the norm in X. Indeed , if is an 
eigenfunction of the operator — A , i.e 

—Aipn = XnTn in O, u = 0 on dfl 


it is immediate that 


^+oo ^+oo 1 

Jo Jo 


3) The introduction of $ is only possible when A is a Hilbert space. If, for instance, 
we work with the semilinear equation 

Ut — Au + f{u) =0 in R'*' x H, u = 0 on R“'' x dfi 

and we try to apply Liapunov’s result with X = L^{Tl) , we shall be very limited in 
our range of application. Indeed in order for the operator F defined by 

{F{u)){x) = f{u{x)), a.e. in D 

to satisfy the condition 

||-?^(m)||x < £\\u\\x for ||m||x small 
it is necessary (and sufficient, of course) that / satisfy the global condition 

|/(s)|<e|s|, Vs e R. 

As a consequence, F cannot be tangent to 0 at the origin, except if f = 0. The situation 
is very different if A = C'o(H); in this case, in order for the operator F to satisfy the 
condition 

||-F(m)||x < s\\u\\x for ||m||x small 
it is sufficient that f satisfy the local condition 

|/(s)| < e|s|, for all s small enough. 

In particular, if / is a function of class and /'(O) = 0, F is tangent to 0 at the origin. 
Considering for instance the equation 

Ut — Alt = \u\'^~^u in R“'' x H, u = 0 on R+ x dO. 

The original Liapunov technique does not give the stability of the 0 solution when 
working in L^(H). The method will work if we replace by some Sobolev space 

of type H™‘{Tl), but then we need some growth conditions on the nonlinearity, impos¬ 
ing extraneous limitations on p. If A = C'o(H), we obtain easily the stability of the 0 
solution for any p > 1, cf. Pror)osition l5.3.1l 
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5.3 Exponentially damped systems governed by PDE 

5.3.1 Simple applications 

In this paragraph, we show how the stability theorem [5.1.21 can be applied to partial 
differential equations. 

a) We first consider the semilinear heat equation (12.1b ; 

Ut — Am + /(m) = 0 in R'*" x 12, m = 0 on x dVt 

where 12 be any open set in M.^ with Lipschitz continuous boundary 912 , and / : 
R —> R is a function of class with 

/(O) = 0 and /'(O) > -Ai(12). 

We have the following simple result: 

Proposition 5.3.1. Under the above hypotheses, the stationary solution u = 0 of (12.1b 
is exponentially stable in X = (7°(12). More precisely : for each 7 G (0, Ai(12) + 
/'(O)), there exists R = Riy) such that for all x G X with ||a;|| < R, the solution u of 
(ED such that m(0) = x is global and satisfies 

V2 > 0, ||M(i)|| < 

with M independent of y and x. 

Proof. We have shown in Theorem l3.2.2l that the contraction semi-group T^if) gener¬ 
ated in C°(12) by the equation 

Ut — Xu = 0 in R"*" X 12, m = 0 on R'*' x 912 

satisfies (15.2b with 6 = Ai(12) and some M > 1. It is therefore sufficient to apply 
Theorem l5. 1 .21 with T{t) = To(f), since for / S (7^(R),F(m) = f{u) — 

/'(0)m satisfies (15.3b with a = 0 and p arbitrarily small. □ 

b) Similarly we can consider the semilinear wave equation (12.4b 

utt — Xu + jut + f{u) = 0 in R"*" x 12, m = 0 on R'*' x 912 

where 12 is a bounded open set in R^ with Lipschitz continuous boundary 912, and / 
is a function of class C^\ R —R with 

/(0) = 0 and /'(O) >-Ai(12). 

satisfying the growth condition (12.5b . We obtain the following result: 

Proposition 5.3.2. Under the above hypotheses, the stationary solution (m, v) = (0,0) 
o/(E4ll is exponentially stable in X = TTg (12) x (12) in the following sense: for each 
(5 > 0 small enough, there exists R = R{S) such that for all x G X with ||a;|| < R, the 
solution u of (EH such that (m(0), Ut{0)) = x is global and satisfies 

V 2 > 0, ||m(2)|| < M{S)\\x\\e~^*^ . 


(5.8) 
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Proof. It follows from Proposition 13 .3.1 1 that the contraction semi-group To{t) gener¬ 
ated in X = Hq{D) X L^(n) by the equation 

utt — Au -I- f'{0)u + ^ut = 0 in R+ x fl, u = 0 on R+ x dO. (5.9) 

satisfies (15.2b with some M > 1 for any 5 > 0 small enough. In order to apply 
Theorem 15.1.21 with T{t) the semi-group generated by (15.9b . all we need to check is 
that the function F(u,v)) = —(0,/(u) — f'{0)u) satisfies (15.3b with a = 0 and rj 
arbitrarily small. But this is immediate since the function ip{s) = /(s) — f'{0)s is 
o(|s|) near the origin and, by (12.5b we have |(/3(s)| < (7(|s|’') for s large. Therefore for 
each d > 0 arbitrarily small, we have |(/3(s)| < (i|s| + (7((i)|s|’', globally on R. The 
result then follows immediately from Sobolev imbedding theorems. □ 

5.3.2 Exponentially stable positive solutions of a heat equation 

In this paragraph, we consider the semilinear heat equation 


Ut — Am + f{u) = 0 in R"*" x Tt, m = 0 on R^ x dTt 


where Tt be any open set in R^ with Lipschitz continuous boundary dO , and / is a 
function of class C^: R —>• R with / convex on R+, /(O) = 0 and 


f'iO) < -Ai(fl). 


We have the following simple result: 

Proposition 5.3.3. Under the above conditions, assuming that f{s) > 0 for some 
s > 0, there exists a unique solution ip > Q of 


Ap + f{p) = 0 in El, p = 0 on dO. 


(5.10) 


In addition, p is asymptotically (even exponentially) stable in C{El) n Hq{EI). 

Proof. If a G l°°(El) we denote by Ai(—A + al) the first eigenvalue of — A + al in 
the sense of Hq (17). First of all if (15.10b has a positive solution p and we set 



p(x) - i^ix) 


if p{x) ^ "fix) 


q(x) = f'ipix)) if p(x) = ip(x). 


By strict convexity we have 


q{x) > p{x) 


everywhere in El. In particular 


\x{-A + qI) > 0 
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On the other hand if ip ^ ip, then ip — ip is an eigenfunction of (—A + ql) with 
eigenvalue 0. This contradiction means that tp = ip and therefore (p is unique. In 
addition since by strict convexity we have 

f'{p>{x)) >p{x) 

everywhere in fl we have in particular 

Ai(-A + f{(p{x))I) > 0 

as soon as a positive solution (p exists. Therefore we have uniqueness and exponential 
stability of ip as soon as it exists. 

To prove the existence of (p, first we deduce from the hypotheses on / that 
3so > 0, /'(s) > /'(so) >0, Vs > So. 


In particular 

lim /(s) = lim F(s) = 

s^+oo s—>-+oo 

where 

F{s) = [ f{a)da. 

Jo 

Therefore 


inf F(s) = C > —oo. 

s>0 


For the proof of existence, first we modify (if necessary) / on R by setting 

Vs < 0, /(-s) = -/(s) 

And then F is extended as the primitive of /. This means 

Vs < 0, F{-s) = F(s) 


We introduce 

m = inf{ f [^|Vzp + F{z)]dx, z G iTo(fl) } > Cini > —oo. 

Jn 2 

Since as s —>■ 0 we have 

F{s) - -/'(O)^ 

and /'(O) < — Ai(n), by taking z = epi and letting e —^ 0 we find 

m < 0 

Since any minimizing sequence is bounded in Hq (17) and F is convex up to a quadratic 
term, there exists, as a consequence of compactness in L‘^{El) and Fatou’s Lemma, a 
function p G -^0 (^^) such that 

[ + Fi‘P)]dx = m 

j ri 

Setting Ip = \p\ 'Nt also have, since F is even: 


-\Vip\^ + F{ip)]dx = m 


Because m < 0, of course ip ^ Q. ll is then classical to conclude that ip is a positive 
solution of (15.101 1. □ 
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5.4 Linear instability and Bellman’s approach 

In any finite dimensional real Hilbert space X , the hypothesis of Theorem 15.2.11 is 
sharp. Actually if / = L is linear and has an eigenvalue s := si + is 2 with Si, S 2 real 
and Si > 0, let 

L{ifi + i(p2) = s{(pi + iip2) 

with (fii, (p 2 real vectors and (tpi,ip 2 ) ^ (0, 0). Then the real vector-valued function 
u{t) = e®^*[cos(s2f)(/3i — sint{s 2 t)(f 2 ] 
is a solution of (15.61 1 because the function 


z{t) = 


is a solution of the extended equation of (15.6b on the complexification of X and L being 
a real endomorphism on X, z{t) = e‘‘'^Tp and u = ^{z{t) — z(t)) are solutions of the 
same equation. But now we observe that 

.kir. . ./cTTSi. , , 

u {—) = {—lYexp{ - )t‘(0) 

S2 S2 

and therefore u cannot converge to anything at all as t goes to infinity. 

In the next paragraph we collect some instability results proved in m in the Hilbert 
space framework. 

5.4.1 The finite dimensional case 

Let X be a finite dimensional normed space, and / £ C^{X, X) a vector field on X. 
Let a £ X be such that /(a) = 0. By Liapunov’s theorem (Theorem 15.2.lb . if all 
eigenvalues of Df[a) have negative real parts, a is an asymptotically Liapunov stable 
equilibrium solution of the equation 

u'= f{u{t)), t>0. (5.12) 

This result is sharp since in the opposite direction we have 

Theorem 5.4.1. (R. Bellman, iHK) Let a G X be such that f{a) = 0 and assume 
that at least one eigenvalue of Df (a) has a positive real part. Then a is an unstable 
equilibrium solution of (15.12b . 

Proof. Let p > 0 be the minimum of real parts of the eigenvalues of Df{a) having 
a positive real part and choose an integer K to be fixed later. The Jordan reduction 
theorem implies in particular the existence of an upper triangular matrix T with zero 
diagonal terms and coefficients all equal to 0 or 1 such that 

K 

—M=D+T 

V 

where M is the matrix of Df{a) in a certain basis of X and D is a complex diagonal 
matrix. Then 


M = L + R 
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where L = is a diagonal matrix and R = is a matrix with all coefficients 
having moduli smaller than . Let us identify X with H = with the usual Hilbert 
norm and the associated real inner product. It is clear that the coefficients of L, in 
other terms the diagonal coefficients of M, are in fact the eigenvalues of in Df{a). In 

addition under this identification we have ||i?|| < - —rr—. Let P : H —> H denote 


the projection operator on 


K 


Y := 0 Ker{L-XI). (5.13) 

Re{X)>0 

If u is any bounded solution of (I5.121 i. for t > 0 setting u = a + v we have 

j^i\Pv\^ - 1(7 - P)vr) = 2[iPv, v') - ((/ - P)v, u')] 

= 2[{Pv, Lv + Rv + g{v)) — ((/ — P)v, Lv + Rv + ^(u))] 
where g{v) = f{a + u) — /(a) — Df{a)v satisfies 

g{v) = o(v) 

Since L < 0 on [L]-*- = (/ — P)H we have: 

-((/ - P)v, Lv) = -(L(I - P)v, (I - P)v) > 0. 

On the other hand by definition of rj we have 

Vw G Y, {Lw, w) > r]\w\'^ 


In particular we find 


2{Pv,Lv) = 2{LPv,Pv) > 2r]\Pv\^. 


And therefore 

2[{Pv, Lv) - ((/ - P)v, Lv)] > 2r]\Pv\^. 

On the other hand we have the easy inequality 

2[{Pv,Rv) - {{I-P)v,Rv)] > -4||i?|||z;|2 > 

and since g{v) = o(v), there exists e > 0 such that if |u| < e, we have 

2{Pv,g{v)) - 2{{I - P)v,g{v)) > + ]{! - P)v\'^) 

Choosing K = 8 dim X and combining the above inequalities we find 

- 1(7 - P)v\^) > v{\Pv\^ -]{!- P)v\^) 

whenever |u| < e. Now assuming that a is Liapunov-stable in X, let us select u(0) = 
vq G X such that 


iPuol > |(7-P)uo| 


(5.14) 
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and |tio| small enough so that 


Vt > 0, |t;(i)| < e. (5.15) 

For instance, vq might be any ’’small” vector of Y. As a consequence of the above 
computation it follows that 

Vf>0, {\Pvit)\^-\iI-P)v\^)>e^\\Pvo\^-\iI-P)vo\^)). (5.16) 

This is clearly absurd since (15.141) . (15.151) and (15.161) are incompatible. The proof of 
Theorem l5.4.1l is complete. □ 


5.4.2 An abstract instability result 

The main result of this Section is a natural inhnite-dimensional extension of Theorem 
15.4. ll to the special case of sel-adjoint linearized operator. 


Theorem 5.4.2. Let H be a real Hilbert space with inner product and norm respec¬ 
tively denoted by (•, •) and \ ■ \, L a (possibly unbounded) self-adjoint operator such 
that 


3c > 0, L + c/ > 0. 


(L + (c + 1)/) ^ is compact 


ML) 


i„t 

ugh,u^o I^uI 


< 0 . 


(5.17) 


Assume that there exists a Banach space X G H with continuous imbedding with norm 
denoted by || .|| for which f : X —> H is a locally Lipschitz map with /(O) = 0 and 
such that 


lim 

uGX\{0}, ||u|| 


l/(w)l 


= 0 . 


(5.18) 


Then if X contains all eigenvectors ofL, the stationary solution 0 of 


u' + L{u) = f(^u) 


(5.19) 


is unstable in X. 

Proof. Let P : H —> H denote the projection operator on 

^ Ker{L - \I). 

A<0 

As a consequence of (15.171 ) we know that dim{H~ ) < oo. If m is any bounded solution 
of (|5T^ . for t > 0 u is differentiable with values in H and we have 

±{\Pu\^-\[I-P)u\^) 

= 2[{Pu,u') - {{I - P)u,u')] 

= 2[{PuJ{u)-Lu) + 2{{I-P)u,Lu-f{u))]. (5.20) 

Since L > 0 on [H~]^ = (/ — P)H we have: 


((/ - P)u, Lu) = {L{I - P)u, {I - P)u) > 0. 


(5.21) 
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On the other hand by (I5.171 i we know that 

3?7 > 0,Vw G H~, (—Lw,w) > r]\w\‘^ 

In particular we find 

2{Pu,-Lu) = 2{-LPu, Pu) > 2 p|Pup. (5.22) 

As a consequence of (15.181) . there exists £ > 0 such that if ||u|| < £, we have 
2(Pu, f{u)) + 2((/ - P)u, f{u)) > -plup = -r?(|Pup + |(/ - P)u\^). (5.23) 
Combining (15.201 ). (15.211 ). (15.221 ) and (15.231) we find 

_ 1(7 _ P)u\^) > _ 1(7 _ p)„|2) (5 24) 

whenever ||it|| < e. Now assuming that 0 is Liapunov-stable in X, let us select u(0) = 
uq G X such that 

\Puo\ > |(/ - P)uo\ 

and ||ito|| small enough so that 

Vf > 0, ||u(f)|| < £. (5.25) 

As a consequence of (15.241) we obtain as previously (15.161 ). clearly incompatible with 
(15.25b . Consequently if X contains all eigenvectors of L, the choice 

uq = rjip; Lip = Xip, A<0 hlHi/jH0 

shows by contradiction that 0 is not Liapunov-stable in V. The proof of Theorem l5.4.21 
is complete. □ 

Remark 5.4.3. One might wander why the condition —)■ 0 is required as u —>■ 0 

|m| 

in the sense of X instead of H. Let us consider the example H = L'^{Lt) where fl is a 
bounded open subset of and 

G n : Vu G L^(0), {f{u)){x) = g{u{x)) a.e. in O. 


In this case, the condition 


l/(«)l 


0 as ImI —!► 0 


implies / = 0. Indeed if /(O) = 0 and /(c) 7 ^ 0 we can consider Uu, = cxu> with w an 
arbitrary open subset of LI, so that 

\u^\ = \c\\uj\^; \f{u^)\ = \f{c)\\uj\i 

If |a;| —?► 0 we have by construction 0 and therefore 

NKO \u\ \c\ 

On the other hand if X C L°° with continuous imbedding, the condition 

itoi„t!tWl = o 

ll«IKo |m| 

is equivalent to the natural assumption lim = 0 . 

s-J-O 
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Remark 5.4.4. The instability result in X is only of interest when the existence of 
at least local (and preferably global) solutions for small initial data in X is fulfilled. 
Otherwise Theorem l5.4.2l might just mean failure of existence in X. 

Remark 5.4.5. The proof of Theorem l5.4.2l actuallv implies a stronger instability prop¬ 
erty, namely 

eigenvector of L, 3e„ —>• 0 : sup ||u„(f)|| > a > 0 

t>o 

where Un is the sequence of solutions of ( 15.191 ) such that m„(0) = EnP- This appears 
much stronger since tends to zero in any reasonable norm while the norm of X 
just needs to fulfill (15.181) . 

5.4.3 Application to the one-dimensional heat equation 

Consider the one - dimensional semilinear heat equation 

ut — Uxx + f{u) = 0 in R'*' X (0, L); u{t, 0) = u{t, L) = 0 on R'*' (5.26) 

where / is a function: R —R. Any solution u of this problem which is global and 
uniformly bounded on R+ x (0, L) converges as f —)► +oo to a solution (p of the elliptic 
problem 

ip€H^iO,L), -p^^ + f(p)=0. (5.27) 

Proposition 5.4.6. If ip is a solution of (15.271) which is stable as a solution of (15.261) . 
then p has a constant sign on (0, L) =: Q. 

Proof. Indeed, if p is not identically 0 and vanishes somewhere in (0, L), the function 
w := p' has two zeroes in (0, L) and satisfies 

w G n 77o(0,L), -Wxx + f {t)w = Q in ( 0 ,L). 

In particular if0<Q;<^<Lare such that w{a) = w(/3) = 0, tu ^ 0 on 
{a, 13) and if we set uj = {a, 13), we clearly have Ai(a;; —A + f'{p)I) — 0 where 
Ai(a;; — A + f'{p)I) denotes the first eigenvalue of —A + f'{p)I in the sense of 
Hq{uj). We introduce the quadratic form J and the real number rj defined by 

Wz G Hg{n),J{z) := [ {\Zx\'^ + f'{p)z'^}dx 

Jn 

r] = Inf {J{z),z G Hq{TI), ( z^dx = 1} 

JQ 

Let us also denote by f the extension of w by 0 outside oj. Because 

^(0= f {\Cx\^ + f{p)e}dx= [ {\Cx\^ + f{p)e}dx= [ {\Zx\^ + f{p)z^}dx = 

«/ $*2 J UJ J UJ 

we clearly have 

77 = Ai(0;-A + /'((^)J) <0. 

Assuming 77 = 0 means that a multiple A^ = of C, realizes the minimum of J and 
therefore is a solution of 

V' G C2([0, L]) n iTi(0, L), -fx. + f{<p)f^ = 0. 


This is impossible since if is not identically 0 and however vanishes throughout (0, a) 
for instance. Therefore 77 < 0, and p is unstable. □ 
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5.5 Other infinite-dimensional systems 

The following generalization of theorems 15.4. Il and l5.4.2l is not difficult. 

Theorem 5.5.1. Let H be a real Hilbert space with inner product and norm respec¬ 
tively denoted by (.,.) and |.|, La (possibly unbounded)linear operator such that 

3c > 0, L + c/ > 0 

R{Lt{c+1)I)=H. 

Assume in addition that we have a decomposition H = X (BY with dim{X) < oo and 

X C D{L), LX CX; Y = X-^, L{Y n D{L)) CY, L>0onY. 

Let f : H —> H be a locally Lipschitz map such that /(O) = 0 and such that there 
exists a Banach space V <Z H with continuous imbedding with norm denoted by ||.|| 
for which 

\f{n)\ 


lim 


= 0 . 


uGmol.lkiKo \u\ 

Then ifV contains all eigenvectors of L, the stationary solution 0 of 

u' + Lu = f(u) 

is unstable in V as soon as L has at least one eigenvalue with negative real part and 
eigenvector in X. 

As a typical application of Theorem l5.5.1 1 we can consider the abstract second order 
evolution equation 

u" -Bu -B Au = f{u) (5.28) 

where A is a self-adjoint operator with compact resolvant on a real Hilbert space H 
such that A -f ml > 0 for some m > 0 . Introducing V = D{{A -B {m -B 1)1)^) we 
can set 

% = ¥ xH, D{L) = D(A) X V 

and 

VU = (u, v) € D{L), L{u, v) = (—u, Au -B v) 

Then (15.28b takes the form 

U' + LU = F(u) = {0,f{u)) 

By considering {A„}„gN the nondecreasing sequence of eigenvalues of A eventually 
repeated according to their multiplicity and observing that 


H = y xif= 0[(A-A„)-i(O)]2 


riGN 


it is not difficult to check the hypotheses of Theorem l5.5.1I Hence we can state 

Corollary 5.5.2. Under the above conditions, if A has a negative eigenvalue, and if 
f, W are such that 

\fin)\ 


lim 


= 0 


«GW\{0}.||«||w-*-0 \U\ 

the solution (0, 0) is unstable in the sense ofV :=W x H as a solution of (15.28b . 
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By the same method as in section l5'.4.2l we deduce easily the following consequences 
of Corollarv l5.5.2l : 

Corollary 5.5.3. Let 12 be as in the introduction, / € and ip € C(fi) fl Hq(D) 

a solution of the elliptic problem 

—Ap + f{p) = 0 in Tt; p = 0 on dVl 

such that 

M-A + f'{p)I)<G 

then (^5,0) is unstable in fl Hq(D)] x L^[fX) as a solution of the hyperbolic 

problem 

Utt + Ut — Au + f{u) = 0 in x 12; u = 0 on x 912. 

Corollary 5.5.4. Consider the one - dimensional semilinear wave equation 
Utt + Ut — Uxx + f{u) = 0 in R"*" x (0, L); u{t, 0) = u{t, L) = 0 on R“'' (5.29) 
where f is a function: M ^ R. If p is a solution of the elliptic problem 

pGH^{0,L), -pxx + f{p) =0 

such that {p, 0) is stable in Hq{0, L) x L^(0, L) as a solution of (15.291 1, then p has a 
constant sign on (0, L). 








Chapter 6 


Gradient-like systems 


6.1 A simple general property 

Let S{t) be a dynamical system on {Z, d). We denote by the set of equilibrium point 

of S(t) i.e. 

J- = {x € Z, Vf > 0, S{t)x = x}. (6.1) 

Theorem 6.1.1. Let uq G Z be such that the trajectory S(t)uo has precompact range 
in Z. The following properties are equivalent 

oj{uq) C T, (6.2) 

V/i > 0, d{S(t + h)uo, S{t)uo )—^0 as f —>■+oo, (6.3) 

3a > 0, V/i G [0, a], d{S{t + h)uo, S{t)uo )—^0 as f —>■+oo. (6.4) 

Proof, i) (I6.4l i implies (I6.2l i. Indeed assume (16.41) and let x G oj(uo). There exists 
tending to +oo for which 

lim S{tn)uo = X. 
n—>-oo 

Therefore by continuity of S{h) 

V/i > 0, lim S{tn + h)uo= lim S{h + tn)uo = S{h)x. 

n—ftyo n—ftyo 

As a consequence of (16.41 ) we have on the other hand 

V/i G [0,a], lim S{tn + h)uo = x. 

n—>-oo 

By comparing the two previous formulas we find 

V/iG[ 0, a], S{h)x = X. 

Then a trivial induction argument gives 

V/iG[0,a], Vn G N, S{na + h)x = x. 

This obviously implies (16.21) . 

ii) (16.2b implies (16.3b . Indeed assume that (16.3b false. Then for some h > 0 there is 
an e > 0 and a sequence tending to +oo for which 

Vn G N, d{S{tn +h)uo,S{tn)uo) > e. 
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We can replace the sequence by a subsequence, still denoted tn, for which S{tn)uQ 
converges to a limit x G X. As a consequence of (16.21 1 we have x G X. By letting 
n tend to infinity in the above inequality, since S'(f„ + ft.) = S{h)S{tn) and S{h) is 
continuous we obtain 

d{S{h)x,x) > e. 

This contradicts (I6.21 i. Hence (I6.21 i implies (16.31) and this concludes the proof. □ 


6.2 A minimal differential criterion 

In this section we assume that Z is a closed subset of some Banach space X. 

Corollary 6.2.1. Let uq G Z be such that the trajectory S{t)uo has precompact range 
in Z. Assume in addition that 

S{t)uo =: u{t) G 

Then if 

t-t+ 

3a > 0, J ||M 

we have dOi . 

Proof. It is sufficient to observe that 

rt-\-cx 

M'(s)fis|| < j ||u'(s)||ds —>■ 0 as i —>•+oo. 

Hence (16.41) is fulfilled, and by Theorem l6.1.1l this implies (16.2b . □ 

Corollary 6.2.2. Let uq G X be such that the trajectory S{t)uo has precompact range 
in Z. Assume in addition that 


'(s)||(is —>■ 0 as t —>■+00 


(6.5) 


S{t)u^=-.u{t)GWtl(M+,X). 


Then if for some p > 1 


u' G LP(M+,X) 


(6.6) 


we have (IQ) . 

Proof. Indeed in this case we have 


(■t+i 




r-t+l 


.'(s)||P(is) —!► 0 as f —>■+ 00 . 


□ 


6.3 The case of gradient systems 

Let > 1 and F G We consider the equation 

u'{t) + XF{u{t)) = 0 (6.7) 

and we define 

£’ = {2GR",VJ^(z)=0}. 
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Corollary 6.3.1. Any solution u{t) of (16.71 defined and bounded on R'*" satisfies 

lim dist{u(t),£} — 0. 

In other terms we have a;(w(0)) C £. In addition, if for each c, the set £c = {u £ 
£, F(u) = c} is discrete, then there exists u* G £ such that 

lim u(t) = u*. 

Proof. We consider the dynamical system generated by (lO) on the closure of the 
range of u. It is obvious here that the set T of fixed points of S (t) is precisely equal 
to £ defined above. Multiplying by u' in the sense of the inner product of and 
integrating we find 

\\u'{t)\\^dt = F{u{0) — F{u{t). 

Hence since u is bounded we obtain u' G L^(R~^,X) with X = . By Corollary 

16.2.21 we have a;(M(0)) C £. Moreover F{u{t)) is non-increasing along the trajectory 
since 

j^F{u{t)) = -\\u\t)r 

Hence F{u{t)) tends to a limit c as f becomes infinite and therefore w(rt(0)) C £c- The 
rest is a consequence of Proposition l4.2.1l □ 

Remark 6.3.2. By using Lemma fl .2.21 applied to the function ||u'(f)|p it is easy to 
prove that u'{t) tends to 0 at infinity. One might wonder whether u{t) is always con¬ 
vergent. In 2 dimensions, it was conjectured by H.B. Curry 041 and proven by J. Palis 
and W. De Melo llT^ that convergence may fail even for a C°° potential F. 



6.4 A class of second order systems 


Let F, £ be as in Section l63] We consider the equation 


u"{f) E u (f) -f XF{u{t)) = 0. 


( 6 . 8 ) 


Corollary 6.4.1. Any solution u{t) of (16.SI defined and bounded on K+ together with 
u' satisfies 


lim |lu'(f)|| 

t^+OO 


lim dist{u{t),£\ = 0 

t—>-+CXD 


In other terms we have w(m( 0), u'(0)) G £ x {0}. In addition, if for each c, the set 
£c = {u G £, F{u) = c} is discrete, then there exists u* G £ such that 


lim u(t) = u*. 

Proof. We consider the dynamical system generated by (I 6 . 8 I 1 on the closure of the 
range of U = {u, u'). Here the set F of fixed points of S{f) is made of points (y, z) G 
X for which the solution u of (16.81) of initial data {y, z) is independent of t. 
Consequently F = £ x {ff}. Multiplying by u' in the sense of the inner product of R^ 
and integrating we find 

J^{\\w{t)r+F{um = -\w{t)r 
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hence in particular 

£ \\u\t)rdt = F(u(0) - F{uit) + i(||u'(0)f - llu'Wf). 

Hence since u is bounded we obtain u' G L^(R+,X) with X = . Moreover 

differentiating the equation we have 

u'" + u" + X^F{u{t))u' = 0. 


By multiplying by u" in the sense of the inner product of and integrating we find 



W'itWdt 



{V^F{u{t))u\u"{t))dt + ^{\\u"{0)r 




Since u" is bounded by the equation, it follows immediately that u” G L^(R^,X), 
therefore 17' = (u',u") G XxX).Bv Corollarv l6.2.2l we have oj(u(0). u'fO)) C 

S X {0}. In particular u'(t) tends to 0 as f becomes infinite. Moreover i||u'(f)||^ + 
F(u(t)) is non-increasing along the trajectory and therefore tends to a limit c as t 
becomes infinite. Finally a;(u(0), u'(0)) C Sc x {0}. The rest is a consequence of 
ProDosition l4.2.1l □ 


6.5 Application to the semi-linear heat equation 

Let fl and / be as in Section 12741 and let X = Co(D). Throughout this section we 
assume that fl is bounded and we define 

S = {uGXn i7d(H), - Am -f f(u) = 0}, 

Vip G X n Hg, F(ip) = - I \Vip\'^dx+ j F{ip)dx 
2 Jn Jn 

with ^ 

F{u) = f f{s) ds, Vu G R. 

Jo 

Moreover let Sc = {u G S, E{u) = c}, for c G R. We shall prove 

Theorem 6.5.1. let u be a global solution of (12. Il l which is bounded in X for t > 0. 
Then we have the following properties 

(i) E(u{t)) tends to a limit cas t ^ +c»; 

(ii) Sc f 0; 

(iii) dist(u{f), Sc) 0 as t ^ +oo, where dist denotes the distance in X H Hq{D). 

Proof. The smoothing effect of the heat equation implies (cf. e.g. Il60l for a proof 
based on the theory of holomorphic semi-groups) that for each e > 0 and a G [0,1), 

U{ u{t)} is bounded in 

t>e 

In particular, Ut>o{^(^)} precompact in X and Ut>i{r*(f)} is precompact in Hq{D). 
Let us denote by Z the closure in Xni7d(^) of m(R+). E is continuous on Ani7d(f2), 
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hence on {Z, d) where d is the distance in X fi H^. In addition by precompactness the 
topologies of X n Hq{LI) and Ld{Ll) coincide on Z. An easy calculation shows that 
for t> 1, we have 


f j uf{T,x)dxdT + E{u{t)) = E{u{l)) 

Ji Jn 


Hence By Corollary 16.2.21 we have uj{u{0)) C £. Since E{u{t) is nonincreasing the 
result follows as in the previous examples. □ 


6.6 Application to a semilinear wave equation with a 
linear damping 


Let Lt and / be as in Section 1231 and consider the wave equation (12.41) . Throughout 
this section we assume that SI is bounded. Keeping the notation and the hypotheses of 
Section lZSl Introducing 


£ = {u e Hq,—Au + f{u) = 0}, 
and £c = {u G £, E(u, 0) = c}, for c € R. We can state 

Theorem 6.6.1. Assume 7 > 0 and that the growth condition (lOT l is satisfied with the 
strict inequality: r < 2/(iV — 2) if N > 3. Let (ip, fi) G X := Hq x L^, and let u be 
the corresponding maximal solution of (1^ with u(0) = ip and Ut{ 0 ) = ip. Assume 
that T{ip, Ip) = +CX) and 

SVLp{\\{u{t),Ut{t))\\xfi> 0 } < +00. 

Then we have the following properties : 

(i) E{u{t), Utif)) tends to a limit cas t ^ +00; 

(ii) £c f 0; 

(iii) ||Mt(f)||L2 —>• 0, as f —>■+00; 

(iv) dist{u{t), £c) —^ 0 as f > +c», where dist denotes the distance in H^. 

The proof of Theorem 16.6.11 relies on a general compactness lemma due to G.F. 
Webb ia : 

Lemma 6.6.2. Let X be a real Banach space and T(t) a contraction semi-group on 
X satisfying 

\\TmLix)<Me-'^\ Vf>0. (6.9) 

where M,a are some positive constants. Let H G L'''°“(R“'', 26 ) and let K be a 
compact set in X such that H{t) G K , a.e. on R'*’. Then the function V : ^ X 

defined by 

V(t) =T{t){ip,ip)f T{s)H{t—s)ds 
Jo 

satisfies: V (R"*") is precompact in X. 
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Proof. We have V{t) = where 


W{t)= [ T{s)H{t-s)ds. 
Jo 


Since T{t){(p,ip) —0 in X as t —>■ +oo, there is a compact subset Ki of X such 
that V')} It is therefore sufficient to prove that there is a compact 

subset K 2 of X for which 

[j{Wit)} c K 2 . 

t>0 

Let e > 0, and according to (I6.9L let r be such that 

/ oo 

II'?"('S)IIl(x)£Is < 


For f > T, we have 


consequently, 


\\W{t)- f T{s)Hit-s)ds\\x <s 
Jo 

\J{W{t)}cK 3 + B{0,e) 


( 6 . 10 ) 


t>T 


with 


K. 


y{ f\is)Hit-s)ds}. 

t>r -^0 


Observe that the map (s, x) 1 —>■ T{s)x is continuous: [0, +c») x X —>■ X. As a conse¬ 
quence, U = y T{t)K is compact in X. Hence, F = t ■ conv([/) is precompact 

0<t<T 

in X. Since K'i C F. K-i is nrecomnact in X. Bv (I 6 .IOI 1 . we can cover y{w"W} 

t>r 

by a finite union of balls of radius 2e. On the other hand, W G <^([0, -foo), X), hence 
y {fF(f)} is compact and can also be covered by a finite union of balls of radius 

0<t<T 

2e. Finally we can cover y {fF(f)} by a finite union of balls of radius 2e. Since e > 0 

t>o 

is arbitrary, y {H^(f)} is precompact, and the conclusion follows. □ 

t>o 

Proof of Theorem l6.6.1l We define an unbounded operator AC on X by 

D{AF) = {{u, v) G X, Au G and v G Hq}; 


A'^{u, v) = {v, Au — 'yv),y{u, v) G D{A'^). 

It is easily seen that A'*' is m-dissipative on X. Asa consequence of Proposition B. 3 .T1 
the contraction semi-group T{t) generated by AC on X satisfies (16.91 1. 

Now set U{t) = {u{t),ut{t)) and H{t) = (0, —/(^(t)), for f > 0. Clearly rt is a 
solution of (I 2 . 4 I 1 on [0, r] if, and only if U G C([0, r]; X) and ?7 is a solution of the 
equation 

U{t) = + f T{t — s)H{s)ds, VfG[0,T]. 

Jo 
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Now we recall the energy identity 


7 / f u^{t,x)dxdt + E{u{t),ut{t)) = E(ip,'ip) 

Jo Jci 


with 


■= \ ( +i f \Jj{x)\^dx+ j F{if{x))dx. 

E is continuous on X, hence on (Z, d) where Z is the closure of u(R+) in X. The 
energy identity shows that E{u{t),ut{t)) is non-increasing. The set of stationary points 
of the dynamical system is easily identified as £ x {0}. On the other hand the function 
H : K.+ —^ X defined by H{t) — (0,—/(«(<)) for f > 0 is such that iJ(K.+ ) is 
precompact in X. (This comes from the strict condition; r < 2/ (TV — 2) if > 3.) 
Applying Lemma l6.6.2l we obtain compactness of bounded trajectories in X. Then the 
topologies of = Hq X and Y = LJ x H~^ coincide on Z and an easy calculation 
using the equation now shows that 

u' = {uuutt) e L"(]R+,r). 

Indeed the energy identity gives Ut S L^). On the other hand the growth as¬ 

sumption on / is easily seen to imply 

V(u,u)G2f, f'{u)v€H~^ 


with 

\\nu)v\\H-^ < Kii + \\ur^,j\\vh^. 

By multiplying the equation by utt in the sense of H~^ and integrating in t we find 

/ + [ uutdx 

Jo ^ J 0 

+ [{fiu),Ut)H-Ao+ / {-J^u,Utt)H-^ds = / {ut,f{u)ut)H-^ds. 

JO JO 

Hence, using the identity 


/ {Au,utt)H-^ds = / \\Xut\\H-ids + [{Au,ut)H-^] 
Jo Jo 

= [ IWtWlds+ [{Au,ut)H-Ao 

JO 


we derive easily 


< Cl-f 02 / \\ut\\lds. 


Then the conclusion follows as in the previous example. 


□ 


Remark 6.6.3. Under the conditions of Prot)osition l2.5.11 the conclusions of Theorem 
16.6. H are valid for any solution u of (12.4b . 
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Liapunov’s second method and 
the invariance principle 

7.1 Liapunov’s second method 

As explained in Section l5.2.3l the Liapunov stability theorem for equation (15.61) can be 
proved by exhibiting a positive definite quadratic form which decreases exponentially 
along the trajectory if the initial data are close enough to the equilibrium under consid¬ 
eration: such a function is called a Liapunov function. Sometimes it is possible to find 
directly such a function without calculating the fundamental matrix of the linearized 
equation, and this is the principle of the so-called ‘direct” or second Liapunov method. 
This method can often be reduced to the following simple criterion 

Proposition 7.1.1. Let V € be such that V{u) tends to -foo as ||m|| —>■ -l-oo 

and let a G be such that 

Vu^a, {V'iu)Jiu)) <0. (7.1) 

Then we have f{a) = 0 and in addition 

- 'iu G V{u) > V{a) 

- a is an asymptotically stable equilibrium point of the equation u' = f(u). 

Proof. Since V is continuous and V{u) tends to -foo as ||m|| —)► +oo, then there exists 
b G such that V{u) >V(b) for all u G R^. Clearly we have V'{b) = 0 and now 
(ED imply that b = a. 

Once again by (17.11) V is non-increasing along the trajectories, therefore all trajec¬ 
tories are bounded. Given such a trajectory u , let G uj{u{0)) and let z be the solution 
of 

z = fiz) z( 0 ) = p 

Since V{u{t)) tends to a limit Z as f —>■ -|-c», we have by (14.2b 

Vf, V{z{t)) = I 

and then 

Vf, {V'{z{t)),f{zm = jV{z{t))=0. 
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In particular Vt, z{t) = a, hence ip = a and since z is constant we have /(a) = 
f{z{0)) = z'(0) = 0. The stability of a follows easily from Corollary 14.4.31 Indeed 
for any trajectory u we have 


^V(u(t) = {V'(u(t))J(u(t)) 


therefore either u(t) = a or ■^V{u{t) < 
V{a) = liml7(u(f)) < 17(u(0)) . 


0. Whenever u(0) a we deduce that 

□ 


Example 7.1.2. Let us consider the system 


u' = —u 


cv 

1 + a|t;| ’ 


v' = —V 


du 

l + /3|u| 


where a > 0, /3 > 0 and sup{|c|, |d|} < 1 which has the form (15.61) with / Lipschitz 
but not differentiable at the origin except when a = /3 = 0 . Setting 


V{u,v) = 


we find easily that 

V(m, v ) ^ (0,0), (F'(u,u)), /(u,u)) = -2(u^ + 

< -2(1 - sup{|c|, |fi|})(M^ + < 0 

Therefore (0,0) is the unique equilibrium point and is globally asymptotically (here 
exponentially) stable. In the special case c — d > Q and a — fi = 1 , assuming 
Mo > OjUo > 0 the solutions of the above system with initial data (uo,vo) remain 
non-negative for all times and coincide with the solutions of 


/ , 

u =—u-\ -; V =—V 

1 + v 


cu 

1 + u 


which is known as the Naka-Rushton model for neuron dynamics in the short term 
memory framework. 


7.2 Asymptotic stability obtained by Liapunov functions 

Consider the nonlinear wave equation 

utt — Au + g{ut) =0 in R'*' x SI, u = 0 on R'*" x dU (7.2) 

where SI is a bounded domain of R^ and g satisfies the conditions 

3a > 0, g{v)v > a\v\^, Vu S R (7.3) 

3C > 0, |g(t;)| < (7(|u| + |t;|T'), Vu G R (7.4) 

with ; 

7>1 and if iV > 3, 7 < (A^ + 2)/(iV - 2). (7.5) 

For the sake of simplicity we consider classical solutions of (17.21) for which dif¬ 
ferentiations are plainly justified. We obtain the following result of global asymptotic 
stability : 
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Theorem 7.2.1. Let 

u G n Him n (R+, n (R+, 

be a solution of (O- Then we have 

I {\\7uf + u^}{t, x)dx < M ( I \Vu{0, x)\'^dx, f \ut{0, x)\'^dxl (7.6) 

Jn Jn ' 

where 5 > 0 depends only on a, C and 7 and M is bounded on bounded sets. 

Proof. We denote by (•, •) the inner product in if {El), by | • | the corresponding norm 
and by || • || the norm in i7g (li). In addition the duality pairing on H~^{Tl) x Hl{D.) 
will be denoted by (•,•). Now we dehne 

= ||M(i)f + \utm + e{u{t),utit)) 

where e > 0 is at our disposal. For e small enough, is comparable to the usual 
energy and we obtain ; 

d f 

+ lutitm = {utt + Lu,Ut) =-2 / g{ut)utdx 

J Q 

^{u{t),ut{t)) = |zit(<)P + - IhWf - J^g{u')udx. 

Therefore: 

=-2 [ g{ut)utdx + e\utm - e\\u{m - £ [ g{ut)udx. (7.7) 
™ JQ Jn 

It follows from (I7.31 l and (I7.41 i that 


|'?(t^)| < 2C\v\ for lul < 1 

< 2C(vg{v)y for |u| > 1. 

In particular for each v G we have by setting (3 = and denoting as || • ||/3 

the norm in L^{D) 


||5(t>)||/3<2C||u||;3 + (2C)^ 



vg(v)dx 


<Ci\\v\\p+C 2 



1 


Since the condition 7 < (TV + 2)/(A^ — 2) yields /3 = = 

implies 


^ = (-a+ smm - s\\u{t)r+mwurniuM 


- / g{ut)utdx + C 2 £( / utg{ut)dx)] ||u(f)|| 


(2*)', (O 


(7.8) 


By reordering the terms and using Young’s inequality with exponents 7 + 1 and (3 we 
deduce from (I7.81 l : 


^ < (-f + e)\u,m + {Ke^ - e)\\u{tW + {C2ey+^u{t)\m■ 
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Since 2E{t) = ||M(t)||^ + \ut{t)\'^ is a nonincreasing function ot f > 0, we can first 
choose e > 0 small, depending on E{0), such that 

Vf>0, + (7.9) 

This shows that E{t) —^ 0 exponentially, uniformly on bounded subsets of iTg (^) ^ 
L^{Ll). Then for each initial condition , we can find Tq > 0, depending on i7(0), such 
that E{t) < 1 whenever t > Tq. Now for t > Tq, we have ( 17.91 ) with e independent of 
E(0). Hence (17.61 ) follows with S independent of £'(0). □ 


In section 1531 we saw that even in the nonlinear case , the existence of an eigen¬ 
value s of Df{a) with Re{s) > 0 implies the instability of a. On the other hand, in 
the marginal case Re{s) = 0 (for instance when s = 0), a can still be asymptotically 
stable, as shown by the next examples. 

1) A typical example of such a situation is the first order ODE 


u' = —|u|^ ^u, t>0 


with p > 1. The solutions uO of (17.101 ) are given by the formula 


u{t) = 


_ sgn{uo) _ 

{(p- l)i + 


(7.10) 


(7.11) 


It is clear from (17.1 lb that 


as f —>■ +00 for every uq ^ 0. Analogous, but somewhat artificial parabolic example 
can be given. Let us consider now some second order examples. 

2) First we consider the equation ( with c > 0,p > 1.) 

u” + u + c\u\P-\' = 0, t>0. (7.12) 

We set ip{t) = then 

p'{t) = -2c\u'\P+^ > -2c(w" +m'2)(p+i)/ 2 = _2cp(f)(P+i)/2 
and as in the previous example we deduce 


pit) > 


1 

[p(0)]^ +c(p-l)t 


2 

P-1 


Hence the energy tends to 0 at most like t as t —>■ +oo. In fact we have 


Proposition 7.2.2. For each solution u of (17.12b we have 

Vf > 0, {v?{t) + u’'^{t)}^ < C{u{0),u'. 


(7.13) 
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Proof. We set: 

$E(t) = u^{t) + u'^{t) + e\u(t)\^~^u{t)u (t) 

Then: 

= -2c\u'\P+^ + e\u\P-\pu'^ + uu") = -2c\u'\P+^ + e[p\u\P-\'^ - |m|p+i 

-c\u'\P-^u'\u\P-^u) < -2c\u'\P+^ + e{-{l/2)\u\P+^ + C\u'\p+^}, 
where C depends only on u(0), m'( 0). For e > 0 small enough, we therefore obtain 
< -{e/2){\u\P+^ + |u7+i} < 

Clearly, (17.141 1 implies (17.131 1 for e small enough. □ 

3) Finally , by using the method of proof of Theorem l7.2.1l one can prove 


Theorem 7.2.3. Assume that g S C^{R) satisfies the conditions 
3a > 0, g{v)v > Vw € R, 

3C > 0, |(?(t;)| < (C(|t;| + |uP),Vu e R, 
with : 1 < p < 7,7 satisfying (|73T |. Then for each solution 

u G Li:,(R+, 77' n Him n (R+, n (R+, 

0/(O we have 

f {\\7u\^+ u^}{t,x)dx < m( f |Vu(0, a:)p(ia;, f \ut{0,x)\^dx\u^ (7.15) 

Jn ^Jn Jn ' 

Idea of the proof. Let 

= iiu(f)f + \u'm+ 

By adapting the proof of Theorem l7.2.1l and Proposition l7.2.2l we establish 

valid for e > 0 small enough and some <5 > 0 depending on the initial energy. 
Remark 7.2.4. It is not known whether (I7.151 l is optimal when for instance 

g(v) = c\v\P~'^v, c > 0, p > 1. 

A very partial result in this direction (lower estimate comparable to the upper decay 
estimate raised to the power |) can be found in ll50l in the case TV = 1, relying on an 
argument specific to dimension 1. 
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7.3 The Barbashin-Krasovski-LaSalle criterion for asymp¬ 
totic stability 

After Liapunov, the stability theory has been pursued mainly by the russian school 
which was also involved in control theory of ODE under the impulsion of major russian 
experts such as L. S. Pontryaguin. In this context, interesting contacts have been estab¬ 
lished between the russian school and american mathematicians such as J.K. Hale and 
J.P. LaSalle. The exchanges between J.P. LaSalle, E. A. Barbashin and N.N. Krasovskii 
led to the now well-known invariance principle, and LaSalle in his papers is quite clear 
about the influence of the russian school on his own research. To illustrate the progres¬ 
sion of ideas, we start with a simple and convenient result about asymptotic stability. 

Theorem 7.3.1. Let f G and consider the differential system dSHi. Let a G 

be such that f{a) = 0 and U be a bounded open set with a G U such that 

(i) For any x close enough to a, the solution u of (EHi with u{0) = X is global and 
remains in U. 

(ii) 3V G C^(R^) such that 

VuGU, {V'iu)Jiu))<0. 

(Hi) The set u G U, {V'{u), f(u)) = 0 contains the range of no trajectory of (15.6b 
except the constant trajectory a. 

Then a is a strict local minimum ofV, it is the only equilibrium point in U and a is an 
asymptotically stable equilibrium point of (ES. 

Proof. Given a trajectory u of (15.6b with u(0) close enough to a so that u remains in 
U, let (fi G uj(u) and let z be the solution of 

z' = fiz) ^(0) = ip 

Since V{u{t)) tends to a limit Z as f —>■ -foo, we have 

Vf > 0, V{z{t)) = I 

In addition Vf > 0, z(f) G C/and (17'(z(f)),/(z(f))) = ^E(z(f)) = 0. 

In particular as a consequence of (iii) we have Vf > 0, z{t) = a, hence (p = a. So 
u{t) converges to a as f —>■ oo. Moreover if m(0) f a, by (iii) there is some T G R’*' 
for which {V'{u{T)), f(^u{T)) < 0 and then 17 (m(0) > V(^a). Therefore a is a strict 
local minimum of V and the conclusion now follows from Corollary 14.4. 31 □ 

Example 7.3.2. Let us consider the system 

u' = u; v' = —u — g{v) + c 
where c G R and g is increasing with g{0) = 0. Setting 

V{u, v) = {u — c)^ + 

we find easily that V(m,u) G R^, ly' {u,v)), f{u,v)) = —2g{v)v < 0. Taking for 
U any ball centered at (c, 0) , conditions i) and ii) are obviously fulfilled. Then if a 
trajectory {u,v) satisfies {V{u,v)), f{u,v)) = 0 , from —2g{v)v = 0 we deduce 
u = 0, hence u' = 0 and by the second equation u = c. Einally (c, 0) is the only 
equilibrium and is globally asymptotically stable as a consequence of Theorem l7.3.1l 
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Example 7.3.3. Let us consider the system 

u' = v; v' = J~^ {—psinu — kv + c) 

where c > 0 and J, p, k are positive with c < p. This represents the motion of a robot 
arm with one degree of freedom submitted to a constant torque. Setting 

V(u, v) = —v^ + p(l — cos u) — cu 


we find easily that 

V{u,v) G {V'{u,v)), f{u,v)) = —< 0 

We claim that the hypotheses of Theorem l7.3.1l are satisfied when a = argsin ^ and 
a = {a, 0). Indeed from the equation above it follows that the function V{u{t),v{t)) 
is constant if and only if {u{t),v{t)) = (/3,0) and psin/? = c. Moreover, setting 
F{u) = p{l — cosu) — cu we see immediately that 

F'{a) = psina — c = 0, F''{a) = pcosa > 0 

Therefore a = (a, 0) , is a strict minimum of V, and is consequently a stable equi¬ 
librium by Corollary 14.4.31 Since a is an isolated solution of this equation, the only 
possibility is /3 = a. By Theorem 17. 3. II we conclude that a is asymptotically stable. 
The same property holds true for the other equilibria of the form (a + 2fc7r, 0). 


7.4 The general Lasalle’s invariance principle 

Let {Z, d) be a complete metric space and {S(t)}t>o a dynamical system on Z. 

Definition 7.4.1. A function^ G C{Z,]&) is called a Liapunov function for {S{t)'\t>o 
if we have 

$(S'(<)z) < $(z), VzGZ, Vf>0. (7.16) 

Remark 7.4.2. By using the semi-group property of S{t), it is immediate to see that 
$ is a Liapunov function for {S'(f)}t>o if, and only if for each z G Z the function 
1 1 —>■ <l>(S'(f)z) is nonincreasing. 

The following result is known as LaSalle’s invariance principle. 

Theorem 7.4.3. (cf MZl?) Let ^ be a Liapunov function for {<S'(f)}t>o. o-nd let z & Z 
be such that (J {S'(f)z} is precompact in Z. Then 
t>o 

(i) c= lim ^{S{t)z) exists. 

(ii) $(y) = c,VyG uj(z). 

In particular: 

Vy e uj{z),Vt>0, $(S'(f)y) = $(y). 








7.5. APPLICATION TO SOME DIFFERENTIAL SYSTEMS IN 


69 


Proof, (i) ^{S{t)z) is nonincreasing and bounded since is precompact. 

t>o 

This implies the existence of the limit c. 

(ii) If y € w(z), there exists a sequence —>■ +oo such that S{tn)z —>■ y. Hence 

^{S{tn)z) —>■ $(?/) and this implies $(?/) = c. 

The last property is now an immediate consequence of the invariance of uj{z) (the- 
orem l4.1.^ i)). □ 

Remark 7.4.4. Practically, Theorem l7.4.3l is used to show the convergence of some 
trajectories of {5'(f)}t>o to an equilibrium. Therefore the following definition and 
theorem are basic. 

Definition 7.4.5. A Liapunov function $/or {S'(f)}t>o called a strict Liapunov 
function if the following condition is fulfilled : Any z € Z such that $(S'(f)z) = $( 2 :) 
for all t > 0 is an equilibrium of {S{t)}t>Q. 

Theorem 7.4.6. Let ^ be a strict Liapunov function for {S'(f)}t>o, and let z € Z be 
such that (J{5'(f)z} is precompact in Z. Let £ be the set of equilibria of {S{t)}t>o. 
t>o 

Then 

(i) £ is a closed, nonempty subset of Z; 

(ii) d(S(t)z,£) 0 as t ^ +c», i.e. uj{z) C £. 

Proof. By continuity of S(t),£ is closed. By Theorem l4.1.8l ffl. uj{z) f 0. Now let 
y G w(z). The last assertion of Theorem l7.4.3l gives 

^{S{t)y) = ^{y), Vf>0 

and therefore, since T* is a strict Liapunov function, y is an equilibrium : in particular 
we have (i) and a;(z) C £. Then Theorem l4.1 .SK iiil implies (ii). □ 

Remark 7.4.7. Theorem |7A6] means that the set of equilibria attracts all trajectories 
of {S'(f)}t>o. 

Corollary 7.4.8. Under the hypotheses ofTheorem \7.4.6\ let 

c= lim ^{S{t)z) and £c = {x G £,^{x) = c}. 

Then £c is a closed, nonempty subset of Z and d(^S{t)z,£c) 0 as t ^ +00. If in 

addition £(. est discrete, there exists y G £c such that S(f)z ^ y as t ^ +00. 

Proof. Since £ is closed and $ is continuous, £c is closed. The rest of the corollary is 
a consequence of Theorems 17.4.3117.4.61 and 14.1.81 (ii). □ 


7.5 Application to some differential systems in 

Theorem l7.4.3l Theorem l7.4.6l and Corollary 17.4.81 allow to recover easily the results 
of chapter 6 on gradient systems and second-order gradient-like systems with linear 
dissipation. But they show their full power in more complicated situations in which 
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calculations implying convergence to 0 of the time-derivative become less natural. As 
a typical example we can consider the equation 


u''{t) + g{u'{t)) + VF{u{t)) = 0. (7.17) 

where F S R) and g : is a continuous function such that 

Vu G R'^ \ {0}, (p(u),u)>0. 


Corollary 7.5.1. Any solution u(t) of (17.17b defined and bounded on R'*' together with 


u' satisfies 


lim |lu'(f)|| = lim 


dist{u{t),£} = 0 


with 

£ = {z&M.^,VF{z) = Q}. 


If in addition for each c, the set £,. = {u £ £, F(u) = c} is discrete, then there exists 
u* € £ such that 

lim u(t) = u*. 


Proof. We consider the dynamical system generated by (3) on the closure of the range 
of U = (u,u'). Here the set F of fixed points of S[t) is made of points (y, z) G 
R^ X R^ for which the solution u of (17.17b of initial data (y, z) is independent of t. 
Consequently F = £ x {If}. Multiplying by u' in the sense of the inner product of R^ 
and integrating we find 


j^{\\\u'{t)rFF{um = -{g{u'),u')<Q 

hence 

^{u,v) := i||uf + F(u) 

is a Liapunov function. On the other hand if $ is constant on a trajectory {u{t),u'{t)) 
we have u' = 0. Hence $ is a strict Liapunov function and the result follows. □ 

As an example of application of Corollary 17. 5. II the equation 


u” + u' + — u = 0 


already considered in Section 4.2 provides a good illustration. Here the set of equilibria 
has only points solutions : (-1, 0), (0, 0) and (1,0). Note that here and more generally 
under the hypotheses of Corollary 17. 5. 11 the t-derivative of the Liapunov function van¬ 
ishes at some point to only if u'{to) — 0. Then it follows easily that energy conserving 
trajectories are made of equilibria. In the next example the condition u' = 0 does not 
follow immediately, but as a consequence of the connnectedness of trajectories; 

Example 7.5.2. Let us consider the scalar equation 

u" + au^u'+ — u = 0 (7.18) 


m = 



= —av?u'^ 


where a > 0. Let 
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Since E is non-increasing, {u, u') are bounded and we are in a good position to apply 
the invariance principle. Indeed let it be a solution of (17.181) for which E is constant, 
then uu' = 0 hence v? is constant and then, by connectedness, u is constant. So 
u' = It" = 0. As in the previous example,the stationary equation — u = 0 has only 
three solutions : -1,0, 1. So that we have convergence of all solutions, although the 
t-derivative of the Liapunov function vanishes also at points to for which it(fo) = 0 
and the equation is not a special case of Corollary 17. 5. II 


The next example shows that sometimes, the invariance principle provides some 
information which is not so easy to recover by more elementary methods. 

Example 7.5.3. Let us consider the coupled system of second order scalar ODE: 


It" -I- It' -I- Alt -I- ci; = 0, 
v" + Xv + cu = 0, 


A > 0 and c ^ 0 with < A^. Let 

E{t) = £{u, It', V, v') = i [u'^ -I- it'^ -I- + i;^)] + cuv. 


(7.19) 


Since E is non-increasing, (it, u, it', v') are bounded and we are in a good position to 
apply the invariance principle. Indeed let (it, v) be a solution of (17.19b for which E is 
constant. Then 2ii'^ = 0 implies it' = 0, hence it is constant and it" = 0. Then by the 
hrst equation v = —-^it is also constant. Finally since by the hypothesis < A^ , the 
stationary system Ait -b ci; = cit -I- Au = 0 has no non-trivial solution, we conclude 
that u = V = 0 and therefore (0,0, 0,0) is asymptotically stable. Because the system 
is linear and hnite-dimensional, by taking a basis in it follows immediately that the 
norm of the fundamental matrix tends to 0 as f tends to inhnity, and using the semi¬ 
group property it follows that convergence is exponential. The general theory designed 
by Liapunov in his seminal paper (1892) shows the existence of a quadratic form $ on 
satisfying the identity 


(with T] > 0 )for any solution Y = {u, v, it', v') of (17.19b which means that the older 
method of quadratic energies must allow to recover directly that (0, 0,0, 0) is asymp¬ 
totically stable, with quantitative information about the decay rate. Since the form can 
be computed on a basis of = iQ monomials in (it, v, it', !>'), the challenge is 

now to find one of the strict quadratic Liapunov functions (they form a non-empty open 
set in the space of coefficients) by a direct method. It turns out that for any p > 1 and 
for all £ > 0 small enough the quadratic form 


H = £ — evv' + peuu + 


{p + 1)A£ 


(it'll 


uv') 


(7.20) 


is a strict Liapunov function for our system. The calculations are not immediate, espe¬ 
cially if we do not know in advance the formula! Here, LaSalle’s invariance principle 
was very useful since,without the information of asymptotic stability obtained by a 
very simple sequence of calculations, it would have been very difficult to imagine that 
such a function can be devised. 
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7.6 Two infinite dimensional examples 


Example 7.6.1. Let us consider the coupled system of second order scalar ODE: 


u" + u' + Au + cu = 0, 
v" + Av + cu = Q, 


(7.21) 


where A is a possibly unbounded linear operator on a Hilbert space H with norm 
denoted by |. | such that for some A > 0, 

A = A*>XI 


and c ^ 0 with < A^. In addition we assume that the unit ball of is 

compact in H. Let 


E{t) = £{u, u', V, v') 




+ c{u,v). 


We have the formal energy identity : = —\u'\'^. Since E is non-increasing, 

the vector {u,v,u',v') is bounded in D{A^/^) x D{A^/‘^) x H x H. Actually it is 
not difficult to check that if (u(0), u(0), m'( 0), t;'(0)) G D{A) x D{A) x D{A^/'^) x 
= W, then the vector {u,v,u',v') remains and is bounded in W for f > 0 
and the energy identity is rigorously satisfied. Then the trajectory is precompact and 
if (u, v) be a solution of (17.211) for which E is constant. Then u' — 0, hence u is 
constant and u" — 0. Then by the first equation r; = — ^ is also constant. Finally 
since by the hypothesis < A^ , the stationary system Au + cv = cu + Av = 0 
has no non-trivial solution, we conclude that u = v = 0 and therefore the solution 
tends to (0, 0,0, 0) . In fact, the system generates a uniformly bounded semi-group 
in X D{A^I'^') x H x H and it is then easy to conclude that (0,0, 0,0) is 

asymptotically stable. For the exact nature of the convergence we refer to 0 


Example 7.6.2. Consider the nonlinear wave equation 


Utt — Uxx + giut) =0 in X H, m = 0 on R.'*" x dVl (7.22) 


where H = (0,L) is an bounded interval of R and p in a non-decreasing locally 
Lipschitz continous function on R which satisfies gifi) — 0 and does not vanish 
identically in any neighborhood of 0. Then for any solution u of (17.221 ) such that 
(w(0),ut(0)) G D Hq(0,L) X Hq{0^L) = W, the vector {u,ut) remains and is 
bounded in W for f > 0 and we have 

^ [ {u^{t,x) + ul{t,x))dx =-2 [ g{ut)ut{t,x)dx 
™ Jn Jn 

By using the fact that for every regular solution v of the usual string equation, Vt(t, x) 
is 2L-periodic with mean-value 0, the invariance principle now shows that (rt, ut) tends 
to (0,0) in Hq{0, L) X L^(0, L) as t tends to infinity . 

Remark 7.6.3. The analog of the last example is valid in higher dimension in a much 
more general context, relying on the theory of monotonicity in Hilbert spaces and the 
concept of almost periodic functions. Since these methods fall outside the scope of this 
text, we refer to EH mi for the statements and proofs of the general results. 








Chapter 8 

Some basic examples 


In this chapter, we consider a few special cases in which asymptotic behavior can 
be studied completely by simple direct methods. These examples will serve later as 
models to undersand more complicated systems. 

8.1 Scalar first order autonomous ODE 

In this section we consider the simplest possible differential equation 

u' + fiu)=0, t>0 (8.1) 

The asymptotic behavior of bounded trajectories is obvious as shown by the following 
result 

Theorem 8.1.1. Let f € M). Each global and bounded solution u{t) of 

(IQ on R'*' tends to a limit c with f{c) = 0. 

Proof. If for some t > 0 we have /(itjr)) = 0, then u{t) = u{t) for all t and the 
result is trivial. If f{u{t)) never vanishes on K+, it keeps a constant sign and u{t) is 
monotone on R+. Since by hypothesis u{t) is bounded on R+, it follows immediately 
that u{t) tends to a limit c as f —>■ +oo. The equation shows that u'{t) tends to —/(c), 
and we conclude that /(c) = 0. □ 

8.2 Scalar second order autonomous ODE 

We now consider the slightly more complicated case of the equation 

u" + g{u') + f{u) = 0, t>0 (8.2) 

where f,g:M. —> R are locally Lipschitz continuous such that 

Vc s R \ {0}, g{v)v > 0. 

The term —giu') can be viewed as a dissipation while —f{u) represents a restoring 
force. We will show that convergence or divergence of the general solution of equation 
(18.2b depends on the strength of the dissipative term |p(u)| for small values of the 
velocity v. As a consequence of Corollary 7.5.1, (18.2b generates a gradient-like system. 
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8.2.1 A convergence result 

Theorem 8.2.1. Assume that f, g are as above and in addition, for some e G (0, 1] and 
S > 0, we have 

V?; G K, > 5inf{l, (8.3) 

Then ifu G is a solution of (18.21 . we have 

lim {|M'(t)| + \u{t) - c| = 0, 

t—>-+oo 


for some cG f ^{0}. 

Remark 8.2.2. A typical example of function g that satisfied hypothesis (18.31) is g{s) = 
|s|“s with a G [0,1). 

Proof. First, since the system is gradient-like, we have 

u}{uo,ui) C /“^{O} X {0}. 

By connectedness we have either oj(uo,ui) = {a} x {0} for some a G /~^{0} and 
the result is established, or 


w(uo,Mi) = [a, 6] X {0}, (a < b). 

In this case, we set c := As a consequence of the definition of uj(uo, ui), there 
exists a sequence (f„) of positive numbers such that 

lim tn = -foo, u(tn) = C. 
n—f+oo 

For any n G N, there exists 0 such that 

u(t) G [a, 6], Vt G (8.4) 

We claim that for all n G N large enough, we can take Sn = -l-c» in (18.41) . Indeed, let 

0n = Inf{t > t„,, u{t) ^ [a, b]} 

and assume < +c». Then we have 

Vf G [tn, On], u'ff) + g{u'{f)) = 0. (8.5) 

We may assume that n is large enough to imply |M'(f)| < 1 on oo[, so that from 
(18.51 we deduce as a consequence of (18.31 

'itG[tn,0n], |w'(f)| < {(1 - £)l5(f - -f . (8.6) 

In fact, if there is s G [f„, On] such that u'{s) = 0, then u'{t) = 0 for all t G [tn. On] 
and (18.61 is obviously satished. Otherwise u' has a constant sign, then 

= (^ _ l)u"|u'|^-2 sign (m') 

at 

= {l-e)g{u')[u'r^ 

> (1 — £)5. 


(8.7) 
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By integrating (18.71) over {t € [tn, Sn]), we get (18.61 ). Now from (18.61 ) we deduce 
by integration 

|m'(s)| ds < yt€[t„,9n]. 

For n large enough, the right-hand side is less than = |6 — c| = |a — c|. 

Therefore there exists ng S N such that we have 

Vn > no, u{t) G J, yt G [tn, +oo[ and 

\u{tn)-u{t)\ < ^\u'{tn)\^, ytG [f„,+00[. 

£0 

Since u(f„) = c for all n G N, we deduce 

Vn > no, Vt G [f„,+oo[, \u{t) - c\ < (8.8) 

£0 

Since u' [tn) —^ 0 as n —>■ -foo, it is clear that (18.81) implies 

lim \u{t) — c| = 0. 

t —¥-\-00 

Therefore J = {c} and this contradicts the hypothesis J = [a, h] with a < b. The 
proof of theorem r8.2.11 is completed. □ 

8.2.2 A non convergence result 

Theorem 8.2.3. Assume that there exists a,b G M. with a < b and a positive constant 
C such that 


/(s) <0, Vs < a 
/(s) = 0, Vs G [a, b] 

/(s) >0, Vs > & 

\g{v)\ < Cv^, V|t;| < 1. (8.9) 

Then for every bounded non constant solution of (IQ) . there exist a sequence tn —> 
+c» such that u(f„) < a for all n and a sequence 9n —> +oo such that u{9n) > b 
for all n. 

Remark 8.2.4. A typical example of function g that satisfied hypothesis (18.9b is g{s) = 
|s|s. 

In the proof, we have to use the following lemma. 

Lemma 8.2.5. Let v G C'^ (K.’*', M) satisfying 

v'{0) > 0, v"{t) > -Cv'itf, \/t G M+, 

where C > 0 is a constant. Then v is nondecreasing and lim v{f) = +c». 

t—>-+oo 
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Proof. It is clear that v'{t) > 0 for t small enough. Let 

T = sup{r > 0, v'[t) > 0, Vf G [0,t)}. 
For all t G [0, T), we have 


dVv'{t)^ - ■ 

By integrating over (0, t), we get 


Vf G [o,r), v'{t) > 




( 8 . 10 ) 


If T < + 00 , we obtain that v'{t) > 0 in a right neighborhood of T which contradicts 
the definition of T. Then T = +oo and (IS.lOl l becomes 


Vf G [0, +oo), v'{t) > 1 . 


By integrating this inequality, we get the last part of the lemma. 


□ 


Proof of theorem l8.2.3[ Since the system in {u, v) is gradient-like we have 

lim |m'| -I- dist(M(f), [a, b]) = 0. 

t^ + OO 

Assume that u{t) > a for t > Iq. We must prove that u is constant. We distinguish 
two cases : 

- If u’{t) > 0 on [fo, +oo), then u is nondecreasing and tend to c G /~^({0}). So 
f{u{t)) = 0 on [to, +oo) and we have 

u” + g{u') = 0, on [fo, +oo). 

If u' = 0 on [to, +oo), then u is constant. Otherwise, there exists ti > to such that 
u'iti) > 0. Applying lemma l8!2.5l to v{t) := u{t + ti), we get a contradiction. 

- If there exists ti > to such that u'{ti) < 0, therefore (since u{t) > a when 

t > to) 

u" + g{u') = —f{u) < 0 on [fo, +oo). 

In particular, u" < —g{u') < Cu^ on [fo, Tcxc) and then v{t) := —u{t + ti) verify 
v'{0) > 0, v"{t) > -Cv'itf, Vf G M+. 

Applying lemma [8'.2.5l to v, we get a new contradiction. □ 


8.3 Contractive and unconditionally stable systems 

In this section, {Z, d) denotes a complete metric space and we consider a dynamical 
system {S'(f)}t>o on {Z, d). The main result is as follows. 
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Theorem 8.3.1. Assume that the system {S(t)}t>o B unconditionally stable in the 
following sense 

Ve > 0,3^ > 0, y) G X X X, d{x, y) < 6 => sup d{S(t)x, S(t)y) < e. 

t>o 

( 8 . 11 ) 

Let T be given by (lO) . Then ifuo G X generates a precompact trajectory under S{t) 
and if uj(uo) fl 7 ^ 0, the trajectory S(t)uo converges to some limit a G X as t ^ 00 . 

Proof. Let a G w(uo) D X. Given any e > 0 and choosing (5 > 0 so that (18.1 lb is 
fulfilled, by definition there is r > 0 for which 

d{S{T)uo,a) < S 


Then we have 

Vf > T, d{S{t)uo, a) = d{S{t — t)S{t)uo, S{t — r)a) < e. 


□ 

Remark 8.3.2. Actually the above proof shows the following more general result; if 
uq G X generates a precompact trajectory under S(t) and if u!{uo) n X contains a 
stable equilibrium point a, the trajectory S(t)uo converges to a G X as t ^ 00 . 

A classical class of unconditionally stable systems is the class of contractive sys¬ 
tems; 

Definition 8.3.3. A dynamical system {S{t)}t>o on (Z, d) is said to be contractive if 

y{x,y) G X X X,yt > 0,d{S{t)x,S{t)y) < d{x,y) (8.12) 

An obvious consequence of Theorem l8.3.1l is the following 

Corollary 8.3.4. Assume that the system {S{t)}t>o is contractive. Then if uq G X 
generates a precompact trajectory under S{t) and if oj{uo) fl 7 ^ 0, the trajectory 
S{t)uo converges to some limit a G X as t ^ 00 . 

More generally, we have 

Corollary 8.3.5. Assume that the system {S{t)}t>o is such that for some M > 1 

y{x, y) G X X X,yt > 0 , d{S(t)x, S{t)y) < Md{x, y) (8.13) 

Then if uq G X generates a precompact trajectory under S(t) and if ui{uf) fl 7 ^ 0, 

the trajectory S{t)uo converges to some limit a G X as t ^ 00. 

Theorem l8.3.1l especiallv applies to gradient-like systems. 

Definition 8.3.6. A dynamical system {S{t)}t>o on (Z, d) is said to be gradient-like if 
whenever uq G X generates a precompact trajectory under S(f), we have uj(uo) C X. 

Corollary 8.3.7. Assume that the system {S{t)}t>o is gradient-like and uncondition¬ 
ally stable. Then ifu^GX generates a precompact trajectory under S{t), the trajec¬ 
tory S{t)uo converges to some limit a G X as t ^ 00 . 
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Remark 8.3.8. If we consider the ODE 

u" + u = 0 

written on as a system 

u' = V, v' = —u 

it is easy to check that any trajectory starting from Uq = {uq^vq) ^ ( 0 , 0 ) is non- 
convergent. Here S {t) is an isometry group on , hence trivially contracting. What 
happens here is that the system is not gradient-like. More precisely, whenever Uq = 
(mq, Uq) ^ (0, 0), we have oj{Uo) fl = 0 since = {0} and the norm of S{t)Uo is 
constant. 

As a basic application of theorem [8. 3. 41 Let N > 1 and F G C'^(K^) be convex. 
We consider the equation (I6.71 l 


u'{t) + XF{u{t)) = 0 


We obtain 

Corollary 8.3.9. Assume that £ = {z G ]R^,VF(z) = 0} ^ 0. Then any solution 
u(t) of (16.71) is bounded on R'*' and converses, as t ^ oo to some limit a G £ = Iz G 
R^,VF(z) = 0}. 

Proof. We already showed that the dynamical system S{t) generated by (16.7b on the 
closure of the range of u is gradient-like with set of equilibria £. Under the hypothesis 
that F is convex, it is easy to check that the operator VF G (R^, R^) is monotone, 
which means 


V(u, u) G R^ X R^, {VF{u) - VF{v),u -v)>0 
Then if {u, v) are 2 solutions of (16.7b . we have 

Vf > 0, ^l|w(i) - v{t)f = -2{VF{u{t)) - XF{v{t)),u{t) - v{t)) < 0 

Hence the system S{t) is contractive in the usual norm. In particular, since any a G 
£ = {z G R^, VF(z) = 0} is a solution of (16.7b independent of t, the function 

t !->• ||M(f) — a|| 

is non-increasing and all trajectories are bounded. Finally Corollary 18.3.71 gives the 
result. □ 

Remark 8.3.10. A much more general convergence result holds true for the equation 

0 G u' -f 9 $(m) 

where i9$(u) is the (possibly multivalued) subdifferential of any proper convex Isc 
function with arbitrary domain on a Hilbert space H, cf. Brack ED. In general only 
weak convergence is obtained, cf 0 . Besides, the asymptotic behavior of precompact 
trajectories of nonlinear contraction semi-groups has been the object of intensive study 
in the seventies, cf. e.g. ESlISSlIlll. 






8.4. THE FINITE DIMENSIONAL CASE OF A RES ULT DUE TO ALVAREZ 79 


8.4 The finite dimensional case of a result due to Al¬ 
varez 


In this section, we consider the equation (I 6 . 8 I 1 

u”{t) + u{t) + 'VF{u{t)) = 0 

where > 1 and F G C'^(K^) is convex. In contrast with the gradient system (I6.71 i. 
the system generated by (I 6 . 8 I 1 is gradient-like but generally non-contractive. However 
we have a convergence result similar to Corollary 18.3.91 which is a special case of a 
more general weak convergence theorem due to Alvarez, cf. Q. 

Corollary 8.4.1. Assume that £ = {z G ,yF{z) = 0} ^ 0. Then any solution 
uit) of (16.8b is elobal, bounded on R+ and converses, as t ^ 00 to some limit a G 
£ = {zG R^,S7F{z) = 0 }. 

Proof. From our Hypothesis it follows that F is bounded from below. First we con¬ 
sider a local solution u of (| 6 JJ on some interval [0, L) Given any positive T < L, the 
identity 

£ \\u'{t)rdt+= Fium - F{u{t))+ 

shows that u' G L°“(0,T;R^), therefore the solution is global and uniformly Lip- 
schitz. In addition u' G X) with X = We already showed that if all 

solutions U = (u,u') are bounded, the system S(t) generated by (16.8b is gradient¬ 
like and the set of fixed points of S(t) is F = £ x {0}. We now show that in fact 
u is bounded and the numerical function (p(t) = ||u(f) — a|p has a limit at infinity 
whenever a G £. Indeed a straightforward calculation shows that ip G with 

p” Fp' = -2{XF{u(t)), u{t) - a)) + 2\\u'{t)f < 2||M'f = hG L\R+) 


Writing this inequality as 

{e^ip'Y < e*h(t) 

provides 

p'{t) < e~*p'{0) + f e^~*h{s)ds := H{t) + e~*(p'(0) = K{t) 
Jo 


Now we have 


pi pi pt pi pi 

/ H{t)dt = / e^~*h{s)dsdt = / e'^h{s) / e~*dtds 

Jo Jo Jo Jo Js 

= [ — e~^)ds < ( h{s)ds 

Jo Jo 


Thus H,K G and since p > Q, the function xpit) := p{t) — K(s)ds 

is bounded with non-positive derivative. It tends to a limit at infinity and so does 
p. In particular u is bounded, and since S{t) is gradient-like, the omega-limit set is 
contained in F. Picking (a, 0) G uj{Uo) , the limit of p at infinity is 0 and we end up 
with convergence of u to a and u' to 0 


□ 




Chapter 9 

The convergence problem in 
finite dimensions 


9.1 A first order system 

In this section we consider the first order gradient system 

u' + Vipiu)=0 (9.1) 

where (p : —>■ K is assumed to be C^, and we set 

5 = {a e = 0}. 

As we saw in Section l631 any bounded solution of (19.11) approaches the set 5 as t goes 
to infinity. The question is then to determine whether or not it actually converges to a 
point in S. The next result shows that this is not always true. 


9.1.1 A non convergence result 

Theorem 9.1.1. Let k be a positive integer and let us consider 


^{x,y) = f{r,9) = 


_ 1 
e (TTsp' 

0 



4k^r* 


4fe2r4_|_(l_r2)2fc + 2 


sin(0 


1 

(l-r2)'= 


) 


ifr < 1 

ifr > 1 
(9.2) 


where we use the polar coordinates {x,y) = (r cos0, r sin0). Then there exists a 
bounded solution u of (igrT i whose ui-limit set is homeomorphic to S^. 


Proof. For = 2, by setting u = {x, y) equation (19. li becomes 


The system (19.3b becomes 


x' + |f(a:,j/)=0, 

y' + = 0, 


r' + %{r,e) = 0 , 

= 0 , 


(9.3) 


(9.4) 
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We define 

Let tq G (0,1) and let r be the local solution of 


j.' _ 2fcr(l-r^)'^+|‘ ^ n 

^ 4fc2r"‘+(l-r2)2'=+2® 


r(0) = ro 
Clearly, r is global and satisfies 


Vf S (0,+oo), 0 < r{t) < 1, and lim r{t) = 1. 

t—>-oo 


Now if we impose that 


(9.5) 


(1 - 

then a straightforward calculation shows that (r, 6) is a solution of (19.41 1. Hence, the 
solution (r, 6) verifies 


lim r(t) = 1, lim 0(t) = oo. 

t—^oo t—^oo 

Clearly, the w-limit set of the trajectory u = (r cos 0, r sin 9) of ( I9.31 l with ip given by 
(19.2b and (r, 0) satisfying (19.5b is the entire circle {(r, 0) / r = 1}. □ 

Remark 9.1.2. We recall that a function / S , R) is in the uniform Gevrey 

class Gi+ 5 (R^,R) if there exists a constant M = M{f) >0 for which 

Vm G ||i9™/||z,oo < 


where 


N 

\m\ := 

i=i 


is the length of the differentiation index m. It is natural to conjecture that, written 
in cartesian coordinates, ip G Gi+i outside any ball centered at 0 and therefore pp G 
^ 1 +^ (R^, R) for any p G ^ 1 +.^ (R^, R) which vanishes in a small ball around 0 and is 
equal to 1 ouside the ball of radius e < 1. If the conjecture is valid, this reinforces to the 
stronger regularity class Gi+ 5 (R^,R) C G°“(R^,R) with (5 = ^ the non-convergence 
result from J. Palis and W. De Melo which stated the existence oip G C°° (R^, R) 
for which there is a bounded solution u of ( 19.1b whose w-limit set is homeomorphic to 
S^. As (5 tends to 0 the space Gi+ 5 (R^, R) approaches the space of analytic functions 
Gi (R^, R) , showing that the next result is optimal if we look for a regularity class in 
which convergence of bounded trajectories is always true. 


9.1.2 The analytic case 

In 17011211, S. Lojasiewicz proved the following result which implies that the ’’bad” 
situation of Theorem l9 .1.11 cannot happen for analytic functions. 

Theorem 9.1.3. (Lojasiewicz Theorem iTQl [77]/ ) Let p : R^ —> R be an analytic 
function. Then for all a G S, there exists Ca > 0, aa > 0 and 0 < 9a < ^ such that : 

\\yp{u)\\ > Ca\piu) - p{a)\^~^‘' Vu G R^ ||u - all < (Ta. (9.6) 
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Remark 9.1.4. In the sequel, 9a will be called a Lqjasiewicz exponent of (p at point a. 
Each 9' < 9a is also a Lojasiewicz exponent of ip at point a, associated to a possibly 
smaller radius a < a a- Moreover when considering 9' and reducing a if needed, the 
constant Ca can be replaced by arbitrarily large constants, in particular by 1. This was 
the choice made by Lojasiewicz in his pioneering paper. On the other hand, in the cases 
where an optimal (= largest) 9 can be reached for instance by a direct calculation, it 
may happen that the choice c = 1 is irrelevant. For instance if TV = 1 and if{u) = ev?, 
we have || V(/3(u)|| = 2e|u| so that in particular 

\\Vp{u)\\=2e-^p{uf-^^ 


In this case the optimal value 0 = ^ is associated to a maximal constant cq which tends 
to 0 with the parameter e. Similar examples can be built with any super-quadratic 
power function. 

Remark 9.1.5. If a ^ 5, the inequality becomes trivial since p is of class C^. 

Theorem 9.1.6. (Lojasiewicz Theorem ^TQl [771/ ) Assume that p satisfies (i931 l at any 
equilibrium point a and let u G be a solution of dO- Then there exists 

a G S such that 

lim ||m(/) - all = 0. 

i^+oo 

Moreover, let 9 be any Lojasiewicz exponent of p at point a. Then we have 


0{e for some 5 > 0 if 9 = 

0(/-®/(i-2»)) if0<9 < i. 


In particular if p is analytic , all bounded solutions of (i9Tl l are convergent. 
Proof. We define the function z by z{t) = p{u{f)). Then 


(9.7) 


z'{t) = -\\Vp{u{f))f, V/> 0 . (9.8) 


So z is nonincreasing. Since u is bounded and p is continuous, it follows that K = 
lim p{u{f)) exists. Replacing p hy p — K may assume iV = 0. If z(io) = 0 for 

t—foo 

some to > 0, then z{t) = 0 for every t > to, and therefore, u is constant for t > to. 
In this case, there remains nothing to prove. Then we can assume that z{f) > 0 for all 
t > 0. 

Define T := uj{u). Theorem I4.1.8l iil implies that T is compact and connected. Let 
a G T, then there exists —>■ -foo such that u(tn) —> a. Then we get 


lim p{u(tn)) = p{a) = K = 0. 

n—>-+oo 

On the other hand, p satisfies the Lojasiewicz inequality ( 19.61 ) at every point a G S. 
Applying Lemma fl .2.6! with W = X = M.^, E = p and Q = Vpwe obtain, 

3cr,c> 0, 39 G (0, i]/ [dist(u,r) < cr =)> ||V(/5('u)|| > c\p{u)\^-^] . 

Now since T = by Theorem l4.1.8l iiil. there exists T > 0 such that dist(M, T) < 
a. Then we get for all f > T 


||Vv?(u)|| > c|(/5(u)|i ^ 


(9.9) 
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By combining (I9.81 l and (I9.91 l. we get 

z\t) < Vi > T. (9.10) 

In the case 9 G (0, ^), by integrating (19.101 1 over (T, t) we find 

z{t) < --j — < , Vi > T. 

(z(T)2e-i + (1 - 20)c2(i - r))T^ 

Now since 

||u'(i)f = -z'(i) 

we have 

||M'(s)|p(is = z(i) — z{2t) < Cit~ . 

Applying Lemma [l.2.5l to v(t) := ||u'(i)||, we get 

\\u'{s)\\ds < C 2 t~. (9.11) 

By Cauchy’s criterion, a := lim u{t) exists and 

Vi > T, ||M(i) — a|| < C2t~ . 

On the other hand, if 0 = i. the application of Lemma 11.2.41 to v(t) := 
gives the exponential decay. To conclude the proof, we remark that at the end, the 
global Lojasiewicz exponent used to prove convergence can be replaced by any local 
Lojasiewicz exponent of tp at a. □ 

Remark 9.1.7. Since the Lojasiewicz theorem is actually local, it suffices to assume 
that tp is analytic in a ball where the solution stays for all i. 

9.2 A second order system 

We now consider the gradient-like system 

m" + u'+ VT>(m) = 0 (9.12) 

where <1> : ^ R is assumed to be C^, and we set 

S = {aGR^, V$(a) =0}. 

9.2.1 A non convergence result 

The non-convergence result of Curry - Palis - De Melo (cf. Theorem 19.Til l has been 
extended to (19.12b by Veron ll79l (see also ||8]|65l)- More precisely 

Proposition 9.2.1. Given any (p G 1 < k < oo, there is a $ G 

such that each solution of (19.1b is at the same time a solution of (19.12b . 

Proof. The statement is readily satisfied for ^ = (p — |Vi^|2/2. □ 

Corollary 9.2.2. There exist $ G and a bounded solution u of (19.12b 

whose Lo-limit set is homeomorphic to S^. 

Proof. Take (p as in Theorem 19. 1.1 1 Then (19.1b has a bounded solution u whose w- 
limit set is homeomorphic to S^. By Proposition 19.2.11 u is also a solution of (19.12b 
for some smooth $, which proves the corollary. □ 
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9 . 2.2 A convergence result 

Theorem 9.2.3. Assume that 4) is analytic and let u G be a solution 

of (19.121) . Then there exists a G S such that 

lim |lu'(t)|| + ||M(i) - all = 0. 

t—¥-\-00 

Moreover, let 9 be any Lojasiewicz exponent of p at a. Then we have for some constant 
C>0 

\\u{t) - a\\ < , if0<9<^ 

\\u{t) — a|| < C exp(—i5t), for some i5 > 0 if 9 = 

Proof. Let E(t) = 5 l|w'(i)|P + $(u(t)). We have 

= {ufu') + {V^u),u') 

= {u" + V<I>{u),u') = -\\u'it)\\^ 

From Theorem l4.1.8l iil we know that uj{u, u') is a non-empty compact, connected set. 
We also know that limi_>+oo |1 m^|| = 0 and u}(u, u') C S x {0} (see corollary 16. 4. Ik 
Let r = {a/ (a, 0) G u!(u, m')} and K = limt_>oo E{t). As in the proof of theorem 
19.1.61 we may assume K = 0 and for all a € L, $(a) = 0. Then we introduce 

H{t) = l||u'(f)f + 4>(u(f)) +e(V4>(u(f)),u'(f)) 

where £ is to be fixed later. Therefore 

Hft) = -||u'f+e(V$(u),u")+e(V2$(w)w',w') 

= -||u'f + e(V$(u), -u' - V$(u)) + e(V2$(u) • u', u') 

= -||u'f - e||V$(u)f - e(V$(u),u') +e(v2$(u) • ufu'). 

Since u is bounded we have 

e{N^<S>{u)-u\u') < Ciellu'f. 

Thanks to Cauchy-Schwarz and Young inequalities we have 

Therefore selecting e < eo we find 

Hft) < -(l-C'2e)h'f-|||V$(u)f 

< -|(||u'f+ ||V$(u)f). (9.13) 

Then H is nonincreasing with limit 0, we have in particular H is nonnegative. As in 
the proof of the Theorem 19.1.61 we can assume that H(t) > 0 for all f > 0. On the 
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other hand, since $ is analytic then by using Lemma fl.2.6l once again as in the proof 
of Theorem l9.1.6l there exist 9 G (0, i], T > 0 such that for all f > T we get 

> C3(||«'f+ ||V$(u)f+ 

> (9.14) 
Combining the inequalities (19.131 1 and (I9.141 i we find 

H'{t) < 

\f 9 G (0, ^), intergrating this differential inequality we get 

H{t) < Cet~T^. 

When 9 = 5 , we find that H decays exponentially. 

Now from (19.131 1. we get 

p2t 2 

/ (llu'f+ ||V$(«)f)ds<-i?(f). 

Jt ^ 

The proof concludes exactly as in Theorem |9.1.6l □ 


9.3 Generalization 


The goal of this section is to give a general framework which covers the results of 
section |9.1.21 and [9. 2. 21 as well as some new examples. For this end, we consider the 
differential equation 

u{t)ET{u{t))=0, f>0, (9.15) 

v^hereT gC{R^;R^) . 

Theorem 9.3.1. Let u G (Ci(R+;R^) be a bounded solution of the differential equa¬ 
tion (19.15b . Assume that there exists a function E G C'^(R^), (3 > 1, 9 G {0, 1) and 


c,ci,T>0 such that 

P{1-9)<1, (9.16) 

£{u{t)) > 0 for every t > T, (9.17) 

{N£{u{t)),T{u{t))) > c\\N£{u{t))\f ||J'(M(f))|| for every t> T (9.18) 
II Vi?(m(<))|| > ci£{u{t)Y~^ for every t > T (9.19) 

for every a G one has : V£{a) = 0 ^ J^(a) = 0, (9.20) 

Then there exists a G such that lim u(t) = a. 

t—¥CO 

If moreover, £ satisfies for some C 2 > 0 

||•7^('t^(f))|| > C 2 £{u{t)y~^ for every t > T, (9.21) 


Then, as t ^ 00 , 




for some <5 > 0 if (3 = 


l-/3(l-6) 


ifp>j^. 


\\u{t) - all 


(9.22) 





















86 CHAPTER 9. THE CONVERGENCE PROBLEM IN FINITE DIMENSIONS 


Proof. We apply Lemma fl.2.3l with X = R.^ and H{t) = £{u{t )). Let m be a solution 
of (19.15b which is continuously differentiable, then, by the chain rule, 

= -{V£iuit)),u'it)) = {V£{u{t)),X{u{t))). 

By using (19.18b . (19.19b and equation (19.15b we get for all f > T 


> c\\v£{umfmum\ 

> ccf£{u{t))^^^~^'>\\u'{t)\\. (9.23) 


This is condition (11.4b with rj := 1 — /3(1 — 9) (thanks to (19.16b r] > 0.) 

It follows that the function t i—^ £{u{t)) is nonincreasing. Now if £’(u(fo)) = 0 for 
some fo > T, then £{u{t)) = 0 for every t > Iq, and therefore, by conditions ( 19.18b . 
(19.20b and the equation (19.15b the function u is constant for t > to. In this case, there 
remains nothing to prove. Hence we can assume £{u{t)) > 0 for all t > T. This is 
condition (11.3b . By applying Lemma fl.2.3l we deduce the convergence result. Now we 
will prove the decay estimate (19.22b . From (19.23b we deduce for alH > T 


By integrating this last inequality we get 


/ OO 

||u'(s)|| ds 


< 


cr]c{ 


.£{u{t))T 


By using hypothesis (19.21b and equation (19.15b . we get 

[£{u{m^ = £[u{t)f-’^ < -II J-(«(f))|| = l||u'( 

C2 C2 


(9.24) 


(9.25) 


(9.26) 


Combining (19.24b and (19.26b . we obtain 

J^[£{u{t)f] < -77ccfc2[f’']^. 

Solving this differential inequality (we have to distinguish two cases = 1 or > 
1), we obtain the estimate 


£{u{t)y 




Combining this estimate with ( 19.25b . the claim follows. 


□ 


In the next subsections we discuss several applications of our abstract results. 
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9.3.1 A gradient system in finite dimensions 

We start by applying our abstract results to the gradient system 

u'{t) + Nip{u{t)) = 0, 

where ip G The system is a special case of (19.151 1 if we take T = Vip. 

The function p is nonincreasing along u. Now if rt is a bounded solution of the above 
gradient system and since p is continuous, it follows that pao = lim p{u{t)) exists. 

If we define £ by 

£{v) = p{v) - Poo 

we see that hypothesis ( 19.171 ) is satisfied for all f > 0. If tp is real analytic , then it 
satisfies Lojasiewicz inequality (19.61 ). Therefore by applying lemma fr.2. 61 with W = 
X = E = p, Q = Xp and T = u!{u) we get 

3r > 0, 3c> 0, 30 e (0, i]/ ||V(/5(u(f))|| > c|(p(u(f)) - \/t>T. 

Now it easy to see that all hypotheses of Theorem 19.3. II are satisfied (here /3 = 1). 
Then there exists a G such that lim u(t) = a and the estimate 

t—^OO 

r 0(e-5*) if0=i 

\\u(t) — a|| = •( 

" ^ ^ \ (9(f-«/(i-2«)) if0< i. 

We thus recover the result of Section l9.1.2l 

9.3.2 A second order ordinary differential system 

Let $ G and consider the second order ordinary differential system 

u''{t) + u'{t) + VT>(M(f)) = 0. (9.27) 

This system is equivalent to the first order system (19.151 ) if we define T : —>■ 

by 

Now let u G 1L^’°°(R+,R.'^) be a solution of (19.27b . We define the energy of this 
system 

E(t) = ^llu'(m^ + ^(n(t)). 

We know that the function E is nonincreasing and Eoo = limt_>.oo E{t) exists. It is 
also well known that uj{u,u') is compact connected subset of $“^({0}) x {0} (see 
corollary 16. 4. lb . Let e > 0, and define £ : R by 

£{u,v) := + ^(u) - Eoo+ e{V^{u),v)RN, u, v G R^, 


f V$(m) \ f V'^^{u)v \ 

V j )' 


so that 


X£{u, v) 
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Fix R> 0, and let M := sup||^||<^_|_]^ || V^<i>(u)||. Choose e G (0,1) small enough so 
that [M + i)e < i. Then for every u,v€ satisfying ||u|| < i? we obtain 

{'S/£{u,v),R{u,v))k2n 

= llull^ -£(V^$(u)t;,w)KN +£(t;, V$(u))Riv +£||V$(ii)|p 

> (l-M£-|)|Hp + |||V$Hf 

> a'diuf+ ||V$(u)f). (9.28) 

Since u'(f))] = —{S7£{u,v),R{u{t),u'{t))) < 0, Then the function 

t I —> £ {u{t), u'it)) is nonincreasing. Thanks to the fact that u' —> 0 as f —^ c», it 
follows that lim £{u{t), u'{t)) = 0. Then £ satisfy hypothesis (19.171 1. Moreover, 

>- + oo 

||Vf(u,u)|| + ||J-(u,u)|| < Cdlull + ||V<i>(u)||). (9.29) 

By combining (19.281 1 and (19.291 1. we obtain that 

{V£{u,v),R{u,v))k^n > cd|V£:('u,u)||||J'('u,u)||. 

This is condition (19.181 1 with /3 = 1. On the other hand, if V£’(a, b) = 0 then by (19.281 1 
we have b = 0 and V$(a) = 0, then J^(a, b) = 0, hence (19.201 i. 

Now if we assume that $ is analytic, then £ is also analytic and satisfies Lojasiewicz 
inequality (19.6b . Therefore by applying lemma fr.2.61 with W = X = E = £, 
Q — V£ and T = uj{u, u') we obtain 

3T>0, 3c> 0,36>G (0,i]/ \\V£iu{t),u'{t))\\ > c£{u{t))^-'^. (9.30) 

Then hypothesis (19.19b is satisfied. On the other hand, by using (19.29b we get 

||J-(u,u)|| = Hull + ||u +V<i>(u)|| > idlull + ||Vcl>(u)||) > ^||Vf(u,u)||. 

Combining this last inequality with (19.30b we obtain that hypothesis (19.21b is satisfied. 
Therefore by Theorem 19. 3. T1 lim (u(<), = (a, 0) exists. We thus recover the 

t—¥oo 

result of Section l9.2.21 

In ll58l . also the case of nonlinear damping was considered. The damping, however, 
should not degenerate in the sense that near 0 the damping is in principle linear. The 
case of degenerate damping which is the object of the next section has been considered 
by L. Chergui in Il26l . 

9.3.3 A second order gradient like system with nonlinear dissipa¬ 
tion 

Let $ G C'^(R^, R) and consider the second order ordinary differential system 

u"{t) + g{u'{t)) + V^{u{t)) = 0, (9.31) 

where g G C(R^,R^) satisfying 

{ 9 (,v),v) > c||u||“+^ (9.32) 

||5(t^)ll < C'||u|r+i (9.33) 

and a > 0. 
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Theorem 9.3.2. We suppose that 

30e]O,i], Vae .5, 3aa >0/||V$(w)|| > |$(m) - $(a)|l-^ yuGB{a,aa). 

(9.34) 

Assume that a G [0, and let u € M^) a solution of ( I9.311 I. Then there 

exists a G S such that 


lim (||M(t)|| + ||u(t) - all) = 0. 

t—f + OO 


We also have 


\\u{t)-a\\=0{t i- 2 e+-(i”«) ) 


Proof. First of all, we define the energy of this system 


E{t) = ^\\u'{t)r+^u{t)). 


We know that the function E is nonincreasing and E^o = liint_>oo E{t) exists. It 
is also well known (see corollary 17.5. Il l that uj { u , u ') is compact connected subset of 
(V$)“^({0}) X {0}. In order to apply Theorem 19. 3. II we must write equation (19.311 1 
as a first order system (19.151) . This is the case if we define T : —>■ by 


T{u,v) := ( .s ^ , u, V G K^. 

V 5(t^) + V4>(M) J 

Let e > 0, and define £ : —>■ M by 

£{u,v) = i||t;f+ 4>(M)-L;oo+£||V$(u)||“(V4>(M),'u)Riv, u,vGR^ 

so that 

+ e|| V(^(a)||“V2$(u) • V + ea\\N 4){u)\\°^-^ {N (j){u), v)N^<I>(u) ■ Nfiu) 
\/t(u,v) - ,; + e||V(/)(u)||“V$(u) 

Let B C X K.^ be a suffiently large closed ball which is a neighbourhood of the 
range of {u, u'), then we have 


||.F(u,u)||<C'i(||t;|| + ||V$(u)||); (9.35) 

||Vf(u,t;)||<C2(||z;|| + ||Vcl>(u)||). 

Now choosing e G (0,1) small enough and by using Young inequality together with 
hypotheses (19.32b and (19.33b . we get 

{V£{u,v),E{u,v)) > C3(||ur+2 + ||V$(a)|r+2) > C4(|M| + ||V$(a)||)“+2. 

(9.36) 

Combining these three last inequalities we obtain 

{N£{u,v),Eiu,v)) > C 5 \\N£{u,v)r+^ \\Eiu,v)\\. (9.37) 


This is (19.18b with B = a+1. Since ^[£{u{f),u'{t))] = —{'V£{u,v),lF{u{t),u'{t))) < 
0, then the function t i—^ £{u{t), u'(t)) is nonincreasing. Thanks to the fact that 
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u' —^ 0 as t —7> oo, it follows that lim £(u(t), u(t)) = 0. Then £ satisfy hypoth- 

i^+oo 

esis (19.17b . Now if Vf (a, b) = 0, then by (19.36b 6 = V$(a) = 0 which imply by 
( 19.35b that T{a, b) = 0. This is hypothesis (19.20b . On the other hand by using Young 
inequality we get 

£{u,vy-^ < Cedlt^ll + ||V$(u)|| + !<!>(«) - i^ood-®). 

We also have 

||J-(u,z;)||>cr(||u|| + ||V<l>(u)||). 

Combining this two last inequalities together with the Lojasiewicz inequality ( 19.34b . 
we get 

||J'(u(f),u'(f))|| > c’£{u{t)y-^. 

This is (19.21b . Since a G [0, then /3(1 — 6*) = (a + 1)(1 — 0) < 1, then (19.16b is 
satished. Theorem l9.3.2l is proved. □ 


















Chapter 10 

The infinite dimensional case 


In 1781 . L. Simon completed the fundamental one dimensional result of Zelenyak ll83l 
and Matano lT^ by showing that the pioneering work of S. Lojasiewicz can be extended 
to some infinite dimensional context, among which the semi-linear parabolic equations 
with analytic generating function in any space dimension. The objective of this chapter 
is to clarify to which extent the Lojasiewicz method can be generalized to infinite di¬ 
mensional systems. Throughout this chapter, we consider two real Hilbert spaces V, H 
where V C H with continuous and dense imbedding and H', the topological dual of 
H is identified with H, therefore 

V CH = H' cV 
with continuous and dense imbeddings. 

Definition 10.0.3. ITe say that the function E G {V, K.) satisfies the Lojasiewicz 
gradient inequality near some point (p G V, if there exist constants 9 G (0, ^], c > 0 
and cr > 0 such that for all u G V with ||u — < cr 

\\DEiu)\\v' > c\Eiu) - E{p)\^-^. (10.1) 

Remark 10.0.4. 1) The Lojasiewicz gradient inequality is trivial if (p is not a critical 
point of E. 

2) The number 9 will be called a Lojasiewicz exponent (of E at p). 


10.1 Analytic functions and the Lojasiewicz gradient 
inequality 

One might wonder if Lojasiewicz gradient inequality is valid for any analytic function 
on an infinite dimensional Banach space. However, even if V = 77 it is not the case. 
Actually, if (77, (•, •)) is a Hilbert space and E is defined by F{u) = {Ku, u) with 
K = K* > 0 and compact, then E does not satisfy the Lojasiewicz gradient inequality. 
More precisely 

Proposition 10.1.1. Let 77 = P(N) and F : 77 —>■ R 7>e the continous quadratic 
(hence analytic ) functional given by 

OO 

F(^uo, Ml, ...M„,...) := ^ 

7=0 
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where (efe)fegfij is a real sequence satisfying > 0 and lim Sk = 0. Then F satisfies 

k—foo 

no Lojasiewicz gradient inequality. 

Proof. Defining {ei)j = 5ij, an immediate calculation shows that 
Vf > 0, F{tek)=fek-, \VF{tek)\ = 2tek. 

In particular for each 0 > 0 we have 

|Vf (fefc)l _ rs 0,20-1 

For any 0 > 0 small, choosing t small enough and letting k tend to inhnity we can see 
that the Lojasiewicz gradient inequality fails in the ball of radius t. □ 

More generally, in |[54l . we considered a real Hilbert space {H, (•, •)), a linear 
operator A such that 

A€L{Hy, A*=A 

and the associated quadratic form $ : H —> R defined by 

Wu G H, ^{u) = ^{Au,u). 

In this context a characterization of continous quadratic forms for which the Lojasiewicz 
gradient inequality is valid was obtained and expressed by the following statement 

Theorem 10.1.2. The following properties are equivalent 

i) 0 is not an accumulation point ofsp(A). 

ii) For some p > 0 we have 

Vu e ker(A)-^, \\Au\\h>p\\u\\h- 

Hi) $ satisfies the Lojasiewicz gradient inequality at the origin for some 9 > 0. 
iv) <i> satisfies the Lojasiewicz gradient inequality at any point for 9 = ^. 


For a general nonlinear potential F, one might wander if the equivalent properties 
above for A = D^F{a) are sufficient to obtain a Lojasiewicz gradient inequality near 
a. The proposition below shows that it is not the case. 

Proposition 10.1.3. Let H = P(N) and F : H ^ IS. be the analytic functional given 
by 

°° li;,|2fc+2 

Then F satisfies no Lojasiewicz gradient inequality. 


Proof. First we note that D^F{0) = 0, hence sp{D'^F{0)) = {0} and in particular 0 
is isolated in sp{D‘^F{Q)). Defining = 5ij, an immediate calculation shows that 


Vf > 0, F{tek) = 


^ 2fc+2 


(2fc + 2)!’ 

In particular for each 0 > 0 we have 


|VF(fefc)| = 


^2fc+l 


(2A: + 1)!' 


F{tei 


M-fl 


|VF(fefc)| 


= c(0,fc)fi-(2fc+2)e. 


Choosing k large enough gives a contradiction for t small. 


□ 
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In this example, the difficulty comes from the fact that dimker(i9^F(0)) = oo. 
Assuming the equivalent properties of Theorem ll0.1.2l and dimkerfU^f fOii oo 
is equivalent to the semi-Fredholm character of D^F{0) (cf. Theorem 1 1.3.31 1 In the 
next section 110.21 we shall show that this condition is sufficient in a rather general 
framework, in particular V will not be assumed equal to H in view of applications to 
semilinear PDF. 


10.2 An abstract Lojasiewicz gradient inequality 

The purpose of this section is to give sufficient conditions on E for the inequality 
(IIP. lb to be satisfied. Let E G and ip G V such that DE{(p) = 0. Up 

to the change of variable u = p + v and the change of function G{v) = E{p + 
v) — E{p), we can assume whithout loss of generality that p — 0, E{0) = 0 and 
DE{0) = 0. Although the formulation of the Lojasiewicz gradient inequality requires 
only U G K.), one way of proving it requires U G (^^(UjR). In fact the operator 

A := D‘^E{0) plays an important role. 


We start with the following very simple result 
Proposition 10.2.1. Assume that 

A G L{V, V') is an isomporphism. 

Then the Lojasiewicz gradient inequality is satisfied near 0 with the exponent 0 = ^ : 
there exist two positive constants a > 0 and c > 0 such that 

||m||v <a =P \\DE{u)\\v' > c|£’('u)|5. 

Proof. It is easy to see, using Taylor’s expansion formula, that for ||M||y small enough 
we have 

\E{u)\ < C\\u\\^y. (10.2) 

On the other hand, since DE{u) = Au + o{u), we have 

u = A~^ DE{u) + o{u), 

and therefore for any given e > 0 we can find ( 5 (e) > 0 such that if ||u||y < ( 5 (e) then 

\\DE{u)\\v,>\\A-T^\\u\\v-e\\u\\v. 

Choosinge := eo := ||A“^||“^/2, we obtain for ||M||y < (5(eo) 

\\DEiu)\\v'>eo\\u\\v. (10.3) 

The result follows by combining (110.2b and (110.3b . □ 

Remark 10.2.2. Since A = D‘^E{0) is symmetric, then if A is semi-Fredholm and 
d = dimker(A) = 0, by corollarv ll.3.6I A is an isomorphism. Hereinafter we assume 
that d > 0. We denote by H : U —> ker(A) the projection in the sense of H. 
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Proposition 10.2.3. Assume that A := D^E{0) is a semi-Fredholm operator and let 

AT :V V' 

u I—^ Ilu-\-DE{u). 

Then there exist a neighborhood of 0, W^i(O) in V, a neighborhood of Q, W2{0) in V' 
and a map 'P : VF 2 ( 0 ) —> Wi{Q) which satisfies 

A^('I'(/)) = / V/G 1^2(0), 

T!{Af{u))=u Vu G IFi(O), 

\\'^if)-'^i9)\\v<Ci\\f-g\\v' V/,gGW^ 2 ( 0 ), Ci > 0. (10.4) 

Proof. The function Af is and DAf{0) = If + D‘^E{0) which by coroll arv ll. 3. 6l is 
an isomorphism from V to 14'. We have just to apply the local inversion theorem. □ 

Let {(pi,(p 2 , ■■■Td) denote an orthonormal basis of ker(yl) relatively to the inner 

d 

product of H. For ^ G small enough to achieve G W2(0), we define the 


r(e) = F;(vi/(5]0^,))- (10.5) 

3 = 1 

Let 1F2(0) be the open neighborhood of 0 in such that 

d 

^g11^2(0)^^^,</p, G W2(0). 
f=i 

The function T is in ll2(0). Let us define also 

ILi(O) = {u€ lLi(0)/n(u) G 1 L2(0)}. 


d 


Proposition 10.2.4. Let u G IFi(O) and let ^ G ll2(0) such that n(M) = 
W 2 (0). Then there are two constants C, K > 0 such that 


^^3^3 

3=1 


G 


||Vr(0||M<i <C||i9i?(w)||y., (10.6) 

\Eiu)-TiO\ <K\\DEiu)fy,. (10.7) 

Proof. For any fc G {1, • • • d} we have the formula 


9r 

d^k 


d_ 

ds 




s=o — 


Now we claim that for all ^ G 112(0) 


{DEi^TiY, ^3T3)Tk). 

i=i i=i 

( 10 . 8 ) 


d 


HE 


ar 

d^k 


iOPk- DEiT>{YC3T3))\\v' <C2mm'i'iY^^^Mv'- 


3 = 1 


3 = 1 


(10.9) 
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d 

In fact by using (llO.Sl l. remarking that DE[^ipj)) G ker A we obtain 

i=i 

d op d 

k^l j^l 

d d d 

= \\Y<DEi'f{J2^^ j(pj)),D'Ii{YijTj){Tk) - Tk > TkWv- 

k=l j=l j=l 

Now by using Cauchy-Schwarz inequality, and the fact that D'l/ (0) (Cu) = u, the claim 
follows. On the other hand, since E is C^, there exists C 3 such that 

\\DE{u)-DE{v)\\v'<C3\\u-v\\v V(u, u) G kki(O). (10.10) 

Then by using (110.4b . (110.9b and (IIP.10b we obtain 


||vr(c)|lM. < c.WDEiTfiY^^^fiMv 

i=i 

= C^WDEiT'iUiumiv' 

= C 4 \\DE{'T{U{u))) - DE{u) + DE{u)\\v' 

< C4DE{u)\\v + C'3C4||4'(n(u)) - u\\v 

= C4\\DE{u)\\v' + C3C4||4'(n(u)) - T>{Uu + DE{u))\\v 

< C4\\DE{u)\\v'+C5\\DE{u)\\v' 


hence (110.6b . On the other hand 


\E{u)-m\ 


= |T;(u)-T;(vi-(n(u)))| 

= 1^ ^{Eiu + ti'Ii{U{u))-u)]dt\ 

= I [ {DE{u + t{T>{Il{u))-u)),'I/{Il{u))-u)dt\ 

Jo 

< ||'l'(n(u)) - u||y / \\DE{u + t{'Ii(n{u)) - u)\\v'dt 

Jo 

< [f (||i9ii;(u)||4^. + fC'3||4/(n(u))-u||y)df]||'k(n(u))-u||y 
Jo 

< C3\\DE{u)\\v'\\'fmu)) - T-(n(u) + DEiu))\\v 


hence (110.7b . 


□ 


Theorem 10.2.5. Assume that A := D^E{0) is a semi-Fredholm operator and let 
d = dim ker A. Assume moreover that 

(Hl)d > 0 and there exists O C R'^ open, and h G C^{0, V) such that 0 G h{0) C 
{DE)~^{0) and h : O —> h(0) is a dijfeomorphism. 

Then there exist two positive constants cr > 0 and c > 0 such that 


\\DE{u)\\v' > c|£'(u)| 2 . 


||u||y < cr 














96 


CHAPTER 10. THE INEINITE DIMENSIONAL CASE 


Proof. We have by using (I10.9t (choosing a smaller W2(0) if necessary) 

d 

< C'r||Vr(C)||. (10.11) 

i=i 


If u G W'r(O) such that DE(u) = 0, then M(u) = 
'I'(n(u)). Moreover by using (IIO. 6 I 1 we have Vr(^) 

d 

Hu = 


f=i 


!!(«) which implies that u = 
= 0 where ^ G fF 2 ( 0 ) with 


d 

On the other hand let ^ G W2(0) with Vr(^) = 0. Then ^jPj) £ W^r(O) and 

i=i 


d d 

DE{'T{'^^j^j)) = 0 by using (IIP.lit . So n(T'(^ 
i=i i=i 


d 

Conse- 

i=i 


d d 

quently rh 0^,) G W-i(O) and DEi^iY^ 0^,)) = 0. 

i=i i=i 

Finally we have: 


d 

{u G VFi(O), DE{u) = 0} = e e W2{0) and Vr(e) = 0}). (10.12) 

Now we introduce the d—dimensional manifold 


7 = h{0) 


with O and h as in (HI). Let 

O = h-^{{u G W^i(O), DE{u) = 0}). 

Clearly O is an open subset of K.'^ and 0 G h{0). 

We now have 

d 

7 := h{d) c{ue WdO), de{u) = 0 } c e e w^m)- 

f=l 

Since the extreme terms are d—dimensional open manifolds, they must coincide lo¬ 
cally. Therefore, changing if necessary Wi(0) and IT2(0)) to smaller open sets, we 
obtain 

d 

7 = {u G ILi(O), DE{u) = 0} = C e ^^ 2 ( 0 )}. (10.13) 

i=i 

Now by comparing (110.12b and (110.13b . we get 

r(0 = o, v^gil2(o). 


The proof of Theorem ll0.2.5l follows immediately by using this last equality in (110.7b . 

□ 
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In the next theorem, we will prove inequality like (110.11 under hypotheses of ana- 
lyticity of E and DE. We consider a Banach space Z such that ker A G Z and Z C H 
with continuous and dense imbedding. 

Proposition 10.2.6. Assume that A := D^E{0) is a semi-Fredholm operator. Let 
£ := n + A Then W := C~^{Z) is a Banach space with respect to ||tu||w = 
and C £ L{W, Z) is an isomporphism. 

Proof. Using corollary 11.3.61 we know that C : V — V is one to one and onto. 
Since W GV and by the definition of W we also have £ : W —> Z is one to one 
and onto. Obviously we have £ £ L{W, Z) because ||£m||z = ||w||tv for all u &W. 
Now we prove that W is a Banach space. Let (wn) be a Cauchy sequence in W, then 
{C{wn)) is a Cauchy sequence in the Banach space Z. Denote by z its limit. {C{wn)) 
is also a Cauchy sequence in V, so (ui„)) is also a Cauchy sequence in V. Denote by 
w its limit, since £ £ £(U, V), then Cw = z. The claim is proved. Banach’s theorem 
gives the fact that £“^ £ L{Z, W). □ 

Theorem 10.2.7. Assume that A := D^E{0) is a semi-Fredholm operator and that 
N := ker A C Z. Assume moreover that : 

(H2) E : U ^ is analytic in the sense of definition M.4. l\ where U G W is an open 
neighborhood o/O, that DE{U) G Z and DE : U —> Z is analytic. 

Then there exists 9 £ (0,1/2], cr > 0 and c > 0 such that 

\\u\\v<a^\\DE{u)\\v'>c\E(u)\^-^. 


Proof. For the proof we need the following result. 

Lemma 10.2.8. Then there exist a neighborhood o/O, Vi (0) in W, a neighborhood of 
0, 1 / 2 ( 0 ) in Z and an analytic map T*! : V2(0) —> 1^(0) which satisfies 


Af(vki(/)) = / V/ £ U2(0), 

T>i{Af{u))=u Vu£Vi(0), 

in U2(0)nW2(0) 

ll^(/) - ^(5)lk < C[\\f-g\\z V(/,5) G U 2 ( 0 ) nlU 2 ( 0 ), 


Proof. We first establish that 


Af:W 


Z 


u 


IIm + DE{u). 


(10.14) 


is a diffeomorphism near 0, because DAf{0) = 11 + A = £ £ L{W, Z) is an 
isomorphism (see r)roposition ll0.2.6t and the classical local inversion theorem applies. 
Therefore we can find a neighborhood Vi(0) of 0 in kU and a neighborhood ¥ 2 ( 0 ) of 
0 in Z such that Af : Vi(0) —V2(0) is a diffeomorphism. Finally it is clear that 
vki = in 1 / 2 ( 0 ) n 11 / 2 ( 0 ). By Theorem ll.4.9l we have 4/1 is analytic in 1 / 2 ( 0 ). □ 

End of proof of Theorem By using the chain rule (Theorem I 1 .4!6l) . since E : U — > R, 
DE : U — > Z and 4* : 1 / 2 ( 0 ) H 112(0) —> 1/(0) are analytic , the function F defined 
in (110.5b is real analytic in some neighborhood of 0 in 

Applying the classical Lojasiewicz inequality (Theorem 19.1.3b to the scalar analytic 
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function F defined on some neighborhhod of 0 in by tbe formula (I10.5t . we now 
obtain (since (1 — 6*) G (0,1)): 

< |r(c)r-^ + |r(c) - E{u)\^-^ < ^||vr(e)|U. + |r(0 - Eiu)\^-^. 

'^0 

(10.15) 

By combining (110.6b . (110.7b . (110.15b we obtain 

\E{u)\^-<^ < ^\\DE{u)\\v'+K^-^\\DE{u)\\l^}-^\ 

Then since 2(1 — 0) > 1, there exist cr > 0, c > 0 sucb that 

\\DE{u)\\v' > c\E{u)\^~^ for all m G y such that ||u||y < cr. 
Theorem ll0.2.7l is proved. □ 


10.3 Two abstract convergence results 

This section is exceptionnally devoted to an abstract situation in which a trajectory of 
some evolution equation is known independently of any well-posedness result for the 
corresponding initial value problem. In particular there is no underlying continuous 
semi-group to rely on and we cannot apply directly the simple results of chapters 4 and 
6 . However, by performing essentially the same kind of calculations as those needed 
to apply the invariance principle, we end up with a “gradient-like” property which is 
the starting point for the Lojasiewicz method to be applicable. Our results contain 
as special cases the semi-linear examples of section 10.4 (for which the semi-group 
framework could be applied as an alternative method) but they can also be used for 
strongly non-linear problems as soon as a solution with the right regularity properties 
is known, even if the well-posedness is either false or presently out of reach. 

Let V and H be two Hilbert spaces such that V is a dense subspace of H and the 
imbedding of V in 77 is compact. We identify H with its topological dual and we 
denote by V the dual of V, so that H cV with continuous imbedding. 

Let E G R). We study the following two abstract evolution equations; the 

first order equation 

u'{t) + VE{u{t)) = 0, f > 0 (10.16) 

and the second order equation 

u”{t) + u'{t) + VE{u{t)) = 0, f > 0 (10.17) 

Theorem 10.3.1. Let u G C^(]R+, V) be a solution of (110.16b . and assume that 

(i) Ut>i{u(f)} is compact in V; 

(ii) E satisfies the Lojasiewicz gradient inequality near every point ip G 5 := G 
y NE{p) = 0}. 

Then there exists p G S such that 

lim ||M(f) - (^||y = 0. 

>- + oo 
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Moreover, let 9 be any Lojasiewicz exponent of E at p. Then we have 


IkW - tWh = 


0{e forsome5>0if9=^, 

0(t-®/(i-2®)) if0<9< i. 


(10.18) 


Proof. We define the function z by z{f) := E{u{t)) for all f > 0. Since u G 
(7^(18.+ , V) and E G C^{V, K.), then by chain rule, z is differentiable and 

z'{t) = -\\u'[t)\\'^H, Vf>0. (10.19) 


Integrating this last equation and by using (i), we get u' G L^(]R+; H). Now, since 
the range of u is precompact in V, and u is uniformly Holder continuous on the half¬ 
line with values in H, it is also uniformly continuous with values in V and u' = 
—VE(u(t)) is uniformly continuous with values in V'. Then by applying Lemma fl.2.21 
to the numerical function Wu'itWv , we obtain that u'{t) tends to 0 in V as t tends 
to infinity, hence also in H by compactness. We conclude that oj(uo) C S. Moreover, 
since the function z is bounded and decreasing, the limit K := lim E{u{t)) exists. 

t—fco 

Replacing E hy E — K we may assume K = 0. 


If z(fo) = 0 for some to > 0, then z(t) = 0 for every t > to, and therefore, u is 
constant for t > Iq. In this case, there remains nothing to prove. Then we can assume 
that z{t) > 0 for all t > 0. Define T := uj{uo). It is clear that T is compact and 
connected. Let ip G T, then there exists —>■ +oo such that \\u{tn) — ip\\v —> 0. 
Then we get 

lim E{u{tn)) = Eiyp) = K = 0. 

n—>-+oo 

On the other hand, by assumption (ii), E satisfies the Lojasiewicz gradient inequality 
(IIP. lb at every point (p G S. Applying Lemma fl .2.61 with W = V, X = V', and 
Q — VE we obtain, 

3(7, c> 0, 36 G (0, i]/ [dist(M, L) < cr => \\\7E{u)\\v' > c\E{u)\^~^] . 

Now since L = uj{uq), by Theorem l4.L8l iiii. there exists T > 0 such that dist(M, L) < 
(7 for all t >T. Then we get 

Vf>T ||VL;(u)||y/> ( 10 . 20 ) 

By combining ( 110.19b and (110.20b . we get 

zft) < -c2(z(f))2(i-«), Vf > T. (10.21) 

The end of the proof is identical to that of Theorem 19. 1.61 we obtain the convergence 
of u{t) in H and the convergence in V follows by compactness □ 

Theorem 10.3.2. Let u G C'^(R+,F) H be a solution of (110.17b and 

assume that 

(i) Ut>i{u(f), M'(f)} is compact in V x H; 

(ii) if K \ V' denotes the duality map, then the operator K o E"(v) G C{V) 
extends to a bounded linear operator on H for every v G V, and K o E" : V —>■ C{H) 
maps bounded sets into bounded sets; 
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(Hi) E satisfies the Lojasiewicz gradient inequality near every point ip G S := 

{ip e y, VEiip) = 0}. 

Then there exists ip G S such that 

lim llu'llff+ ||u(t)-(/5||y = 0. 


Moreover, let 9 be any Lojasiewicz exponent of E at ip. Then we have 


IkW - tWh = 




,-St\ 


for some S > 0 if 9 = i 


0(t-®/(i-2®)) if0<9< i 


( 10 . 22 ) 


Proof Let 

m ■.= \\\u'mi+E{u{t)). 

By the assumptions on u and E, E is differentiable everywhere and for all f > 0 

Eft) = 


Hence E is decreasing, and by using (z) it is bounded. By integrating the last equality, 
we deduce that u' G H). Since H ^ V' we deduce that h{t) := ||M^(f)||y/ 

is integrable. Moreover by assumption (i) and the equation (110.171) . for almost every 
f > 0 we find 

\h'{t)\<2\\u'{t)\\v,\\u"{t)\\v,<C 

Hence the function/i is Lipschitz continuous and integrable which implies lim h{t) = 

t—foo 

0. Since u' is compact with values in H we deduce 

lim \\u'(t)\\H = 0. 


Let {ipjf) G uj{u,u'), and let C M+ be an unbounded increasing sequence 

such that lim {u{tn),u'(fn)) = {ip, '*/')■ Obviously we get i/j = 0. On the other hand, 

n—>-oo 

since ||u'||rr —> 0, it follows that 


lim sup + s) — tpllrr = 0. (10.23) 

sG[0,l] 

Actually the same is true with values in V. In fact, assuming the contrary, there is 
<5 > 0 such that 

Vn G N, sup ||u(f„ + s) — > i5. 

sG[0,l] 

Then we can find a sequence (s„) C [0,1] such that 

VuGN, \\u{tn + Sn) - <p\\v > 

By compactness of u in V, we can find ip G V and subsequences still denoted (f„) and 
(s„) such that 

\\u{tn + Sn) - IpWv - 0 

which imply that — ip\\v > f . Now from (110.231) we deduce that ip = ip, a. contra¬ 
diction. 
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Therefore, lim \7E(u(tn + s)) = uniformly in s £ [0,1]. 

n ^oo 

By equation (110.171 1. 


S7E{ip) = / S7E{(p) ds 


= lim [ VE{u{tn + s)) ds 

= lim [ {—u"{tn + s) — u'{tn + s)) ds 

= lim -u{tn + 1) + u'{tn) - u{tn + 1) + u{tn) 

n—^oo 

= 0 . 


Hence ip G S. Now since £ is bounded and decreasing, the limit K := lim £{t) = 

t—^OO 

lim E{u{t)) exists. Replacing E hy E — K we may assume K = 0. 

t—^OO 

Now let e be a positive real number, and as in ll57ll let us define for all f > 0 

Z{t) = ^Wu'Wj, + E{u) + e{yE{u), u')v' (10.24) 

where (•) ■}v' denotes the inner product in V. We note that Z makes sense as a conse¬ 
quence of hypothesis (i). We have for almost all t > 0: 

Z'{t) = -\]urH + e{-{VE{u), u')v, - WVEiuWy, + ((Vi7(«))', 

Then, using (ii), for almost alH > 0 we obtain for some P > 0 

Z'it) < (-1 + Pe)\\u'fH - e{VE{u),u')v' - e\\VE{u)\\l.. 


Since we have by Cauchy-Schwarz inequality 

we deduce; 

z'{t) < (-1 + PeWfH + IWu'fy, - \\\VE{u)ry,. 

By choosing e small enough, we see that there exists ci > 0 such that for almost all 
t > 0 

Z'{t) < -c^{\\u'fH + \\^E{u)fy,). (10.25) 

Since Z is nonincreasing with limit 0, we have in particular Z is nonnegative. As in 
the proof of the Theorem ll0.3.1l we can assume that Z{t) > 0 for all f > 0. 

Let r = {if / (yi, 0) £ uj{u,u')}. Theorem 14.1.81 iil implies that T is compact and 
connected. Now by assumption {Hi), E satisfies the Lojasiewicz gradient inequality 
(IIP.lb at every point p G S. Applying Lemma Q .2.61 with W = V, X = V, and 
Q — VE we obtain, 

3a, c> 0, 36 G (0, i]/ [dist(rt, L) < cr =» \\\7E{u)\\v' > c\E{u)\^~^] . 

Now by the definition of L and using Theorem l4.L8l iiil. we obtain that there exists 
T > 0 such that dist(it, L) < u for all t >T. Then we get 

\/t>T \\XE{u)\\v' >c\E{u)\^-^. 


(10.26) 
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Using this last inequality together with Cauchy-Schwarz and Young inequalities, we 
get for alH > T 

< C^iWu'Wj, + \NE{u)\\l, + \Eiu)\y^^-^\ 

< C:,{\\u'rH + \\NE{u)rv,} (10.27) 

Combining the inequalities (110.251 1 and (110.271) we find for all f > T 

Z'{t) < 

The conclusion follows easily □ 


10.4 Examples 

10.4.1 A semilinear heat equation 

As a first application we study the asymptotic behaviour of the semilinear heat equation 


{ ut — Au + f{x, u) = 0, {t, x) S R+ X U, 

u(f,-)|an = 0, t€R+, (10.28) 

u{0,x) = uo(x), X G Tl. 

In equation (110.281 ) we assume that U C (N > 1) is a bounded domain. We 
assume that / : U x R —> R is continuously differentiable and if A > 2, we assume 
in addition that 


3C > 0, a > 0 such that [N — 2)a < 2 
and s)| < (7(1 + |s|“) a.e. on U x R 

With this condition on /, the energy functionnal E given by 


(10.29) 


Vm G Tfp (U), i?(M) = i [ \\7u\'^dx+ f F{u)dx, 

2 Jn JQ 

where F{x, s) := f(x, r) dr, is well defined. By using Proposition 1.17.5 page 66 
of ll66l . we know that E is of class on i7o (U) and 

DE(u) = —Au + f{x,u), Vu G Hq{D), 

D^E{u)^ = 

It is well known that D^E{tf) is a semi-Fredholm operator for all G Tfg (U) 

(see examr)le fl.3.7b . Let d = dimker DE{(p). 

Proposition 10.4.1. Assume that hypothesis ( 110.291 ) is satisfied. Let p G Hq{TI) D 
L°“(U) be a critical point of E. Assume also that one of the following hypotheses is 
satisfied: 


d = 0 (10.30) 

d > 0 and there exists (7 C R'^ open, and h G C^{0, V)/ (10.31) 

ip G h{0) C {DE)~^{If) and h : O —> h{0) is a dijfeomorphism; 
f is analytic in s, uniformly with respect to x G El (10.32) 
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Then there exist 9 G (0, and a > 0 such that 

WuGH^in), ||M-v?||ffi(n) < cr ^ ||-Au+/(a;,M)||//-i( 0 ) > \E{u)-E{ip)\^~^. 

(10.33) 

Proof. Let A := D‘^E{(p) and assume that d = 0. Corollarv ll.3.61 gives that A — 
D^E{ip) is an isomorphism from Hq{EI) into H~^{El). To conclude we have just to 
apply proposition 1 10.2.11 Now assume (110.31b holds. To apply Theorem llO.2.51 we 
have just to remark that A is a semi-Fredholm operator. For the proof of (110.33b under 
hypothesis (110.32b . we distinguish two cases : 

Case 1 : N < 3. Let Z = L^{Q), by elliptic regularity ID we get that W :=(!! + 
A)~^{Z) C H'^{El) where If is the orthogonal projection in L^(r2) on A( A) :=kerA. 
The functional E : H^{Vl) fl H^iyi) —> M is clearly analytic since it is the sum of 
a continuous quadratic functional and a Nemytskii operator which is analytic on the 
Banach algebra C (see example[TA2l) By using Proposition ! 1.4. 51 we 

also obtain that DE : W —> Z is analytic. We can apply Theorem IIP. 2. 7! to obtain 
(110331) . 

Case 2 : N > 4. Let p > ^ ^nd Z = U’[El). By elliptic regularity IH, we know 
that W := (If -f A)~^(Z) C W'^’'P(n) which is a Banach algebra since p > ^. The 
end is the same as in the hrst case. □ 

Remark 10.4.2. 1) The result of proposition ! 10.4.1 [ remains true for the general energy 
dehned by : 


1 . ° r Ti'ii rtn r 

T X! / -^ ueH^{n), (10.34) 


where E{u) = f{s) ds, f satishes ( !10.29b and aij satishes the following conditions 


1. Uij e (7^(11), 

2. Qij = aji, and 

d 

3. ^ > 7 IICIP for some 7 > 0 and every ^ S t S K_|_, x € il, 

7i=i 

2) A similar result holds true for Neumann boundary conditions 


The following result is an immediate application of Theorem !10.3.1! using the 
Proposition !!0.4. l1 The smoothing effect of the heat equation implies (cf.||60| ) that 
for each e > 0 and a G [0,1), 

U{ u{t)} is bounded in 

t>e 

as soon as u{t) is bounded in for f > 0. In particular, Ut>o{^(^)} precom¬ 

pact in Lfp (fl) . 
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Theorem 10.4.3. Let u € Hq{D)) be a bounded solution of equation (110.281) . 

Assume that for all tp G S := {ip G H^ifl)/ — Ap + f{p) = 0} we have p G 
and one of the three conditions (110.301 ). (110.311 ) or (110.321) of Proposition \10.4.1\ is 
satisfied. Then 

lim ||M(i) - pWhi = 0. 

Moreover, let 9 be any Lojasiewicz exponent of E at p. Then we have 


\\u{f)-p\\L 2 


0(e for some 6>0if6 = ^, 

0(^-e/(i-2e)) if{)ce< i. 


Remark 10.4.4. It has been shown in lf55l that if d < 1, convergence holds without 
any need of condition (110.31b or (110.32b . However, if d = 1 and convergence occurs, 
in general the convergence can be arbitrarily slow. The hypothesis d < 1 provides 
convergence results in a wide framework, cf. e.g. iS, ED- 


10.4.2 A semilinear wave equation 


As a next application we study the asymptotic behaviour of the semilinear wave equa¬ 
tion 

{ utt + ut — Au + f{x, u) = 0, (f, x) G M+ X 

w(f,-)|an = 0, tGR+, (10.35) 

u(0,a;) = uo{x),UtiO,x) = ui(x), x G fl. 

We let fl C M'^, / G (7^(12 x ]R.;R), the spaces H := Lf{Tl) and V := H{(fl) as 
in Subsection 11 0.4. 11 If TV > 2, then we replace the growth condition (110.29b by the 
following condition; 


3(7 > 0, a > 0 such that (TV — 2)a < 2 
and s)| < (7(1 + |s|“) a.e. on H x R 

Theorem 10.4.5. Let u be a solution of (110.35b such that 


(10.36) 


U(>o{w(T,.), Ut{t, .)} is bounded in Hq{TI) x Lf{ri). 

Assume that for all p G S := {p G Hq{D)/ — Ap + f{p) = 0} we have p G L°°{D) 
and one of the three conditions (110.30b or (110.31b or (110.32b of ProDosition \Tb.4.1\ is 
satisfied. Then 

lim ||Mt|lL 2 -f ||M(f) - pWhi = 0. 

t—foo 

Moreover, let 9 be any Lojasiewicz exponent of E at p. Then we have 


IkW - pWl^ 


0{e for some S>0if9 = ^, 

(7(t-®/(i-2«)) if0<9< i. 


Proof. First (110.36b implies that the Nemytskii operator associated to / is compact: 
H{{fl) —!► Lf{Vt), then by the lemma 16.6.21 the orbit U(>g{u(T,.), Mt(f,.)} is pre¬ 
compact in H^ifl) X Lf{Vl). This is condition (i) of theorem 110.3.21 Moreover, 
the duality mapping K : H~^{Ll) —ddg(il) is given by Kv = (—A)“^u, so that 
KE"{v) = I + {—A)~^f'{v). From this, the growth assumption on / (110.36b . and 
the Sobolev embedding theorem, it is not difficult to deduce that the condition (ii) of 
Theorem ll0.3.21 is satisfied. □ 








































Chapter 11 

Variants and additional results 


In this chapter, we collect, most of the time without proofs a few additional results 
which complement, mainly in the inhnite dimensional framework and often at the price 
of additional technicalities, the simple theory developed in the two previous chapters. 
For the proofs, the reader is invited to read the corresponding specialized papers 

11.1 Convergence in natural norms 

In the last chapter, we obtained convergence to equilibrium for some semi-linear parabolic 
and hyperbolic equations in the energy space. However the rate of convergence to 
equilibrium was specihed in In Il56l . it is shown that the same decay occurs in 

IFq (H) for the wave equation and in L°°{n) with an arbitrarily small loss for the heat 
equation. This loss is most probably artihcial but this becomes only important when 
the Lojasiewicz exponent of tp is exactly known, which is possible only in exceptional 
cases. 

11.2 Convergence without growth restriction for the heat 
equation 

In 1641, the second author gave a proof of the Simon convergence theorem (cf. m in 
the framework of Sobolev spaces instead of (7“ spaces which were used by L. Simon. 
His proof is quite similar to that of our main parabolic result, but uses more compli¬ 
cated spaces. The advantage is that no growth restriction is assumed for the nonlinear 
perturbative term. 


11.3 More general applications 

11.3.1 Systems 

Let V = (iTQ(H))"', H = (L^(H))”, V' = and we define the function 


E : (iTo^(H))" —^ Rby 
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When TV > 2, we assume that 

II s)|| < (7(1 + ||s||“) a.e. on x M (H-l) 

for some (7 > 0 and a > 0 such that {N — 2)a < 2. By using Proposition 1.17.5 page 
66 of ll66ll . we know that E is of class on Hq (17) and 

DE{u) = {-Aui + fi{x,u),--- ,-Aun + fn{x,u)) 

D^Eium = <-Aa + ^(x,u)ei,--- ,-Ae„ + |^(a;,u)C„)VCG (i7o'(f^))"- 

OSi OSn 

It is well known that dimker D‘^E{ip) is hnite for all Lp G (i7g (17))" fl (L°°(17))". Let 
d — dimkerZ7i?((^). 

Proposition 11.3.1. Assume that hypothesis (II 1. Il l is satisfied. Let ip G {Hq{TI))'^ H 
(L°“(17))" be a critical point of E. Assume also that one of the following hypotheses 
is satisfied: 

d = 0 

d > 0 and there exists O open, and h G C^{0, V)/(p G h{0) C {DE)~^{0) 
and h : O —> h{0) is a diffeomorphism; 
f is analytic in s, uniformly with respect to x G El 

Then there exist 9 G (0, and a > 0 such that 

Vu G (i7o'(f^)r, lh-V>IUi(n) < ^ ^ \\DE{u)\\^H-Hn))’> > \E{u)-E{ip)\^-^ 

(11.2) 

11.3.2 Fourth order operators 

Let V = Hq{EI), H = L^(17), V' = and we dehne the function E : 

77^(17) —^ Rby 

VuGHq{E1.), E{u) = ^ I \Auf dx + I F{u)dx 
2 Jo. Jq. 


where F{u) = f{s) ds. When iV > 4, we assume that f{x, 0) G L°°{D) and 
df 

| —(a;,s)| < (7(1 + |s|“) a.e. on 17 x K (1L3) 

as 

for some (7 > 0 and a > 0 such that {N — 4)(a + 1) < + 4. By using Proposition 

1.17.5 page 66 of |[66l, we know that E is of class (7^ on TLq (17) and 

< DE{u),ip >h~ 2 ^h^ = < A‘^u +f{x,u),'ip G 

df 

< D^E{u)^,fi >h-2^h^ = < A2^+^(x.u)^,?)’>ff-2xH2 VV^ G i7o(17). 

It is well known that dimker Li'((/?) is finite for all tp G Hq{D). Letd = dim keri?'(</?). 
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Proposition 11.3.2. Assume that hypothesis (II 1.31 l is satisfied. Let ip G Hq(I 2) n 
be a critical point of E. Assume also that one of the following hypotheses is 
satisfied: 

d = 0 (11.4) 

d > 0 and there exists O open, and h S C^{0, V)/ (11-5) 

p £ h{0) C {DE)~^ (0) and h : O —> h{0) is a dijfeomorphism; 
f is analytic in s, uniformly with respect to x £ El (11-6) 

Then there exist 9 £ (0, and a > 0 such that 

Vm G H^{Et), ||M-<p||ff=(n) < cr ^ \\A^u+f{x,u)\\H-2(n) > \E{u)-E{p)\^-'^ 

(11.7) 

Remark 11.3.3. In virtue of remark [TO.0.41 if p is not a critical point of E, (110.331 1 
is just the consequence of the fact that E £ C^(V, V). In this case we don’t have to 
assume any assumption. 

Proof. The proof of (II 1.71) under hypotheses (1 11 .41 ) and dl 1.51) is the same as in the 
proposition 110.4.1] Now assume that (II 1.61) holds. As in the proof of the proposition 
110.4.11 we distinguish two cases ; 

Case 1 : iV < 3. Let Z = L^(n), by elliptic regularity ||4l, we know that W := (11 + 
A)~^(Z) C H'^(fl) where 11 is the orthogonal projection in L'^[El) on N := ker A. It 
is also clear that iV £L Z = L'^{Et). The functional : Hq{EI) —> K is clearly analytic 
since it is the sum of a continuous quadratic functional and a Nemytskii operator which 
is analytic on the Banach algebra 77^(17) C L°°{D) (see Example 1 1.4. 71 .1 By using 
Proposition 1 1.4. 51 1. we also obtain that DE : W —> Z is analytic. We can apply 
Theorem ll0.2.7l to obtain (110.331 ). 

Case 2 : > 4. Letp > max(2, and Z = LP(EI). By elliptic regularity [|4], we 

know that W :=(!! + A)~^(Z) C lC^’^(r7) which is a Banach algebra since p > ^. 
The end is the same as in the first case. □ 


11.4 The wave equation with nonlinear damping 

In ll2^ . L. Chergui succeeded to generalize Theorem l9.3.2l to the semilinear wave equa¬ 
tion with nonlinear localized damping 

{ Utt + \ut\°'ut - Au + f{x,u) = 0, {t,x) G K+ X 17, 

•)lan=0, 7 gK+, (11.8) 

u{0,x) = uo{x),Ut{0,x) = ui{x), X £ El. 

One of the difficulties to do that is the proof of compactness of the trajectories in the 
energy space. His result has been extended, under natural hypotheses, to possibly 
nonlocal damping terms in El. 

11.5 Some explicit decay rates under additional condi¬ 
tions 

The Lojasiewicz exponent of an equilibrium point is generally difficult to find, even for 
2-dimensional ODE systems. However in some exceptional case, it turns out, for semi- 
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linear problems involving a power non-linearity, to be computable explicitely. This 
was done in |59l with application to the exact decay of the solution when the limit is 
0, and in 03 under a positivity condition of the energy. The last result allows for a 
continuum of equilibria to exist, but only for Neuman boundary conditions. 


11.6 More information about decay rates 

All our convergence results contain an estimate of the difference between the limiting 
equilibrium and the solution. The question naturally arises of the optimality of this 
estimate. Actually, even when the equation has a single equilibrium playing the role 
of a universal attractor of all solutions, the situation can be rather complicated. If 
the decay estimate obtained for instance by Liapunov’s direct method or Lojasiewicz 
method is optimal for all solutions other than the rest point itself, it means that all non¬ 
trivial solutions tend to the quilibrium at the same rate, a circumstance which tends 
to be the exception rather than the general rule. As an illustration, let us consider the 
simple ODE 


u" + u' + = 0 

Apart from the zero solution, it is true (although not completely trivial to prove, cf. 
e.g.ill that here are only two possible rates of decay: as f 2 or as e *. Actually the 
first case corresponds to solutions behaving as those of u" -f = 0 and is shared by 
most solutions, while the ranges of exponentially decaying solutions lie on a separatrix 
made of two curves symmetric with respect to the origin(0, 0) having the rough shape 
of spirals. 

Analogous properties have been established by the first author for the slightly more 
complicated equation u"-|-c|u'I I It I = 0 wherec,a,/3 are positive constants. If 

a > ao := all trajectories are oscillatory up to infinity and tend to 0 at the same 
rate. If a < ao, all trajectories have a finite number of zeroes on [0, 00 ) and there are 
two different rates of decay at infinity . For the details , cf m- 

In a series of papers, the exact decay rate of solutions have been thoroughly studied 
for more complicated second order ODE and for infinite-dimensional abstract problems 
containing semilinear parabolic and hyperbolic equations. We refer to lfT3lfT4l[JTll38l 
1^ for the details. 


11.7 The asymptotically autonomous case 

It is natural to ask whether the convergence results are robust under a perturbative 
source which dies off sufficiently fast for t large. Such results were obtained in ll62ll . 

EH, gqi, insi and im. 

11.8 Non convergence for heat and wave equations 

Non convergence results for parabolic and hyperbolic equations with smooth non- 
analytic nonlinearities were proved by lEa, GSl and |j65l. Although such negative 
results may look natural since 2 dimensional ODE systems already produce such bad 
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phenomena, the question is whether or not the fact that the generating function is scalar 
forces the system to behave like a scalar equation. The answer is negative but the proof 
is non-trivial. 
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